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I.  INTRODUCTION 

Most  vulnerability  lethality  analyses  completed  at  the  Ballistic 
Research  Laboratory  (BRL)  today  require  some  type  of  geometric  description 
of  the  target.  This  geometric  description  is  usually  completed  employing 
the  combinatorial  geometry  (COM-GEOM)  technique.  The  COM-GEOM  technique 
of  target  description  represents  components  of  a target  as  Boolean  combin- 
ations of  twelve  basic  geometric  solids.  These  solids  are  listed  in  Table 
I. 


As  part  of  a continuing  effort  to  increase  the  accuracy  and  reduce 
the  production  time  of  these  COM-GEOM  descriptions,  a set  of  geometry  aid 
programs  has  been  developed.  The  purpose  of  these  aid  programs  is  to 
quickly  give  the  COM-GEOMer  solutions  to  the  bulk  of  geometry  problems 
that  arise  during  the  production  of  a target  description. 

Since  quick  turnaround  time  is  imperative,  these  aid  programs  have 
been  coded  for  the  WANG  2200  mini-computer  system.  Many  of  these  programs 
are  also  coded  and  available  for  the  WANG  700  series  calculator.  The 
BASIC  type  language  used  on  the  WANG  2200  system  can  be  converted  to  other 
systems,  if  necessary. 

This  report  documents  the  aid  programs  which  are  presently  available 
for  the  WANG  2200  mini-computer  system. 


II.  BACKGROUND 

In  the  course  of  producing  COM-GEOM  target  descriptions,  one  soon 
becomes  aware  of  the  urgent  need  for  solutions  to  relatively  simple 
geometric  problems.  These  include  such  problems  as  vectors  normal  to 
other  vectors;  four  points  in  a plane;  the  intersection  point  of  two 
lines,  etc.  The  requirement  for  quick,  accurate  answers  to  such  problems 
led  to  the  development  of  this  set  of  mini-computer  codes.  Indexed  in 
Table  II  are  the  names  and  functions  of  each  code.  These  available  codes 
will  be  presented  in  a user  oriented  manner.  For  each  of  the  geometry 
aid  programs,  the  following  are  included: 

(1)  Descriptive  paragraphs  on  the  function  and  use 

(2)  Memory  requirement 

(3)  Any  restrictions,  limitations,  or  special  features 

(4)  Instructions  for  use  on  the  WANG  2200  system 

(5)  Sample  output  listings  (note:  whenever  feasible  at  least  one 
of  the  sample  outputs  will  be  a case  where  the  answers  can  be 
readily  checked) 

Appendix  A contains  listings  of  The  codes  for  each  program.  All 
these  codes  are  available  at  the  BRL  on  cassette  magnetic  tape. 


Table  I.  Geometric  Solids  Used  In  COM-GEOM  Descriptions 


SYMBOL 

SOLID  NAME 

RPP 

Rectangular  Parallepiped 

BOX 

Box 

RAW 

Right  Angle  Wedge 

ARB 

Arbitrary  Convex  Polyhedron 

ARS 

Triangular  Surfaced  Polyhedron 

ELL 

Ellipsoid  of  Revolution 

SPH 

Sphere 

RCC 

Right  Circular  Cylinder 

REC 

Right  Elliptical  Cylinder 

TRC 

Truncated  Right  Angle  Cone 

TEC 

Truncated  Elliptic  Cone 

TOR 

Torus 

Table  II.  Index  of  the  Geometry  Aid  Programs 


L 


t 

* 


f 


Number  Name 

1 SPHERE 

2 CIRCIR 

3 CIRCLE 

4 RCC 

5 LINECIR 

6 TANCIR 

7 PLANEINT 

8 LINEPLAN 

9 LINELINE 

10  RFARB 

11  3PTARB 

12  NORMVEC 

13  PERPENV 

14  AMT RACK 


Purpose 

Finds  a sphere  defined  by  4 nonplanar 
points  on  the  surface  of  the  sphere. 

Finds  the  intersection  points  of  2 
circles  (2  dimensions). 

Finds  a circle  defined  by  3 noncol linear 
points  on  the  circumference -of  the  circle 
(2  dimensions). 

Finds  RCC  defined  by  3 noncol linear 
points  on  circumference  of  the  base  and 
a desired  height.  Can  also  be  used  in 
same  manner  as  #3  but  in  3 dimensions. 

Finds  the  intersection  points  of  a line 
and  circle  (2  dimensions). 

Finds  tangent  points  on  a circle  from  a 
point  outside  the  circle  (2  dimensions). 

Finds  the  point  of  intersection  of  3 
different  planes. 

Finds  the  point  of  intersection  of  a line 
and  plane. 

Finds  intersection  point  of  2 lines 
(2  dimensions) . 

Finds  an  ARB8  defined  by  (1)  point  on  one 
face;  (2)  rotation  and  fallback  angles 
of  that  face;  (3)  2 coordinates  of  the 
remaining  3 points  on  that  face,  and  (4) 
a desired  thickness. 

Finds  an  ARB8  defined  by  (1)  3 noncollinear 
points  on  one  face;  (2)  2 coordinates  of 
the  remaining  point  on  that  face;  (3)  a 
desired  thickness. 

Finds  (1)  vector  of  desired  length  in 
the  direction  of  a given  vector  and  (2) 
a vector  of  desired  length  perpendicular 
to  the  given  vector  (2  dimensions). 

Finds  (1)  a vector  of  desired  length 
perpendicular  to  2 given  vectors  and  (2) 
the  angle  between  the  2 given  vectors. 

Adds  tracks  to  domestic  vehicles  - listing 
the  solids  and  a region  table. 
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Table  II.  Index  of  the  Geometry  Aid  Programs  (Continued) 


Number 

Name 

Purpose 

IS 

S0LIDR0T 

Rotates  selected  COM-GEOM  solids  about 
any  point  in  the  XY,  XZ,  or  YZ  planes. 

16 

PL0TS0L 

Plots  selected  COM-GEOM  solids  at 
desired  aspect.  No  hidden  line 
algorithm. 

any 

17 

D ARB IN 

Corrects  a "bad"  ARB  and  computes 
inside  one. 

an 

18 

BOXIN 

Corrects  a "bad"  BOX  and  computes 
inside  one. 

an 

19 

RANIN 

Corrects  a "bad"  RAW  and  computes 
inside  one. 

an 

20 

TRCIN 

Computes  an  inside  TRC. 

21 

RECIN 

Corrects  a "bad"  REC  and  computes 
inside  one. 

an 

22 

TEC  IN 

Corrects  a "bad"  TEC  and  computes 
inside  one. 

an 

23 

PARB 

Computes  the  points  of  intersect  of  a 
set  of  planes  taken  3 at  a time. 
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All  of  these  programs  have  been  coded  to  Interact  with  the  user  via 
the  scope  of  the  WANG  2200  system.  Each  input  is  asked  for  in  a step-wise 
manner  and  all  instructions  are  listed  on  the  screen.  The  loading 
instructions  for  these  programs  are  all  the  same. 

To  load  any  of  these  programs  from  cassette  tape: 

(1)  Place  the  tape  in  the  machine  and  rewind 

(2)  Key  LOAD  ’’Program  Name"  CR/LF 

(3)  After  program  is  loaded  in  machine,  key  RUN  CR/LF. 

At  this  point  instructions  and  the  request  for  the  first  inputs  will  be 
displayed  on  the  screen.  For  specific  directions  on  the  input  required 
see  the  instruction  section  of  the  individual  program.  Wnen  finished 
running  any  of  the  programs  perform  the  following  steps: 

(1)  Key  RESET 

(2)  Key  CLEAR  CR/LF 

(3)  Rewind  the  tape  and  remove  it  from  the  machine. 

Even  though  some  of  these  programs  operate  on  two-dimensional 
problems  only,  they  are  useful  because  many  three-dimensional  space 
problems  can  be  reduced  to  the  two-dimensional  case.  For  example,  most 
of  the  boxes  (BOX)  used  in  COM-GEOM  descriptions  employ  vectors  with  only 
2 nonzero  components  and  many  right  circular  cylinders  (RCC)  have  height 
vectors  of  2 (and  often  only  1)  nonzero  components. 


III.  DISCUSSION  OF  THE  PROGRAMS 


Program  Name:  SPHERE 

Description:  This  program  finds  the  sphere  defined  by  any 

four  noncoplanar  points  on  the  surface  of  the 
sphere. 


Memory:  1467  bytes 

Restrictions:  None 

Instructions:  After  loading  the  program,  enter  the  following: 

(enter  data  separated  by  commas) 

(a)  X,Y,Z  coordinates  of  point  1 on  the  surface 
CR/LF 

(b)  X,Y,Z  coordinates  of  point  2 on  the  surface 
CR/LF 

(c)  X,Y,Z  coordinates  of  point  ? >n  the  surface 
CR/LF 

(d)  X,Y,Z  coordinates  of  point  4 on  the  surface 
CR/LF 

(e)  The  results  will  be  printed  out 

(f)  To  run  again  key  CONTINUE  CR/LF  and  program 
returns  to  step  (a). 

Sample  Outputs:  First  sample  run  is  case  which  is  readily  checked. 


INPUT  POINTS 
f'T  X 

Y 

z 

1 

5.0000 

3.0000 

0.0000 

■» 

0.0000 

5.0000 

0.0000 

o 

0.0000 

0.0000 

y.0<)00 

o.nooo 

0.0000 

-5.0000 

SPHERE 

xc- 

J.0G00  YC- 

0.0000 

RADIUS" 

5.0000 

0.0*on 


INPUT  POINTS 
PT  X 

1 203.2300 

2 123.5402 

3 1.0000 

4 8.0000 

SPHERE  XO 

RADIUS 


Y 

65.2340 
321.2540 
2.0000 
6.0000 
3727.2686  YC- 
- 3947.7220 


Z 

-1025.3260 

321.0214 

3.0000 

1.0000 

-1289.6832  ZC-  ' 


270.7310 


. Program  Name:  CIRCIR 

Description:  This  program  finds  the  intersection  points  of  2 

circles  (in  2 dimensions).  This  program  is 
useful  in  checking  for  the  overlapping  of  stored 
missiles  or  ammunition. 


Memory:  1502  bytes 

Restrictions:  None 

Instructions:  After  loading  the  program,  enter  the  following: 

(enter  data  separated  by  comas) 

(a)  The  X,Y  coordinates  of  the  center  and  the 
raditis  of  circle  1 CR/LF 

(b)  The  X,Y  coordinates  of  the  center  and  the 
radius  of  circle  2 CR/LF 

(c)  Results  will  be  printed  out 

(d)  To  run  again  key  CONTINUE  CR/LF  and  program 
returns  to  step  (a). 


Sample  Outputs:  First  sample  output  is  a case  where  circles  are 
tangent.  The  second  sample  run  is  a case  where 
one  circle  is  inside  the  other.  Note  that  the 
program  gives  a "no  intersection"  result. 


CIRCLE  1 X CENT- 
CIRCLE  2 

CIRCLES  TAflGEfiT  AT 


0.0000  Y CENT- 
15.0000 

X-  5.0000  Y- 


0.0000  R»  5.0000 

0.0000  10."000 
0.0000 


CIRCLE  1 X CENT-  0.0000  Y CENT- 
CIRCLE  2 0.0000 

NO  INTERSECTION 


0.0000  R-  5.0000 

0.0000  10.0000 


CIRCLE  1 X CENT- 
CIRCLE  2 

INTERSECTION  POINTS 


12.3564  Y CENT-  -36.6523  ,R» 

15.2360  39.2365 

X-  57,3018  Y-  1.0763 
-29.6190  5.1776 


59.2600 

56.3210 


3.  Program  Name:  CIRCLE 

Description:  This  program  finds  the  circle  defined  by  3 

noncollinear  points  on  the  circumference  of  the 
circle  (in  2 dimensions).  If  the  3 dimensional 
case  is  required  use  program  4 (RCC). 


Memory:  1093  bytes 

Restrictions:  None 

Instructions:  After  loading  the  program,  enter  the  following: 

(enter  data  separated  by  conmas) 

(a)  The  2 coordinates  of  point  1 CR/LF 

(b)  The  2 coordinates  of  point  2 CR/LF 

(c)  The  2 coordinates  of  point  3 CR/LF 

(d)  Results  will  be  printed  out 

(e)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  returns  to  step  (a). 

Sample  Outputs: 


INPUT  POINTS 
X Y 

5.0000  0.0000 

0.0000  5.0000 

0.0000  -5.0000 

XCENT»  0.0000  YCENT-  0.0000  RAO-  5.0000 


INPUT  POINTS 
X Y 

125.3050  G1 3.2483 

-561.2357  0512.3570 

-00.2350  251.2470 

XCENT-  -4207.7907  YCENT-  3000.1008  RAO-  5060.2272 


4. 


Program  Name: 
Description: 


RCC 


This  program  finds  an  RCC  defined  by:  (1)  3 
noncol linear  points  on  the  circumference  of  the 
base;  (2)  length  of  height  vector  desired. 


Memory: 

Restrictions: 


Instructions : 


3646  bytes 

The  user  must  beware  of  the  direction  of  the 
height  vector.  The  program  will  print  out  the 
two  possible  height  vectors  and  the  user  must 
choose  the  one  he  desires.  Note  that  this 
program  can  also  be  used  to  fit  a circle  (in 
three  dimensions)  through  3 noncol linear  points. 
After  loading  the  program,  enter  the  following: 
(enter  data  separated  by  commas) 

(a)  Key  X,Y,Z  coordinates  of  point  1 on  the 
circumference  CR/LF 
Key  X,Y,Z  coordinates  of  point  2 on  the 
circumference  CR/LF 
Key  X,Y,Z  coordinates  of  point  3 on  the 
circumference  CR/LF 
Key  the  length  of  height  vector  desired 
Results  will  be  printed  out 
Program  returns  to  step  (a). 

Sample  Outputs:  Note  that  if  this  program  is  used  to  fit  a circle 
through  3 points,  use  the  circle  defined  by  the 
center  and  radius  of  the  base  of  the  RCC. 


(b) 

(c) 


(d) 

(e) 

(f) 


CR/LF 


nniT  poi.rrs: 

POIfIT  1 X- 

10.0000 

Y» 

0.0000 

z» 

o.oooo 

POIilT  2 X- 

o.nooo 

Y» 

10.0000 

7. 

o.oooo 

ponrr  o x- 

-10.0000 

Y- 

o.oooo 

> 

0.0000 

LETIOTM  OF  MFIOilT  VfCTOP- 

no 

T::n  papa  inns  of  run  ncc: 
cr::Tm  of  msf.  xc»  o.oooo  yc« 
nx-  0.0000  nv» 

nx-  u.oooc  0Y» 

10. 0000 


I Tim  :T  7IXT0U 

(||t 

.ttilllif.  OF  MSF 


0.0000  zc» 

0.0000  PZ» 

n.'ww  r>7» 


0.0000 
-r, o.owo 
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Jit  7?-  - 


INPUT  POINTS: 


POINT 

1 

X® 

542.3205 

Y» 

321.1250 

7- 

8235.3750 

POINT 

2 

X- 

530.3241 

Y« 

052.1247 

Z» 

-74.2541 

POINT 

a 

X- 

821.0214 

Y- 

-052.3214 

Z» 

302.3714 

LENGTH  OF  HEIGHT  VECTOR-  100 

THE  PARAMETERS  OF  THE  RCC: 

CENTER  OF  BASE  XC-  101.9033  YC-  3000.4030  7.C-  4103.0340 
HEIGHT  VECTOR  DX-  9B.2502  DY-  10.502207-  1.2004 

OR  DX-  -00.2502  OY-  -10.0022  07-  -1.3004 

RADIUS  OF  RASE®  4000.2540 


>. 


Program  Name: 
Description: 


LINECIR 

This  program  finds  the  intersection  points  of  a 
line  and  a circle  (in  2 dimensions).  The  line 
is  entered  using  the  point  slope  method. 


Memory: 

Restrictions: 

Instructions: 


Sample  Outputs: 


1559  bytes 
None 

After  the  program  is  loaded,  enter  the  following: 
(enter  the  data  separated  by  commas) 

(a)  X and  Y coordinates,  DELTA  X,  DELTA  Y, 
for  the  line  CR/LF 

X,Y  coordinates  of  center  and  radius  of 
circle  CR/LF 

Results  will  then  be  printed  out 
To  run  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 


(b) 


(c) 

(d) 


line 

INPUT 

0.0000 

Y- 

0. 

0000 

del  :i- 

1,0000 

0EL 

Y- 

1. 

0000 

CIPXLE  INPUT 

xc- 

0.0000 

YC- 

0. 

oooo 

R- 

i:te;j 

ir.rcTiori 

POINTS 

» 

0.7071 

Y» 

0.7071 

-0.7071 

Y» 

-0.7071 

Lli.T 

INPUT 

X- 

0.0000 

Y» 

10, 

,0000 

nr.L  x- 

0.0000 

on 

Y» 

11, 

,0*00. 

Cir.flF.  INPUT 

xc« 

0.00)0 

YC« 

1. 

.0000 

0* 

l.onnn 


l.onpn 


•ID  INTERSECTION 
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LINE  INPUT 


o.ooon 


DEL  X- 

o.oooo  r 

CIRCLE  INPUT 

XC- 

o.nnoo  yc« 

INTERSECTION  POINTS 

X*  0.0000 

y« 

5.0000 

X«  0.0000 

Y» 

-5.0000 

LINE  INPUT 

X- 

5.0000 

OEL  X- 

0.0000 

CIRCLE  INPUT 

XC- 

o.nnno  vc 

LINE  TANGENT  TO  CIPXLE  AT  X- 

y«  o.ooon 

Y»  n.nnno 

o.nnno  R» 


y»  i n.nnno 

del  y»  li.onnn 

o.ooon  r- 
5. noon  y- 


5.  won 


i.onnn 


n.onoo 


LIME  INPUT  X- 
DEL  X- 

CIRCLE  INPUT  XC- 
INTERSECTION  POINTS 
X-  11)34.5473  Y- 
X»  -545.2092  Y- 


1203.2554 

1257.2175 

235.5540 

3256.2470 

2791.8301 


Y-  -2301.2356 
DEL  Y«  -3691.2547 
YC-  -326.1250  R- 


3214.2560 
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6.  Program  Name:  TANCIR 

Description:  This  program  finds  the  tangent  points  on  a circle 

from  a point  outside  the  circle.  (2  dimensions) 


Memory:  1153  bytes 

Restrictions:  None 

Instructions:  After  the  program  is  loaded,  enter  the  following: 

(enter  data  separated  by  commas) 

(a)  X and  Y coordinates  of  the  point  outside 
the  circle  CR/LF 

(b)  X and  Y coordinates  of  the  center  of  the 
circle  and  the  radius  of  the  circle  CR/LF 

(c)  Results  will  be  printed  out 

(d)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  returns  to  step  (a) . 

Sample  Outputs: 


POINT  IilPtn 

X* 

0.0000  Y« 

o.ooon 

CIRCLE  INPUT 

XCEI.'T* 

n.oooo 

YCL.'IT*  n.onnn  R» 

2.0000 

TA!inr.i;T  pts 

X- 

5. noon  y« 

0.0000 

> 

3.G20G  Y* 

3.M82 

P)I:;t  i::p:t 

> fio.oono  y» 

0.0000 

CIRCLE  INPUT 

:<ce:it= 

0.0000 

YCEflT3  0.0000  R« 

irv-ooo 

TA!!fiL,’,T  PTS 

x« 

n.oooo  y= 

9.7079 

X- 

n.oooo  y« 

-0.7979 
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poi:.t  iiiput  x-  nG.ocoo  y=  q.oooo 

CIRCLE  INPUT  XCENT*  0.0000  YCEIIT*  0.0000  R=  35. 
POINT  IS  ON  CIRCLE 


POINT  INPUT  X-  0.0000  Y=  0.0000 

CIRCLE  INPUT  XCENT*  0.0000  YCEIIT*  0.0000  R=  5C. 

POINT  IS  INSIDE  CIRCLE 


POINT  IIIPUT 

X*  1234. 5G70 

Y« 

2345. G709 

CIRCLE  IIIPUT 

XCENT*  120. 

,3G50 

YCEIIT* 

TANGENT  PTS 

X-  410,7042 

Y* 

33G4.0G44 

X-  -1.0377 

Y= 

27G2.707P 

POINT  INPUT 

X-  102.3040 

Y« 

-320.2509 

CIRCLE  IIIPUT 

XCENT*  0.0000 

YCEMT- 

TANGENT  PTS 

X*  11.RG07 

Y« 

3.257G 

X-  -11.5950 

Y* 

-4.1019 

3125.0140 


0.0000 


n» 


R* 


302 


12. 


OOOO 


oooo 


.3000 


3000 


7.  Program  Name:  PLANE INT 

Description:  This  program  finds  the  point  of  intersection  of 

three  different  planes.  Each  plane  can  be 
entered  using  one  of  three  methods:  (1)  3 
noncol linear  points;  (2)  one  point  on  the  plane 
and  rotation  (azimuth)  and  fallback  (elevation) 
angles  of  that  plane;  or  (3)  coefficients  of 
the  equation  of  the  plane  - that  is  A,B,C,D 
of  the  equation  AX  + BY  + CZ  = D. 


Memory: 

Restrictions: 


Instructions: 


Sample  Outputs: 


3632  bytes 

If  the  three  planes  intersect  in  a line  the 
program  will  print  out  the  message  that  the 
planes  do  not  intersect. 

After  the  program  is  leaded,  enter  the  following: 

(a)  An  indicator  depicting  the  method  of  input 
for  the  first  plane  where: 

Indicator  Type  of  Plane  Input 

1 3 noncol linear  points 

2 Point;  rotation, 

fallback  angles 

3 Coefficients  of  the 

plane 

(b)  Enter  plane  1 by  method  chosen 

(c)  Next,  enter  the  other  2 planes  using  steps 
(a)  and  (b)  above 

(d)  Results  will  then  be  printed  out 

(e)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 

The  first  sample  run  is  a case  easily  checked. 

In  the  second  sample  run  all  three  methods  of 
entering  planes  was  used.  Plane  1 was  entered 
by  the  3 point  method,  plane  2 by  the  point 
and  rotation  fallback  method,  and  plane  3 by 
coefficients  of  the  plane. 
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« 


IIIPIIT  FOR  PLANE  1 


COEFFICIENTS  A- 

1.00000  n- 

0.00000  CB 

o.onnno  d- 

99.00000 

INPUT  FOR  PLANE  2 
COEFFICIENTS  A= 

0.00000  B« 

1.00000  c» 

O.HQOOO  D= 

99.00000 

INPUT  FOR  PLANE  3 
COEFFICIENTS  A- 

0.00000  ft- 

0.00000  c- 

1.00900  n- 

99.00000 

INTERSECTION  POINT 

X«  99.0000 

Y-  99.0000 

Z-  99.0000 

INPUT  FOR  PLANE  1 

Xs  1203.3230  Y- 

922.2570  7- 

320.0000  ROT- 

04.933  F B- 

61.672 

1003.2305 

930.2451 

325.2360 

1145.2300 

707.0169 

390.2360 

COEFFICIENTS  A-  0.04190  B-  0.47203  C-  0.00025 

D»  770.5159 

INPUT  FOR  PLANE  2 

X*  1025.3090  Y= 

012.3697  Z- 

546.2387  ROT- 

i 4?. non  F B» 

20.340 

COEFFICIENTS  A=  0.09071  B-  0.02732  C-  0.34791 

0-  1414.0592 

INPUT  FOR  PLANE  3 
COEFFICIENTS  A- 

12.36900  R- 

-3.26540  C« 

210.36910  D» 

56.23409 

INTERSECTION  POiriT  X-  597.5406  Y*  1596.0051  Z«  -10.9990 


v 


lb 
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Program  Name:  LINEPLAN 

Description:  This  program  finds  the  intersection  point  of  a 

line  and  a plane.  The  line  is  entered  using 
the  point-slope  method.  The  plane  is  entered 
one  of  three  ways:  (1)  3 noncol linear  points; 

(2)  point  in  the  plane  plus  rotation  (azimuth) 
angle  and  fallback  (elevation)  angle  of  the 
plane;  or  (3)  coefficients  A,B,C,D  of  the 
equation  of  the  plane  AX  + BY  + CZ  = D. 


Memory : 
Restrict  ions: 
Instructions: 


Sample  uutputs: 


3360  bytes 
None 

After  the  program  is  loaded,  enter  the  following: 
(enter  data  separated  by  commas) 

(a)  X,Y,Z  coordinates  of  a point  on  the  line 
CR/LF 

(b)  DELTA  X,  DELTA  Y,  DELTA  Z of  the  line  CR/LF 

(c)  Indicator  for  type  of  plane  input  desired, 
where : 

Indicator  Type  of  Plane  Input 

1 3 noncollinear  points 

2 Point,  rotation  angle, 

fallback  angle 

3 Coefficients  of  the 

plane  equation 

(d)  Enter  plane  according  to  the  method  chosen 

(e)  Results  will  be  printed  out 

(f)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 

The  first  three  sample  outputs  are  the  same 
case  but  the  plane  has  been  entered  using  the 
three  different  methods. 
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■ +.-*■  A* 


INPUT  FOR  PLANE 
COEFFICIENTS  A3 

1.00000  B» 

0.00000  C3  0.00000  D3 

123.00000 

LINE  INPUT 
X3  0.0000  Y3 

DEL  X3  2.0000 

0.0000  Z3 
DEL  Y3  2.0000 

0,0000 

DEL  7-  2.0000 

INTERSECTION  POINT 

X3  100.0000 

Y*  ion. oono  Z3  loo.oooo 

INPUT  FOR  PLANE 
X*  100.0030  Y3 
COEFFICIENTS  A3  1 

2.0000  Z»  3.3000  ROT3  0.220  F B3 

.oooon  B3  o.onooo  c=  o.noooo  n-  loo. oono 

0.030 

LINE  INPUT 
X=  0.0000  Y3 

DEL  X3  2.0000 

o.nooo  Z3 

DEL  Y*  2.0000 

0.0000 

DEL  Z3  2.0000 

INTERSECTION  POINT 

X3  ion. ooon 

Y3  100. 0000  7-  100.0022 

INPUT  FOR  PLANE 
X=  100.0000  Y« 
100.0000 
122.0000 

2.0000  Z3 

5.0000 

6.0000 

3.0000  ROT-  0.000  F 5* 

6. oono 

9.0000 

0.230 

COEFFICIEIITS  A*  1.00000  B3  0.00000  C*  0.00000  0=  100.0000 

LiriE  INPUT 

x=  o.ouoo  y-  o.oooo  z=  o.nooo 

PEL  » 0.0000  DEL  Y3  2.D000  DEI.  Z"  2.0D00 


INTERSECTION  POINT 


100.0C00  Y»  100.0000  Z-  100.0000 


< 


f 

K 


» 1 


INPUT  FOR  PLANE 

X*  23.5600  Y»  370.2360  Z-  -657.2350  ROT*  331.100  F R«  -25.200 
120.3245  -23.1500  -65.242.0 

02.1547  -3256.2140  3201.2450 

COEFFICIENTS  A*  0.79202  U--0.43724  C— 0.42593  D-  133.32T. 


LINE  INPUT 

:•>  125.6020  Y=  256.3240  Z*  -235.652.0 

PEL  X3  -2364.0123  PEL  Y3  352, 9.2345  DEL  7. 


INTERSECTION  POINT 


56.2342 

156.5755  Y»  209.2440  7-  -236.  poor, 


K 

« 

t 
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Program  Name : LINELINE 

Description:  This  program  finds  the  point  of  intersection 

of  two  lines  (in  2 dimensions).  The  lines  are 
entered  using  the  point-slope  method. 


Memory:  1108  bytes 

Restrictions:  None 

Instructions:  After  the  program  is  loaded,  enter  the 

following: 

(enter  the  data  separated  by  commas) 

(a)  X and  Y coordinates  of  a point  on  the 
line,  DELTA  X,  DELTA  Y for  line  1 CR/LF 

(b)  X and  Y coordinates  of  a point  on  the 
line,  DELTA  X,  DELTA  Y for  line  2 CR/LF 

(c)  Results  will  be  printed  out 

(d)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 

Sample  Outputs: 


LINE  X 

1 o.ooon 

2 2,0000 
INTERSECTION  POINT 


Y DFL  X DFL  Y 

0.0000  1.0000  l.ooon 

0.0000  0.0000  11 .0000 

X*  2.0000  Y«  2.0000 


LINE  X Y 

’ i 0.0000  0.0000 

2 2.0000  0.0000 

LINES  DO  NOT  INTERS  XT 


DEL  X DEL  Y 

1.0000  1.0000 

l.WVW  1.0000 


LINE  X 

1 120,3692 

? 104.2350 

INTERSECTION  POINT 


Y DEL  X DFL  Y 

-23.5690  1235.0470  -321.2014 

364.25.00  -23.5507  1452.3000 

X*  110.4D71  Y*  -31.0010 
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10.  Program  Name:  RFARB 

Description:  This  program  finds  an  ARB8  with  face  1234  parallel 

to  face  5678  and  a normal  distance  of  any  desired 
length  (T)  apart  (see  figure  below) . Face  1234 
is  defined  by  the  following: 

(1)  one  point  on  the  face 

(2)  two  coordinates  of  the  remaining  three 
points  on  that  face 

(3)  rotation  (azimuth)  angle  for  that  face 

(4)  fallback  (elevation)  angle  for  that  face 


Memory:  2006  bytes 

Restrictions:  Note  that  ARB6s  can  also  be  done  by  making 

point  4 equal  to  point  1.  The  user  must 
beware  that  face  5678  is  on  the  desired  side 
of  face  1234.  If  not,  run  the  program  again 
using  a negative  thickness  T (see  sample 
outputs  1 and  2). 

Instructions:  After  the  program  has  been  loaded,  enter  the 

following: 

(enter  data  separated  by  commas) 

(a)  The  X,Y,Z  coordinates  of  one  point  on  face 
1234  and  the  rotation  and  fallback  angles 
of  face  1234  CR/LF 

(b)  Two  coordinates  of  the  remaining  3 points 
will  then  be  entered  in  the  following 
manner: 

For  each  point,  you  will  be  asked  to 
enter  indicators  and  2 known  coordinates. 
If  the  X and  Y coordinates  are  known  then 
enter  0,0, 1,X  coordinate,  Y coordinate. 


f 

Follow  the  same  scheme 

for  other  coordinates 

wu 

So  enter  the  remaining 

3 points  as 

1,0, 0,Y  coordinate. 

Z coordinate 

CR/LF 

». 

or  0,1, 0,X  coordinate. 

Z coordinate 

CR/LF 

1 

% 

or  0,0, 1,X  coordinate. 

Y coordinate 

CR/LF 

% 

i 


(c)  The  distance  between  face  1234  and  face 
5678  (this  variable  is  called  thickness) 
CR/LF 

(d)  Results  will  be  printed  out 

(e)  To  run  another  case,  key  CONTINUE  CR/LF 
and  prograa  will  return  to  step  (a). 

At  this  point  it  is  usually  advisable 
to  run  another  case  with  a thickness 
which  is  the  negative  of  the  previous 
case's  thickness. 


Sample  Outputs: 


1 10(53.0000 

2 1024.G2G5 

3 109O.4G40 

4 114G.4420 

5 10G7.9GGG 

6 1029.5930 

7 1095.4305 

8 1151.4035 

FACE  A 
1234  0.59332 

5078  9.09832 

THICKNESS- 


Y 

922.2570 

359.0320 

707.1G2G 

787.1570 

925.7288 

3G3.5038 

791.5345 

791.6288 


320.0000 

514.5340 

512.2850 

393.0010 

322.4324 

517.0154 

514.7174 

400.4334 


0.52877  0.34202  1431.66237 
0.62877  0.34202  1433.77487 
7.112 


ROT 

42.00 

42.00 


F B 
20. 0C 
20.00 


j 

PT 

X 

Y 

z • 

. 

1 

1UG3.0000 

922.2570 

320.0000 

[ 

2 

1024.62G5 

859.0320 

514.5840 

J 

3 

1090.4G40 

787.1 62G 

512.2850 

f 

4 

1146.4420 

737.1570 

398.0010 

f 

!) 

1058.0334 

917.7851 

317.5675 

1 

0 

1019.6600 

854.5601 

512.1515 

I 

7 

1035.4974 

732.6908 

509.8525 

1 

c 

«# 

1141.4754 

702.6851 

395.5685 

' 

FACE 

• * 

i i 

n c 

1) 

ROT 

F B 

' 

1234 

0»6‘.<%*32 

J. 52877  0.3420E 

I 1431.66287 

42.00 

20.00 

* 

5678 

0.69832 

0.52877  0.3420E 

! 1424.55087 

42.00 

20.0'' 

- 

THICKNESS- 

■7.112 

\ * 

?T 

u 

A 

Y 

Z 

f *■  * 

1 

123.2650 

23.1200 

55.0000 

h • 

n 

t. 

331.2550 

-G523.1208 

55.0000 

3 

56.2358 

-523.2158 

55.0000 

h 

4 

5.  ;2. 3590 

-453.2147 

55.0000 

i 

5 

123.2650 

23.1200 

75.0000 

|1  ’ 

it 

231.2560 

-5523.1 20" 

75.0000 

r 

7 

56.2558 

-523.2158 

75.0000 

* 

• I 

:»92.5890 

-453.2147 

75.0000 

a 

t 

r\cr 

i* 

/• 

c 

»> 

:r-T 

f 

L> 

12.V 

•1.0000;; 

i-.rnooo  1.00300  ss.'^ono 

n.’oo 

oo 

• 

Ilvs 

Pt 

3t  7 i 

*5.9* ‘600 

j.OWI  1.1  ‘0/ 

j 7r.o-r'nn 

.on 

043  'in 

• 

t:  ic:: 

2) 
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11.  Program  Name: 
Description: 


Memory: 

Restrictions: 


Instructions: 


3PTARB 

This  program  finds  an  ARB8  with  the  face  1234 
parallel  to  the  face  5678  and  a normal  distance 
of  any  desired  length  (T)  apart  (see  figure 
below).  Face  1234  is  defined  by  the  following: 

(1)  3 noncollinear  points  and 

(2)  2 coordinates  of  the  remaining  point 
This  program  is  useful  in  finding  armor 
plates  of  any  desired  thickness. 


2666  bytes 

ARB6s  can  be  done  by  making  the  4th  point 
equal  to  point  1.  The  user  must  beware  that 
5678  is  on  the  correct  side  of  face  1234.  If 
not,  rerun  the  program  using  a negative  T. 

After  the  program  has  been  loaded,  enter  the 
following: 

(enter  the  data  separated  by  commas) 

(a)  The  X,Y,Z  coordinates  of  point  1 CR/LF 

(b)  The  X,Y,Z  coordinates  of  point  2 CR/LF 

(c)  The  X,Y,Z  coordinates  of  point  3 CR/LF 

(d)  The  remaining  point  of  face  1234  in  the 

same  manner  as  the  RFARB  program,  that 
is: 

1,0, 0,Y  coordinate,  Z coordinate  CR/LF 
or  0,1, 0,X  coordinate,  Z coordinate  CR/LF 
or  0,0, 1,X  coordinate,  Y coordinate  CR/LF 

(e)  The  thickness  desired  (i.e.  the  distance 
between  face  1234  and  face  5678) 

(f)  Results  will  be  printed  out 

(g)  To  run  another  case,  key  CONTINUE  CR/LF 
and  the  program  will  return  to  step  (a). 

At  this  point,  it  is  advisable  to  run  the 
same  problem  only  with  the  thickness 
negative  to  the  thickness  of  the  previous 
run. 
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Sample  Outputs 


PT 

V 

A 

Y 

7 

/. 

1 

lon.i.onon 

922.2570 

320.0000 

O 

f. 

i024.G2:r» 

259.0320 

514.6840 

•) 

1090.4040 

787.1626 

512.2850 

4 

1140.4420 

787.1570 

398.0010 

5 

1007.9065 

926.7208 

322.4324 

G 

1029.5930 

863.5038 

517.0164 

7 

1095.4305 

791 .6344 

514.7174 

U 

1151.4005 

791.6228 

400.4334 

FACE 

A 

a C 

n 

ROT 

rn 

ir.4 

3.69332 

0.62877  0.34202  1431.(6281 

42.00 

20.00 

507!* 

0.09832 

0.62077  0.34202  1438.77481 

42.00 

20.00 

THICKIIESS- 

7.112 

PT 

X 

Y 

7. 

1 

99.0000 

123.0000 

456.0000 

2 

99.0000 

450.0000 

734.0000 

3 

99.0000 

569.0000 

912.0000 

4 

99.0000 

5632.0000 

8952.0000 

5 

109.0000 

123.0000 

456.0000 

n 

109.0000 

450.0000 

784.0000 

7 

109.0000 

569.0000 

912.0000 

0 

109.0000 

5632.0000 

8952.0000 

FACE 

A 

D C 

D 

ROT 

FR 

1234 

1.00000 

0.00000  0.00000  99.00000 

0.00 

0.09 

5673 

1.00000 

0.00000  0.00000  109.00000 

0.00 

0.00 

THICKNESS- 

10 

PT 

X 

Y 

7 

i. 

1 

23.0000 

64.0000 

100.0000 

O 

#• 

154.0000 

237.0000 

100.0000 

3 

564.0000 

278.0000 

ion. oooo 

4 

1254.0000 

6587.0000 

ion. noon 

1> 

23.0000 

64.0000 

110.0000 

0 

154.0000 

237.000') 

1 10.0000 

7 

564.0000 

270.0000 

lio.nooo 

0 
* 1 

1254.0000 

6587.0000 

110.0000 

FACE 

A 

R C 

D 

ROT 

FR 

1234 

0.00000 

0.00000  -1.00000  -100.^1000 

******* 

.nr  nr) 

6673 

0.00000 

0.00000  -1.00000  -110.00900 

******* 

-90. *00 

***  !!0Tr.  Slice  FC  IS  -90  ROT  IS  HOT  UNIOIIE*** 
THICKNESS-  -10 
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12.  Program  Name:  NORMVEC 

Description:  This  program  finds  (1)  a vector  of  any  desired 

length  in  the  same  direction  as  a given  vector, 
and  (2)  a vector  of  any  desired  length 
perpendicular  to  the  given  vector  (2  dimensions). 
For  a 3 dimensional  case,  use  program  13,  PERPENV. 


Memory:  802  bytes 

Restrictions:  User  must  beware  the  direction  of  the  perpendi- 

cular vector.  If  the  one  printed  out  is  not 
the  one  desired,  use  its  "negative"  or  run  the 
problem  again  with  negative  length  (see  sample 
outputs  3 and  4). 

Instructions:  After  the  program  is  loaded,  enter  the  following: 

(enter  the  data  separated  by  commas) 

(a)  X and  Y components  of  the  given  vector, 
length  of  vector  in  direction  of  given 
vector,  length  of  perpendicular  vector  CR/LF 

(b)  Results  will  be  printed  out 

(c)  To  run  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 

Sample  Outputs: 


INPUT  VECTOR: 
i rriftTi!«  inn 


LENGTH- 

100.000 

DEL  X- 

0.000 

DEL 

Y« 

100.000 

VECTOR 

III  DIRECTION  OF  INPUT  VECTOR: 

LENGTH* 

10.000 

DEL  X- 

o.ono 

DEL 

Y« 

10.000 

VECTOR  PERPENDICULAR  TO  INPUT  VECTOR: 

LENGTH- 

25.000 

DEL  X- 

-25.000 

DEL 

Y» 

0.000 

INPUT  VECTOR: 

LENGTH- 

141.  <121 

DEL  X- 

100.000 

DEL 

Y* 

100.000 

VECTOR 

IN  DIRECTION  OF  INPUT  VECTOR: 

LENGTH- 

10.000 

DEL  X- 

7.071 

DEL 

Y- 

7.071 

VECTOR  PERPENDICULAR  TO  INPUT  VECTOR: 

LENGTH- 

25.000 

DEL  X- 

-17.077 

DEL 

Y- 

17.677 

30 


INPUT  VECTOR: 

LENGTH-  261.178  DEL  X-  123.456  DEL  Y-  -230.157 

VECTOR  III  DIRECTION  OF  INPUT  VECTOR: 

LENGTH-  12.350  DEL  X-  5.037  DEL  Y-  -10.833 

VECTOR  PERPENDICULAR  TO  INPUT  VECTOR: 

LENGTH-  36.200  DEL  X-  31.D00  DEL  Y«  17.111 


INPUT  VECTOR: 

LENGTH-  261.178  DEL  X-  123.456  DF.L  Y-  -230.157 

VECTOR  IN  DIRECTION  OF  INPUT  VECTOR: 

LENGTH-  12.350  DEL  X-  5.037  DEL  Y-  -10.003 

VECTOR  PERPENDICULAR  TO  INPUT  VECTOR: 

LENGTH-  -36.200  DEL  X-  -31.000  DEL  Y-  -17.111 


31 


— ' -A &ml  iL  2U. 


13.  Program  Name:  PERPENV 

Description:  This  program  finds  (1)  a vector  of  any  desired 

length  perpendicular  to  two  given  vectors  and 
(2)  the  angle  between  the  two  given  vectors. 


Memory:  1373  bytes 

Restrictions:  User  must  beware  of  the  direction  of  the 

perpendicular  vector.  The  program  prints  out 
both  of  the  perpendicular  vectors  and  the  user 
must  choose  the  one  desired. 

Instructions:  After  the  program  is  loaded,  enter  the  following: 

(enter  the  data  separated  by  conus) 

(a)  The  X,Y  and  Z components  of  vector  1 
CR/LF 

(b)  The  X,Y  and  Z components  of  vector  2 
CR/LF 

(c)  The  desired  length  of  the  perpendicular 
vector  CR/LF 

(d)  Results  will  be  printed  out 

(e)  Program  returns  to  step  (a) 

Sample  Outputs : 


iff  *-vtcT?!fttooo  DEL  y-  o.ooon  no.  z»  o.onoo 
nn  v-  o.ooon  del  y»  in. noon  pel  7-  o.nnno 
LENGTHS  VECTOR  !■  15.0000  VECTOR  > in. 0000 

AIIGLE  HETWEEIi  INPUT  VECTORS-  90.000  DEGREES 


VECTOR  PERPENDICULAR  TO  INPUT  VECTORS: 

ox-  it.  oooo  dy-  o.onno  nz»  zs.onon 

OR  DX-  0.0000  DY-  0.0000  D7-  -PC. oooo 


INPUT  VECTORS: 

DEL  X-  155.2340  DEL  Y-  -52.3147  DEL  7-  -91.2400 

DEL  X-  -01.. 0357  DEL  Y-  -35.32K  DEL  7-  .092.35*0 

LEI.OTIIS  VECTOR  1-  191.3554  VECTOR  2-  0°5.1497 
ANGLE  r.ET'.'EEN  INPUT  VECTORS-  121.219  DEGREES 

VECTOR  PERPENDICULAR  TO  INPI*T  VECTORS: 

DX-  -50.2405  DY-  -137.0390  DZ- 

OR  DX-  GO. 2405  DY-  137.0390  D7- 
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-9.5501 

9.5501 
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14.  Program  Name: 
Description: 


AMT RACK 

This  program  adds  tracks  to  domestic  vehicles, 
or  any  vehicle  with  "live"  tracks.  The  two 
figures  below  are  drawings  of  the  track  the 
program  produces.  Included  in  the  output  are 
the  solid  table  and  the  region  combination 
table. 


Memory:  7497  bytes 

Restrictions:  The  following  assumptions  are  made: 

(1)  The  Z-coordinates  of  all  the  road  wheels 
are  equal 

(2)  The  "idler"  wheel  is  to  the  front  of  the 
first  road  wheel 

(3)  The  "drive"  wheel  is  to  the  rear  of  the 
last  road  wheel 

(4)  The  radii  of  the  road  wheels  are  equal 

Instructions:  After  the  program  is  loaded,  enter  the 

following: 

(enter  the  data  separated  by  comm) 

(a)  The  X,Z  coordinates  of  the  first  road 
wheel  CR/LF 

(b)  The  X,Z  coordinates  of  the  last  road 
wheel  CR/LF 

(c)  The  radius  of  the  road  wheels  CR/LF 

(d)  The  X,Z  coordinates  of  the  "idler" 
wheel  CR/LF 
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(e)  The  radius  of  the  "idler"  wheel  CR/LF 

(f)  The  X,Z  coordinates  of  the  "drive"  wheel 
CR/LF 

(g)  The  radius  of  the  "drive"  wheel  CR/LF 

(h)  The  Y min  and  Y max  of  the  tracks  CR/LF 

(i)  The  thickness  wanted  for  the  tracks  CR/LF 

(j)  Results  will  be  printed  out 

(k)  Program  returns  to  step  (a). 

Sample  Output:  The  following  input  data  were  used  to  produce 

this  sample  output  (according  to  the  order 
of  input  required  in  the  instructions  above:) 


(1) 

1500,100 

(2) 

-1600,100 

(3) 

80 

(4) 

2000,500 

(5) 

60 

(6) 

-2200,600 

(7) 

75 

(8) 

0,100 

(9) 

20 

1 BOX 

1539.3963 

0.0000 

0.0000 

512.6554 

0.0000 

439.2405 

FRHT  SLOPE 

-13.0129 

0.0000 

15.1076 

0.0000 

100.0000 

0.0000 

2 BOX 

•1639.3383 

0.0000 

0.0000 

-622.1040 

0.0000 

527.54*3 

REAR  SLOPE 

12.9352 

0.0000 

15.2530 

0.0000 

100.0000 

0.0000 

3 RP*» 

-1639.3383  1539.3963 

0.0000 

100.0000 

0.0000 

29.0000 

TRACK  BOT 

4 ARBR 

2247.2457 

110.0000 

630.0000 

2247.2457 

110.0000 

211. *353 

IDLER  Dir 

2000.0000 

110.0000 

500.0000 

2000.0000 

110.0000 

630.0900 

# 

2247.2457 

-10.0000 

630.0000 

2247.2457 

-10.0000 

211.4353 

2000.0000 

-10.0000 

500.0000 

2000.0000 

-10.0000 

639.0000 

5 ARBS 

-2455.4711 

110.0000 

745.0000 

-2455.4711 

110.0000 

290.7362 

DRIVE  DUI 

-2200.0000 

110.0000 

600.0000 

-2200.0000 

110.0000 

745.0000 

-2455.4711 

-10.0000 

745.0000 

-2455.4711 

-10.0000 

208.7362 

•2200.0000 

-10.0000 

600.0000 

-2200.0000 

-10.0000 

745.0000 

6 RCC 

2000.0000 

0.0000 

500.0000 

0.0000 

100.0000 

0.0900 

TUCK  IDLER 

80.0000 

7 RCC 

•2200.0000 

0.0000 

600.0000 

0.0000 

100.0000 

0.0000 

TRCK  DRIVE 

95.0000 

8 COX 

2009.9962 

100.0000 

559.7262 

-*219.9925 

o.oooo 

115.5*74 

TRACK  TOP 

0.5474 

0.0000 

19.9925 

0.0000 

-100.0000 

0.0000 

9 RCC 

2000.0000 

0.0000 

500.0000 

0.0000 

100.0000 

O.OOOO 

IDLER  BUI 

60.0000 

10  RCC 

-2200.0000 

0.0000 

600.0000 

0.0000 

100.0000 

0.0000 

DRIVE  mil 

75.0000 

REGION  TABLE 

1 

1 -4  0 

Fwrr  slope 

2 

2-5  0 

REAR  SLOPE 

3 

3 -1  -2 

TRACK  BOT 

4 

6 4-9 

TRCK  IDLER 

5 

7 5 -10 

TRCK  naiVE  . 

C 

8 -4  -5 

TRACK  TOP 
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15.  Program  Name:  SOLIDROT 

Description:  This  program  rotates  the  COM-GEOM  solids  (except 

ARS)  about  any  point  in  the  XY,  XZ,  or  YZ  planes. 
Positive  rotation  is  from  positive  axis  to 
positive  axis  - that  is,  positive  rotation  in 
the  XZ  plane  is  from  the  positive  X axis  towards 
the  positive  Z axis.  User  must  beware  of  the 
direction  of  the  rotation. 


Memory: 

Restrictions: 


Instructions : 


7539  bytes 

The  COM-GEOM  solid  ARS  (triangular  surfaced 
polyhedron)  cannot  be  rotated  by  this  program. 
The  large  amount  of  input  required  was  the 
reason  the  ARS  was  not  included.  If  needed 
however,  it  would  not  be  that  difficult  to 
add. 

All  rectangular  parallelpipeds  (RPP)  are 
changed  to  boxes  (BOX)  whether  or  not  the 
rotation  is  such  that  they  remain  RPPs. 

After  the  program  is  loaded,  enter  the 
following: 

(a)  Plane  of  rotation  where 


CR/LF 


(b) 

(c) 

(d) 


(e) 

(f) 


1 ■ XY  plane 

2 ■ XZ  plane 

3 ■ YZ  plane 
Angle  of  rotation  (in  degrees)  CR/LF 
X,Y,Z  coordinates  of  point  about  which 
you  wish  to  rotate  the  solid  CR/LF 
Next  input  the  solid  type  you  wish  to 
rotate,  following  the  instructions  on 
the  scope. 

Results  will  then  be  printed  out 
To  nut  another  case,  key  CONTINUE  CR/LF 
and  program  will  return  to  step  (a). 
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Sample  Outputs:  The  first  four  sample  runs  indicate  the  direction 
of  positive  rotation  for  the  three  planes  of 
rotation  available. 


ANGLE  OF  ROTATION-  00.0000  DEG  IN  THE  XY  PLANE 
PT  AROUND  WHICH  SOLID  HAS  ROTATED  X-  0.0000  Y« 
INPUT  SOLID 

SPH  lO.O'tOO  0.0000  0.0000  1.0000 

ROTATED  SOLID 

SPH  0.0000  10.0000  0.0000  1.0000 


0.0000 


0.0*00 


ANGLE  OF  ROTATION-  00.0000  DEO  IN  THE 
PT  AROUND  WHICH  SOLID  WAS  ROTATED  X- 
INPUT  SOLID 

SPH  10.0000  0.0000  0.0000 

ROTATED  SOLID 

SPH  0.0000  0.0000  10.0000 


XZ  PLANE 

0.0000  Y»  0.0000  7-  0.0000 

1.0000 

1.0000 


ANGLE  OF  ROTATION-  90.0000  DEG  IN  THE  YZ  PLANE 
PT  AROUND  WHICH  SOLID  WAS  ROTATED  X-  0.0000  Y« 
INPUT  SOLID 

SPH  0.0000  10.0000  0.0000  1.0000 

ROTATED  SOLID 

SPli  0.0000  0.0000  10.0000  1.0000 


0.0000 


0.0000 


ANGLE  OF  ROTATION-  90.0000  DEG  IN  THE 
PT  AROUND  Illicit  SOLID  HAS  ROTATED  X- 
INPUT  SOLID 

SPH  0.0000  0.0000  10.0000 

ROTATED  SOLID 

SPli  1.0000  -10.0000  0.0000 


ANGLE  OF  ROTATION-  31.2900  DEG  IN  THE 
PT  AROUND  WHICH  SOLID  HAS  ROTATED  X- 
i;;pirr  SOLID 

RPP  -5G.235S  25.3247  -907.2354 

ROTATED  SOLID 

DOX  -5(1.2350  -1424.3*02  1670.2941 
0.0000  544.5971  330.4595 


YZ  PLANE 

0.0000  Y«  0.0000  Z-  0.0000 

1.0000 

l.onon 

YZ  PLANE 

25.0000  Y»  35.2140  7-  1235.254* 

-250.214*  2255.0000  3214.230* 

P.1.5505  0.0000  0.0000 

0.000*  -440,5570  725.91^ 
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Program  Name: 
Description: 


Memory: 
Restrictions : 


Instructions : 


PLOTSOL 

This  program  produces  scaled  plots  of  selected 
COM-GEQM  solids.  The  following  solids  can  be 
plotted  by  this  program:  RPP,  BOX,  RAN, 

ARBS,  ARB7,  ARB6,  ARBS,  ARB4,  RCC,  TRC  and 
SPH.  Features  of  this  program  are  as  follows: 

(1)  Ability  to  replot  the  same  solids  at  a 
different  aspect  (view). 

(2)  Ability  to  add  solids  as  long  as  the 
limits  on  the  maximum  storage  requirements 
are  not  violated  (see  restrictions  below). 

(3)  As  each  solid  is  entered,  the  user 
acknowledges  that  the  input  is  correct 
before  entering  the  next  solid.  This 
feature  helps  eliminate  input  errors. 

(4)  The  scale  and  azimuth  and  elevation 

angle  are  printed  out  on  each  view  plotted. 

(5)  The  program  gives  the  user  the  opportunity 
to  adjust  the  lengths  of  the  axes  so  they 
will  be  equal.  If  the  axes'  lengths 

are  equal  then  the  scale  printed  out  would 
be  correct  and  circles  will  look  like 
circles. 

32053  bytes 

The  following  are  the  limitations  and  restrictions 
of  this  program: 

(1)  The  sum  of  the  number  of  TRCs  and  RCCs 
to  be  plotted  must  not  be  greater  than 
5. 

(2)  The  maximum  number  of  SPHs  that  can  be 
plotted  is  5. 

(3)  The  storage  requirements  for  the  RPP,  BOX 
ARB8  U 8;  for  the  RAN  and  ARB6  is  6; 
for  the  ARB5  is  5;  and  for  the  ARB4  is  4. 
Using  these  storage  numbers,  the  sum  of 
all  storage  requirements  for  RPPs,  BOXs, 

ARBs  and  RA Ns  must  not  exceed  120.  For 
example,  the  maximum  number  of  ARB4s  that 
can  be  plotted  is  30. 

(4)  This  program  contains  no  hidden  line 
routine. 

(5)  This  program  takes  about  55  seconds  to 
process  a TRC  and  about  35  seconds  for 
a RCC. 

After  the  program  is  loaded,  perform  the 
following  steps,  following  the  instructions 
on  the  screen: 

(a)  Ready  the  plotter 

(b)  Enter  the  number  of  solids  to  plot  CR/LF 

(c)  Enter  the  azimuth  and  elevation  angle  for 
this  view  CR/LF 


Sample  Outputs: 


(d)  Enter  the  solid  type  using  the  code 
displayed  on  the  screen  (1  - RPP,  2 = BOX, 
etc.)  CR/LF 

(e)  Enter  the  parameters  for  the  solid  as  they 
are  asked  for 

(f)  When  finished  entering  the  solid  parameters, 
enter  1 if  data  is  correct,  0 if  not. 

If  data  is  correct  program  processes  the 
data.  If  data  is  incorrect  program  returns 
to  step  (e) . 

(g)  Enter  the  remaining  solids  using  the 
scheme  above  (steps  (d),  («),  and  (f)). 

(h)  When  solids  are  all  entered,  the  program 
asks  if  the  axes  lengths  are  adjusted. 

The  first  time  through  the  program,  the 
user  should  indicate  NO  (by  entering 
zero)  and  two  perpendicular  lines  will 
be  drawn  by  the  plotter.  If  these  lines 
are  of  equal  length  then  the  axes  are 

of  equal  length.  If  these  lines  are  not 
equal,  then  adjust  them  using  the  scale 
adjust  controls  on  the  plotter.  The 
program  will  again  ask  if  axes  lengths 
are  adjusted  and  the  user  continues  to 
say  NO  and  adjust  until  they  are  equal. 

(i)  When  axes  are  adjusted  the  user  should 
change  paper  on  the  plotter  and  key 
CONTINUE  CR/LF  to  plot. 

(j)  Once  the  solids  are  plotted,  indicate 
if  another  view  of  the  same  solids  is 
desired,  if  yes  program  returns  to  step 

(i). 

(k)  If  another  view  is  not  desired,  the  user 
must  next  indicate  if  he  wishes  to  add 
more  solids'.  If  yes,  indicate  the  number 
to  add  and  the  program  returns  to  step 
(d).  If  not,  the  program  returns  to 
step  (a). 

Figure  1 shows  individual  plots  of  the  11 
COM-GEOM  solids  which  this  program  plots. 

These  plots  are  intended  to  familiarize  the 
user  with  the  way  these  solids  look  when  plotted. 
Figure  2 is  a sample  plot  of  S solids  at  90* 
azimuth  and  0*  elevation.  Figure  3 is  a plot 
of  the  same  solids  at  a different  aspect  (125* 
azimuth  and  18*  elevation). 


Cbpf  avcilablo  to  DDC  <••'■''  r >J 

jMnait  fully  legible  rcpro.'x  on 


TRC 


\ 

\ 


SPH 


Figure  1.  Staple  Plots  of  the  Available  COH-GEOM  Solids  (Continued) 
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Figure  3.  Sample  Output  of  the  PLOTSOL  Program 


17.  Program  Name: 
Description: 


Memory: 

Restrictions: 


Instructions: 


DARBIN 

This  program  corrects  an  Arbitrary  Convex 
Polyhedron  (ARB)  which  has  a fourth  vertice 
not  in  a plane  defined  by  the  other  three 
vertices  and  computes  an  inside  ARB  given 
a thickness  for  each  face. 

3432  bytes 

When  solving  for  the  error  denoted  by  the 
message  "POUR  POINTS  NOT  CO-PLANAR  IN  FACE 
1234",  the  program  will  change  the  X,Y  and 
Z of  point  4 to  correct  the  ARB.  To  change 
the  X,Y  and  Z of  point  1,  enter  the  face 
ordinal  numbers  2341  instead  of  1234. 

After  the  program  is  loaded,  enter  the 
following: 

(a)  The  number  of  vertices  CR/LF 

(b)  The  number  of  faces  CR/LF 

(c)  The  X,Y,  and  Z coordinates  of  each 
vertice  CR/LF 

(d)  A four  digit  ordinal  number  for  each 
face  CR/LF.  (Foce  512  would  be 
entered  as  5120)  Figure  4 contains 
face  ordinal  numbers  which  may  be  used. 

(e)  The  option  indicator 

Indicator  Option 


1 compute  inside  ARB 

0 return  to  step  (a) 

CR/LF 

(f)  Thickness  for  each  face  CR/LF 

(g)  The  program  returns  to  step  (a). 

Sample  Outputs: 


l»  jp*:t 


0.0000 

fl.oono 

n.ooon 

o.nooo  • 

5.0000 

0.0000 

Viono  • 

5.  moo 

o.onnn 

5.0000 

0.0000 

O.'V'OO 

o.noon 

o.onnn 

5.0000 

o.onnn 

5.0000 

r».onoft 

•i.onm 

i?34  r»r,7u  r*7u 

s.nono 
I4p,r»  i,?cr» 

a. noon 
4070 

5.0000 

0.0000 

5. 0000 

7i:ir.i;:;R5S  ARE 

;..loor, 

O.ROOO 

0.3000 

-n.4nnn 

0.5000 

u.r.^oo 

icsnc  am 

-0.4000 

0.1 000 

0.5003 

4.7000 

3.1000 

. •(IIih; 

4,7000 

o.i non 

4.4000 

-0.4000 

0.1000 

— 0.4*lij0 

4.0030 

0.5000 

4.7000 

4.0000 

4.7000 

/i.oonn 

4.4000 

-'VO'V) 

4.0^00 
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l . >a,t  I W ■ 


g- 


1.0000 

0.0003 

j.0000 

o.nooo 

5.  oo 00 

0.0000 

5.0000 

5.0000 

0.0000 

5.0000 

■•.00,00 

f)  rr^ri 

0.0000 

0.0000 

5.0000 

0.0000 

5.0000 

5 ! 0000 

5.0000 

5.0000 

5.0000 

5.0000 

0.0000 

5.0000 

1234  507;;  r>7.''.  1/105  1200 
FOUR  points  i!0T  CO-PLANAR 

4570 

IN  FACF  1 

2 G 5 ni!= 

-0.0005 

SOLUTION 

1.0000 

I).  000,0 

0.0000 

o.nooo 

5.000' 

0.0000 

G.Onnfl 

5.0000 

0.0000 

5.0000 

0.0000  , 

O.oooo 

1.0000 

o.oono 

s.ooon 

'.0000 

5. 0000 

5.0000 

5.0000 

5.0000 

5.0000 

5.0000 

O.OONO 

5.0000 

INPUT 

1.0000 

•-.0000 

o.oono 

o.oono 

5.0000 

0.0000 

O.OOOO 

5.0000 

0.0000 

5.0000 

'.'I.OO^o 

j.OOOO 

;.‘T)oo 

0.0003 

;;.onoo 

o.oono 

5,0000 

r».r»oo 

:;.timo 

5. 0000 

5,0000 

5,0000 

o.oono 

5 . 0'OOO 

r.  V-1  5070  2.070 
FOUR  POINTS  NOT 

4501  2051 
CO-PLANAR 

4570 

IN  FACT.  2 

G 5 1 Dl!= 

n.^pflp 

SOLUTION 

0.0000 

0.0000 

0.0000 

o.nooo 

5.0000 

n.oooo 

5.0,000 

5.0000 

0.0000 

5.0000 

0.0000 

0.0000 

0,0000 

o.nooo 

5,0000 

o.nooo 

5.0000 

5.0000 

5.0000 

5.0000 

5.0000 

5.0000 

o.oono 

G.oono 
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FACE  ORDINAL  NUMBERS  WHICH  CAN  BE  ENTERED 


NO.  OF 

NO.  OF 

VERTICES 

FACES 

1234 

3678 

1384 

2376 

1234 

S670 

1430 

2376 

1234 

2363 

1364 

3120 

1234 

3120 

3230 

3340 

1230 

4120 

4230 

4310 

Figure  4.  An  Arbitrary  Convex  Polyhedron  (AR3) 


Program  Name : BOX IN 

Description:  This  program  corrects  a Rectangular  Parallelpiped 

in  any  orientation  (BOX)  whose  vectors  are  not 
normal  and  computes  an  inside  BOX  given  a 
thickness  for  each  face. 


Memory: 
Restrictions : 

Instructions : 


Sample  Outputs: 


2717  bytes 

A corrected  non-normal  vector  maintains  its 

original  length. 

After  the  program  is  loaded,  enter  the  following: 

(a)  The  X,Y  and  Z coordinates  of  the  vertice 
CR/LF 

(b)  The  X,Y  and  Z components  of  the  length  (I) 
vector  CR/LF 

(c)  The  X,Y  and  Z components  of  the  width  (W) 

vector  CR/LF  _ 

(d)  The  X,Y  and  Z components  of  the  height  (H) 
vector  CR/LF 

(e)  The  option  indicator 

Indicator  Option 

1 compute  inside  BOX 

0 return  to  step  (a)  CR/LF 

(f)  The  thicknesses  for  each  face.  CR/LF 

(g)  The  program  returns  to  step  (a). 


ICP'T 


vrvyi 

u.otyi 

0.0000 

10.0000 

0.0000 

o,*oo* 

u.i’i'jOO 

10.0000 

i.nino 

n.oofy’ 

(.''no.) 

io!*roo 

TliICK'KSS  AIT 

;.10:/; 

-o.poon 

0.0000 

-o.4r/  ; 

o.r>ooo 

-o.r^ 

ivtax  s ox 

1.1  Oil . 

-J.."000 

0.1000 

1..1000 

0.0000 

li.  '“•000 

P/.OO-Kl 

0. 7000 

0.0000 

y.onoo 

' 1.000* 

1 0.1*00 
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INPUT 


0.0000 

0.0000 

0.0000 

10.0000 

0.0000 

O.nooo 

0.0000 

10.0000 

0.0000 

1.0000 

0.0000 

10.0000 

HEIGHT  /.HP  LENGTI!  VECTORS  ARE  ! 

IOT  NORMAL, 

ANGLE  = 

34.20 

ARC  SOLUTION 

0.0000 

0.0000 

0.0000 

10.0000 

o.nooo 

o.noon 

10.0000 

10.0000 

0.0000 

o.oono 

10.0000 

0.0000 

1.0000 

0.0000 

10.0000 

11.0000 

0.0000 

1 o.onoo 

11.0000 

1 11.0000 

10.0000 

1.0000 

lo.onno 

l o.onoo 

BOX  SOLUTION 
o.oono 

0.0000 

0.0000 

10.0000 

o.onoo 

0.0000 

0.0000 

10.0000 

0.0000 

0.0000 

0.0000 

in.o/'o" 

THICKNESS  ARE 
0.1000 

' 0.1000 

o.?.ono 

0.2000 

0.0000 

o.oono 

INSIDE  BOX 

0.1000 

0.1000 

0.2000 

0.7000 

o.onoo 

0.0000 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

9.5403 
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19.  Program  Name:  RAWIN 

Description:  This  program  corrects  a Right  Angle  Wedge  (RAW) 

whose  vectors  are  not  normal  and  computes  an 
inside  RAW  given  a thickness  for  each  face. 


Memory:  2872  bytes 

Restrictions:  A corrected  non-normal  vector  maintains  its 

original  length. 

Instructions:  After  the  program  is  loaded,  enter  the  following: 

(a)  The  X,Y  and  Z coordinates  of  the  vertice 
CR/LF 

(b)  The  X,Y  and  Z components  of  the  height  (H) 
vector  CR/LF 

(c)  The  X,Y  and  Z components  of  the  width  (W) 
vector  CR/LF 

(d)  The  X,Y  and  Z components  of  the  depth  (D) 
vector  CR/LF 

(e)  The  option  indicator 

Indicator  Option 

1 compute  inside  RAW 

0 return  to  step  (a)  CR/LF 

(f)  The  thickness  for  each  face  CR/LF 

(g)  The  program  returns  to  step  (a). 

Sample  Outputs: 

INPUT 


a. noon 

u.nono 

o.noon 

o.oonn 

o.nnon 

lo.onno 

j.ntyvi 

lo.nooo 

o.nnon 

lo. noon 

0.0000 

n.ooon 

THICKNESS  m 

0.2000 

0.4000 

0.0000 

0.0000 

-1.0000 

I:. ’SIDE  RAlf 

o.f.oon 

0.4000 

0.2000 

0.0000 

0.0000 

11.2142 

a. noon 

11.0142 

1 1.0000 

o.rooo 

o.nnon 

o.nnon 

r 


f 

f- 


J ■ 
| 


INPUT 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

10.0000 

0.0000 

10.0000 

0.0000 

10.0000 

0.0000 

0.0000 

DEPTH  AND  HEIGHT  VECTORS  ARE  NOT  NORMAL. 

ANGLE  * 

04.2  P. 

ARB  SOLUTION 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

10. 0000 

11.0000 

0.0000 

10.0000 

10.0000 

0.0000 

0.0000 

0.0000 

10.0000 

0.0000 

10.0000 

10.0000 

0.0000 

RAW  SOLUTION 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

10,0000 

0.0000 

10.0000 

0.0000 

9.950’ 

o.oonn 

-o.oor.n 

THICKNESS  ARE 
0.5000 

O.AOOO 

0.3000 

0.2000 

o.ioon 

INSIDE  RAW 

0.3482 

0.4000 

0.4576 

0.0361 

0.0000 

9.361A 

0.0000 

9.4589 

0.0000 

9.4520 

0.0000 

-0.OA52 

20.  Program  Name: 
Description: 


Memory: 

Restrictions: 

Instructions: 


Sample  Outputs: 


TRCIN 

This  program  computes  an  inside  Truncated  Right 
Angled  Cone  (TRC)  given  a thickness  for  each 
surface. 


After  the  program  is  loaded,  enter  the  following: 

(a)  The  X,Y  and  Z coordinates  of  the  vertice 
CR/LF 

(b)  The  X,Y  and  Z components  of  the  height 
vector  CR/LF 

(c)  The  radius  of  the  base  CR/LF 

(d)  The  radius  of  the  top  CR/LF 

(e)  The  thickness  for  each  surface  CR/LF 

(f)  The  program  returns  to  step  (a). 


IMPIJT 


o.onnn 

o.nnnn 

n.onnn 

o.onnn 

n.onnn 

1 o.nnnn 

10.0000 

5.0000 

THICKNESSES  ARF 

-l.nono 

2.0000 

o.ooon 

n:s  i nr  trc 

• ».  0000 

0.0000 

-l .nnno 

o.onnn 

o.nnnn 

o.nnnn 

7.14r»;> 

2.f.«3 

INPUT 

n.onoo 

o.onnn 

n.onnn 

0.0000 

o.onnn 

1 n.onnn 

10.0:. 

1.0000 

Tl  ICK.’OSES  ARE 

l.oono 

2.0000 

i.nonn 

!;AF;;;if:n,  top  radius  of 

ins I OF  TRC  « 

-1.23 

RAM  MS  OF SET  TO 

.unnl  AflO  THICKf.TSS(P) 

« 0 

ifisinr  trc 

0.0000 

o.onnn 

l.nono 

o.onnn 

o.nnnn 

G.r.?rn 

S.Of.TJ 

n.oooi 
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21.  Program  Name:  RECIN 

Description:  This  program  corrects  a Right  Angle  Elliptical 

Cylinder  (REC)  whose  vectors  are  not  normal  and 
computes  an  inside  REC  given  a thickness  for 
each  surface. 


Memory: 
Restrictions : 


Instructions : 


Sample  Outputs: 


B 

2337  bytes 

The  corrected  non-normal  vector  maintains  its 
original  length.  The  inside  elliptical  surface 
does  not  yield  a constant  thickness.  This  is 
especially  true  when  the  ratio  of  the  lengths 
of  the  minor  axis  to  the  major  axis  is  less  than 
0.8.  Therefore,  this  program  computes  an  inside 
REC  with  an  elliptical  surface  which  yields  an 
average  thickness  equal  to  the  thickness  entered. 
After  the  program  is  loaded,  enter  the  following: 

(a)  The  X,Y  and  Z coordinates  of  the  vertice 
CR/LF 

(b)  The  X,Y  and  Z components  of  the  height 
vector  CR/LF 

(c)  The  X,Y  and  Z components  of  the  defining 
the  semi-major  axis  CR/LF 

(d)  The  X,Y  and  Z components  of  the  vector 
defining  the  semi-minor  axis  CR/LF 

(e)  The  option  indicator 

Indicator  Option 

1 compute  inside  REC 

0 return  to  step  (a)  CR/LF 

(f)  The  thickness  for  each  surface  CR/LF 

(g)  The  program  returns  to  step  (a). 


if:rn 

11.0000  0.0000  0.0000 

10.0000  0.0000  0.0000 

THICKfJFSS  ARK 


0.0000  0.0000  10.0000 

0.0000  5. 0000  0.0000 


SI 


1.0000 

2.0000 

-3.0000 

INSIDE  REC 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

7.0000 

13.0000 

0.0000 

0.0000 

0.0000 

0.1021 

0.0000 

INPUT 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

10.0000 

10.0000 

1.0000 

0.0000 

0.0000 

5.0000 

0.0000 

SEMI-MAJOR  AMO  SEMI-MINOR  VECTORS 

NOT  NORMAL.  ANGLE- 

04.23 

SOLUTION 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

10. 0000 

10,0400 

0.0000 

0.0000 

0.0000 

5.0000 

0.0000 

THICKNESS  ARE 

1.0000 

2.0000 

-3.0000 

INSIDE  REC 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

7.0000 

13.0490 

0.0000 

0.0000 

0.0000 

0.1032 

0.1000 
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22.  Program  Name:  TECIN 

Description:  This  program  corrects  a Truncated  Elliptical 

Cone  (TEC)  whose  base  vectors  are  not  normal  and 
computes  an  inside  TEC  given  a thickness  for 
each  surface. 


Memory:  3432  bytes 

Restrictions:  The  corrected  non-normal  vector  maintains  its 

original  length.  The  elliptical  surface  does 
not  yield  a constant  thickness.  This  is 
especially  true  when  the  ratio  of  the  lengths 
of  the  minor  axis  to  major  axis  is  less  than 
0.8.  Therefore,  this  program  computes  an  inside 
TEC  with  an  elliptical  surface  which  yields  an 
average  thickness  equal  to  the  thickness  entered. 

As  a result  of  the  restriction  of  defining  the 
top  of  a TEC  as  a ratio  of  the  base,  the  inside 
TEC  is  a best  fit.  Therefore,  a TEC  with  a 
large  ratio  or  a relatively  large  thickness 
relative  to  the  length  of  the  minor  axis  or  a 
small  minor  axis  relative  to  the  major  axis  may 
have  a large  variation  of  thicknesses  on  its 
elliptical  surface. 

Instructions:  After  the  program  is  loaded,  enter  the  following: 

(a)  The  X,Y  and  Z coordinates  of  the  vertice 
CR/LF 

(b)  The  X,Y  and  Z components  of  the  height 
vector  CR/LF 

(c)  The  X,Y  and  Z components  of  the  vector 
defining  the  semi-major  axis  of  the  base  CR/LF 


Sample  Outputs: 


(d)  The  X,Y  and  Z components  of  the  vector 
defining  the  semi -minor  axis  of  the  base  CR/LF 

(e)  The  ratio  of  base  to  the  top  CR/LF 

(f)  The  option  indicator 

Indicator  Option 

1 compute  inside  TEC 

0 return  to  step  (a)  CR/LF 

(g)  The  thicknesses  for  base,  top  and  elliptical 
side  in  exact  order  CR/LF 

(h)  The  program  returns  to  step  (a). 


INPUT 


1.0000 

0.0000  0.0000 

o.oono 

5.0000 

10. noon 

l!.*.rt0O0 

o.oono  o.oooo 

0.0000 

f>.oonn 

O.oooo 

,3.0000 

ti!Ick::fss  ai:f 

-0.5000 

i.nooo  2.oono 

rsiDF.  ire 

0.0000 

0.0004  -0.5000 

0.0000 

4.75°o 

9.5009 

0.01  O'.) 

0.0000  0.0000 

0.0000 

.3.652" 

o.onoo 

4.0607 

RATIO  OF  MAJOR  AXIS  ,3.4552  RATIO  OF  MINOR  , 

AXIS  5.0101 

TCP  MAJOR  AXIS 

3.2639  TOP  MINOR  AXIS 

0.52"" 

I'.Pi'T 

o.nonn 

0.O000  0.0000 

0.0000 

5.0000 

in.ooon 

lu.0000 

5.0000  (1.0000 

0.0000 

5.0000 

o.onoo 

1.1260 

SEMI -MAJOR  AND 

SOU-MINOR  VECTORS  NOT 

NORMAL, 

AN0I.F  » 63.43 

SOLUTION 

0.0000 

0.0000  0.0000 

0.0000 

5.0000 

10.0009 

11.1003 

0.0000  0.0000 

0.0000 

G.nono 

o.oono 

1.1,350 

TKICI'NFSS  ARC. 

0.1 000 

3.1000  0.6000 

I NS I HE  TFC 

0.01)00 

0.2647  0.1000 

0.0000 

4.0900 

o.onoo 

10.0600 

1.0000  0.0000 

0,0000 

4.1331 

0.0000 

l.nai 

RATIO  OF  riAJtm  AXIS  1.130,?  RATIO  OF  MINOR  AMS  1.1133 
TOP  lA.10.rl  AX  I!?  0.345.".  TOP  MINOR  AX  If*  3.7120 
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23.  Program  Name: 
Description: 


Memory: 

Restrictions: 


Instructions : 


t 


r 

i 


PARB 


This  program  computes  the  points  of  intersection 
of  a set  of  planes  taken  three  at  a time.  Each 
plane  can  be  defined  by  either  three  points; 
azimuth  angle,  elevation  angle  and  point;  or  an 
equation. 

4026  bytes 

This  program  is  limited  to  a maximum  of  25  planes. 
Combinations  of  three  planes  containing  parallel 
planes  or  whose  point  of  intersection  distance 
from  the  origin  is  greater  than  100,000  is  not 
computed  or  printed.  The  plane  equation 
(AX  + BY  ♦ CZ  ■ D)  is  normalized  so  that  a 
parallel  plane  can  be  obtained  by  only  changing 
the  value  of  D.  Furthermore,  the  difference 
between  the  value  of  D in  the  original  equation 
and  the  final  equation  is  the  distance  between 
the  planes. 

After  the  program  is  loaded,  enter  the  following: 

(a)  The  number  of  planes  CR/LF 

(b)  The  type  of  input  indicator  for  each  plane 

Indicator  Input 


0 azimuth,  elevation  and 

point 

1 three  points 

2 plane  equation  CR/LF 
(c)  If  the  indicator  keyed  in  step  (b)  is  0, 

enter  the  following: 


1.  The  azimuth  angle  (about  z axis)  CR/LF 

2.  The  elevation  angle  (about  y axis)  CR/LF 

(d)  If  the  indicator  keyed  in  step  (b)  is  1,  enter 
the  following: 

The  X,Y  and  Z coordinates  for  each  of  three 
points  in  the  plane  CR/LF 

(e)  If  the  indicator  keyed  in  step  (b)  is  2,  enter 
the  following: 

The  A,B,C,D  coefficients  of  the  plane 
equation  AX  ♦ BY  ♦ CY  ■ D CR/LF 

(f)  The  program  returns  to  step  (b)  until  all 
planes  are  entered. 

(g)  The  option  indicator 

Indicator  Option 

0 change  planes  previous 

entered 

1 return  to  step  (a)  CR/LF 

(h)  The  number  of  the  plane  to  be  corrected  CR/LF 
If  zero  is  keyed,  the  program  computes  the 
intersections  of  the  planes  and  returns  to 
step  (g).  If  a number  greater  than  zero  is 
keyed,  the  program  returns  to  step  (b). 
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Sample  Outputs: 


INPUT  FOR  FACE  1 

AZIMUTH  ■ 30.00  ELEVATION  ■ 26.00 

POINT  ■ 120.0000  245.0000  64.0000 

PLANE  EQ.  « 0.77837  X+  0.44939  Y+  0.43837  Z ■ 231.56 


INPUT  FOR  FACE  2 
3-POINTS 

-10.0000  0.0000  0.0000  -10.0000  2.0000  5.0000 

-10.0000  30.0000  58.0000 

PLANE  EQ.  *-1.00000  X+  0.00000  Y+  0.00000  7.  » 10.00 

INPUT  FOR  FACE  3 

PLANE  EO.  » 0.00000  X+  0.55470  Y+  0.83205  Z ■ 126.47 


INPUT  FOR  FACE  4 


AZIMUTH  * 
POINT  » 
PLANE  EO. 
PLANES 
1 2 3 

1 2 * 
1 3 4 

3 / 


O 

i. 


210.00  ELEVATION  « -09.00 
-10.0000  30.0000  58.0000 

* 0.01511  X+  0.00872  Y+  0.99984  Z » 
POSSIBLE  SOLUTION 
-10.0000  1099.0472  -580.6981 

-10.0000  479.8491  54.0739 

132. 2346  146.8891  54.0739 

-10.0000  142.472*  57.0184 


53.10 


INPUT  FOR  FACE  4 

PLANE  EO.  * 0.01511  X+  0.00872  Y+  0.99984  Z - 57.10 


PLANES 

POSSIBLE  SOLUTION 

1 

2 

3 

-10.0000 

1099.0472  - 

580.6981 

1 

2 

4 

-10.0000 

480.3316 

53.0667 

1 

3 

4 

181.9292 

148.4013 

53.0655 

O 

4. 

3 

4 

-10.0000 

143.9933 

56.004/. 
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APPENDIX  A 

Geometry  Aid  Program  Listings 


The  character  set  of  the  typewriter  used  to  produce  these 
program  listings  is  such  that  the  following  corrections  must  be 
made  by  the  user: 

1.  ! must  be  replaced  by  + 

2.  [ must  be  replaced  by  < 

3.  ] must  be  replaced  by  > 


* > 
f 


i 

i 

t. 


a 

r % 


Table  A-I.  SPHERE  Program  Listing 


10  DIM  P(4,3),C(2,3),!>(3) 

20  select  print  005(04) :print  "sphere  program" 

00  PRINT  "Finns  SPHERE  DEFINED  vy  4 UONPLANAR  POINTS":  PRII.7  "ON  SL'RFAC 
I.  OF  THE  SPHERE" 

4.1  FOR  I-1TO  4 

00  SELECT  PRINT  *>OD(u4) tPRIIlT  "POINT", I: INPUT  "X,Y,Z»",P(I,1),P(I,2),P 
(I,1):NEXT  I 

CO  SELECT  PRINT  215(30)  :PRII1T  "INPUT  POINTS" 

7i  print  "pt  x v 7* 

00  FOR  I-1T0  ^:Piii;iTHSIIIG  30,I,P{I,1),P(I,2),P(I,3):NEa7  I 

jtir;  -wm.tm 

no  for  i*i to  i:o(i)-o 
no  FUR  J-1TO  1:C(I,J)-2*(P(4,J)-P(I,J)) 

120  n(i)-P(I)+P{/,.n)!2-P(I,J)!2:tlEXT  JiHEXT  I 

111  FOR  K-1T0  2: IF  C(K,K)"OTNEil  ICO 

140  FOR  I-KTO  2:F»C(I+l,K)/C(K,U):n(I+1)»n(I+1)-F*n(I) 

ISO  FOR  ,1-lTO  3:C(I+1,J)-C(I+1,J)-F*C(KtJ):f!E.7T  ,1:?IEXT  I:OOTO  W 
KO  FOR  I-KTO  2: FOR  J-KTO  3:C(K,0K(K,J)+C(I+1  .jjillFRT  J 
17o  D(K)»D(K)+P(l+1):IF  C(K,K)[l0Tl:E!l  140:UEXT  1 
100  GOTO  .110 


ICO  NEXT  K 
201  IF  C(r»,l)-0TI 

210  Di2)-(n(3)-n 

22i  o(l)-(iHi)-n 
2.1  i n-OjFOR  J-1TI 
240  R-SQR(R) 

230  SELECT  PRINT  .?15(i»):PRIIITUSIi«  270,0(1), 0(2), 1(1) 

2*::  p:  iriTUSiun  ioo.r 

27r:sH:Eni.  ‘X*K****.M*i  %••?*&**.****  Km-w-if  ?,***? 

:is. 


tirn  310:D(3)»D(3)/C(3,1) 

|2i*c|i Izl-nf 3)*C  1,3))/C(1,1) 
0 3:R-R+(P(1,J).D(U))!2:HEX7  J 


.MINT  :PP.IHT  : SELECT  PRIilT  005(04) jPRIIJT  "TO  RUN  AGAIN,  KEY  CONTII! 


n i.xrcuTE" 


..v:  PRINT  "TO  STOP,  ICTY  CLEAR  EXECUTE": STOP  jOOTO  4f 
lie  PRINT  "POINTS  ARE  COPLANAR":GOTO  250 
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Table  A-II.  CIRCIR  Progra*  Listing 


10  SELECT  PRINT  005(64) :PRIMT  “CIRCIR  PROGRAM": PRINT  “FINDS  INTERSECTION 
OF  2 CIRCLES" 

20  INPUT  "X,Y  OF  CENTER  AND  RAO  OF  CIRCLE  1-“,A1,B1,R1 
30  IF  R1]0THEN  50:PRINT  "RADIUS[-0,  TRY  AGAIN" 

40  INPUT  "RADIUS-", R1:G0T0  30 

50  INPUT  "X,Y  OF  CENTER  AND  RAD  OF  CIRCLE  2-",A2,B2,R2 
60  IF  R230THEN  80:PRINT  "RADIUSC-O,  TRY  AGAIN" 

70  INPUT  "RADIUS-", R2:G0T0  60 

80  SELECT  PRINT  215(04):PRINTUSING  90,A1 ,B1 ,R1 : PRINTUSING  100,A2,B2,R2 
90%CIRCLE  1 X CENT-#####.####  Y CENT--#####.####  R -#####.#### 

100XCIRCLE  2 -#####.####  -#####.####  -#####.#### 

110  IF  B1-B2THEN  230 

120  M— (A2-A1)/(B2-B1):B3»(R1!2-R2!2-A1 12+A2I2-B1 12***B2 12)/(2*(B2-B1 ) ) 

130  A4-1+M!2:B4-2*M*(B3-B1 )-2*Al :C4»A1 !2+(B3-Bl) !2-Rl !2 
140  D-B4 ! 2- ( 4*A4*C4 ) : I F D[OTHEN  220 
150  X4r(-B4+SQR(Dj)/(2*A4) 

160  X5«(-B4-SQR{D))/(2*A4) 

170  Y4»M*X4+B3:Y5»M*X5+B3 

180  PRINTUSING  190,X4,Y4:PRIHTUSING  200.XG.Y5 

190%INTERSECTI0N  POINTS  X-#####.####  Y •-#####.#### 

200s  -#####.####  -#####.#### 

210  PRINT  .-PRINT  -.GOTO  20 

220  PRINT  "NO  INTERSECTION": PRINT  :PRINT  :GOTO  20 
230  IF  A1PA2TI1EN  240:  IF  R1-R2THEN  290 .-GOTO  220 
240  IF  ABS[A2-A1 )](R1+R2)THEH  220 
250  IF  ABS(A2-A1)-{R1+R2)THEN  300 

260  X4-(R1 12-R2I2-A1 !2+A2!2-Bl !2+B2I2)/(2*(A2-Al )):X5-X4 
270  Y4«B1+SQR(R1!2-  X4-A1  !2) 

280  Y5-B1 -SQR(R1 !2-(X4-Al j !2) :GOTO  180 
290  PRINT  "SAME  CIRCLE": PRINT  :PRINT  :GOTO  20 
300  X4-A1+R1 :Y4-B1 

310  PRINTUSING  320,X4,Y4:PRINT  :PRINT  :GOTO  20 
320‘,;CIRCLES  TANGENT  AT  X-#####.####  Y --#####.#### 

330  STOP  .-END 
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Table  A-III.  CIRCLE  Program  Listing 


40 

50 

60 

70 

80 


10  PRINT  "CIRCLE  PROGRAM": PRINT  "FINDS  CIRCLE  DEFINED  BY  3 NONCOLLINEA 
R POINTS" 

20  DIM  X(IO).YOO) 

30  FOR  J-1T0  3:PP!NT  "POINT" ,J: INPUT  "X-",X(J) 

INPUT  "Y»",Y(J):NEXT  J 
A«X(1)-X(2  :D»Y(1)-Y{2) 

C*X(2)-X(3):D*Y(2)-Y(3) 

E*(Y(1 ) J2+X(l ) !2-Y{2  !2-X(2) !2)/2 
F«(Y(2) !2+X(2) !2-Y(3) !2-X(3) !2)/2 
90  D1«A*D-C*B 

100  IF  D1  []0'THEN  110:PRINT  "^POINTS  COLLINEAR***":GOTO  30 
110  X1«(E*D-F*B)/D1 :Y1-(A*F-E*C)/D1 
120  R«SQR((X(1)-X1)!2+(Y(1)-Y1)!2) 

130  SELECT  PRINT  211(85) 

140  PRINT  " INPUT  POINTS" 

150  PRINT  "X  Y" 

160  FOR  1*1  TO  3: PRINTUSING  170,X(I),Y(I):NEXT  I 
170%-#####.####  -#####.#### 

130  PRINTUSING  190.X1.Y1.R 

1 902XCENT  *-##### . ####  YCENT-#####.####  RAD*-#####.#### 

200  PRINT  : PRINT 

210  SELECT  PRINT  005(64) 

220  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  EXECUTE" 

230  PRINT  "TO  STOP,  KEY  CLEAR  EXECUTE" 

240  STOP 
250  GOTO  30 
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Table  A-IV.  RCC  Program  Listing 


10  SELECT  PRINT  005(64) : PRINT  "RCC  PROGRAM" 

20  PRINT  "FINDS  RCC  WITH  BASE  DEFINED  BY  3 POINTS  ON  CIRCUMFERENCE" 

30  DIM  A(10.11  ).V(10),X(5),Y(5).7.(5),B(5) 

40  PRINT  "INPUT  THE  3 POINTS" 

50  FOR  1-1  TO  3: PRINT  "POINT". I 
60  INPUT  "X.Y.Z-",X(I).Y(I).Z(I):NEXT  I 
70  INPUT  "LENGTH  OF  HEIGHT  VECT0R\H:N«3 
80  IF  HjOTHEN  100 

90  PRINT  "*****  ERROR  .HEIGHT  [ ■ 0 *****":GOTO  70 
100  U1»X(2)-X(1):U2»Y(2)-Y(1):U3«Z(2)-Z(V 
110  Vl-X  3 -X(2  :V2-Y(3)-Y  2):V3«Z(3  -Z(2, 

120  W1»X(3)-X(1):W2«Y(3)-Y(1):W3«Z(3)-Z(1 

130  A(1 ,1 )*U1 :A(1 ,2)«U2:A(1 ,3)-U3 

140  A(2,l )=V1 :A(2.2)«V2:A(2.3)-V3 

150  Hl-V2*W3-V3*W2:H2«-(V1*W3-V3*W1 ):H3-V1*W2-V2*W1 

160  I 14-SQR (Hll 2+H2 ! 2+H3 ! 2 ) : IF  H4]0THEN  180 

170  PRINT  "*****  ERROR  - PTS  COLLI  NEAR  *****":GOTO  40 
1 80  HI -H*H1  /II4:II2»H2*H/H4:H3«I WI3/H4 
190  A(3,1)»H1 :A(3»2)«H2:A(3,3)*H3 
200  FOR  S-1T0  N 

210  FOR  T-S  TO  N:IF  A(T.S)[]0  THEN  230:NEXT  T 

220  PRINT  "***♦*  ERROR  PTS  COLLI  NEAR  *****":GOTO  40 
230  GOSUD  330 

240  A(S.S)«1/A(S.S):G0SUB  360 
250  FOR  T-l  TO  H:IF  T-S  THEN  270 
260  R--A(T.S):  A(T.S)-0:G0SUB  390 
270  NEXT  TzfiEXT  S 

280  FOR  S-fl  TO  1 STEP  -1:IF  V(S)«S  THEN  320 
290  FOR  J-l  TO  N 

300  B-A(J.S):  A(J,S)-A(J,V(S)):  A(J.V(S))-B 

310  NEXT  J 

320  NEXT  S:GOTO  420 

330  FOR  J-l  TO  N 

340  D-A(S.J):  A(S.J)«A(T.J):  A(T.J)«B 
350  NEXT  J:  V(S)»T:RETURN 
360  FOR  J-l  TO  N: IF  J-S  THEN  300 
370  A(S,J)»A(S,S)*A(S,J) 

380  NEXT  J: RETURN 

390  FOR  J-l  TO  N 

400  A(T,J)-A(T,J)+B*A(S,J) 

410  NEXT  J: RETURN 

420  1,(1 )=(X(2) !2-X(l ) !2+Y(2) !2-Y(l ) !2+Z(2) !2-Z(l ) !2)/2 
430  li  2 -(X(3) !2-X(2) !2+Y(3) !2-Y(2) !2+Z(3) !2-Z(2) !2)/2 
440  t*(3)-lil*X(l)+H2*Y(l )-HI3*Z(l ) 

450  X9-A(l ,1 )*D(1 )+A(l ,2)*B(2)+A(1 ,3)*B(3) 

460  YO-A  2.1  *D  1 +A  2.2  *B  2 +A  2.3  *B(3: 

470  Z,>A(3.1)*U(1)+A(3,2)*D(2)+A(3.3)*B(3) 
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Table  A- IV.  RCC  Progra*  Listing  (Continued) 


480  R9=SQR( (X(l )-X9) !2+{Y(l )-Y9) !2+(Z(] )-Z9) !2) 

490  SELECT  PRINT  215(90) 

500  PRINT  "INPUT  POINTS:" 

510  FOR  1=1  TO  3: PRINTUSING  520,I,X(I),Y(I),Z(I):NEXT  I 
520  '/POINT  f1  X-#####.####  Y -#####.####  Z—Mm.tM 

530  PRINT  "LENGTH  OF  HEIGHT  VECTOR-" ,H: PRINT 
540  PRINT  "THE  PARAMETERS  OF  THE  RCC:" 

550  PRINTUSING  560,X9,Y9,Z9 

560  ^CENTER  OF  CASE  XC «-#####.####  YC --Ml##.####  ZC— 

570  PRINTUSING  580,111  ,H2,H3 

580  WEIGHT  VECTOR  DX— #####.####  DY— ##*##.«##  DZ«-####*.###,4 
590  PRINTUSING  600,-111,-112,-113 

600  % or  dx--#####.##^#  m-MMJW  dz -#####.#### 

610  PRINTUSING  620,R9:PRINT  :PRINT  :PRINT 
620  %RADIUS  OF  CASE—#####.### 

630  SELECT  PRINT  005(64):  GOTO  40 


Table  A-V 


LINECIR  Program  Listing 


10  SELECT  PRINT  000(64) 

20  PRINT  "LINE  CIRCLE  PROGRAM" 

30  PRINT  "FINDS  INTERSECTION  OF  LINE  AND  CIRCLE" 

40  PRINT  "LINE  INPUT:" 

50  INPUT  "X.Y.DEL  X.DEL  Y«",X1 ,Y1 ,D1 ,D2 
60  PRINT  "CIRCLE  INPUT:" 

70  INPUT  "X,Y  OF  CENTER,  RADIUS-", X2.Y2.R 

30  IF  R]0T1IEN  100-.PRINT  "ERROR  RADIUS  [■  0,  TRY  AGAIN" 

90  INPUT  "RADIUS-", R:GOTO  00 
100  SELECT  PRINT  210(00) 

110  PRINTUSING  120.X1 ,Y1 

1 20%LINE  INPUT  X-#####.####  Y— #####.#### 

130  PRINTUSING  140.D1 ,D2 

140%  DEL  X-#####.####  DEL  Y— #####.#### 

100  PRINTUSING  160.X2.Y2.R 

160%CIRCLE  INPUT  XC »-#####.####  YC-#####.####  R •.#####.#### 

170  IF  Dl-OTHEN  300:S-D2/D1:Y3-Y1-S*X1:A-1.+S!2 

130  B-2.*S*(Y3-Y2)-2.*X2:C»X212+(Y3-Y2)!2-R!2 

190  D3»ti!2-4.*A*C:IF  D3[0.TIIEN  360: IF  D3-0.THEN  370 

200  X4»(-B+SQR(D3))/(?.*A):X5*(-B-SQR(D3))/(2.*A) 

210  SELECT  PRINT  210(00) 

220  Y4-S*X4+Y3:Y0-S*X0+Y3:PRINT  "INTERSECTION  POINTS" 

230  PRINTUSING  240,X4,Y4:PRINTUSING  240,X5,Y5 
Zm-MWJM  Y--WW.#### 

250  PRINT  : PRINT 

2C0  SELECT  PRINT  000(G4) 

270  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF" 

280  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 

200  STOP  :GOTO  40 
300  IF  D2-0.THEN  350 

310  D4«R!2-(X1-X2)!2:IF  D4-0T1IEN  330:IF  D4[0TUEN  360 
320  X4-X1 :X5-X1 :Y4-Y2+SQR(D4):Y5-Y2-SQR(D4):G0T0  340 
330  X4-X1 :Y4-Y2:PRINTUSING  300, X4.Y4: GOTO  260 
340  PRINT  "INTERSECTION  POINTS" :GOTO  230 

300  PRINT  "INPUT  ERROR,  SLOPE  UNDEFINED": PRINT  : PRINT  -.GOTO  260 
3C0  PRINT  "NO  INTERSECTION" .-PRINT  :PRINT  :GOTO  250 
370  X4«-U/(2.*A):Y4»S*X4+Y3:PRINTUSING  300,X4,Y4 
330XLINE  TANGENT  TO  CIRCLE  AT  X-#«##.####  Y— #####,#### 

390  GOTO  200 
400  STOP 
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Table  A-VI.  TANCIR  Program  Listing 


10  SELECT  PRINT  005 (G4) 

20  PRINT  "TANCIR  PROGRAM" 

30  PRINT  "FINOS  TANGENT  PTS  ON  A CIRCLE  FROM  A PT  OUTSIDE" 

40  INPUT  "INPUT  POINT  a,Y",X1,Y1 

50  INPUT  "INPUT  CIRCLE  CENTER  X.Y  AND  RADIUS: ".X2.Y2.R 

60  SELECT  PRINT  005(64) 

70  IF  R]OTHEH  P>0: PRINT  "ERROR  RADIUS[«0,  TRY  AGAIN": INPUT  "RADIUS»",R 
GOTO  70 

CO  SELECT  PRINT  215(00):PRINTU5II1G  90.Xl.vi 
i)OrROINT  INPUT  Y— «»#.#### 


100  PRINTUSING  110,X2,Y2,R 

1103CIRCLE  INPUT  KUn-»tm.MP  YCENT 


120  Aa(Xl-X2) !2+(Y1-Y2)!2 
130  IF  SPR(A  [RTIIEU  200 


140  IF  SOR 

A 

-RTliFN  300 

160  X3*X2+| 

!( 

m. 

*(X1-X2) 

ICO  X4-X2+ 

[( 

R!2] 

* X1-X2 

170  Y3-Y2+ 

[ 

R!2| 

* Y1-Y2 

ICO  Y4»Y2+ 

[( 

m) 

|*(Y1-Y2) 

-(R*(Y1-Y2)*SQR(A-R!2)))/A 
+ R*  Y1-Y2  *SQR  A-RI2)  )/A 
+(r*(X1-X2)*SQR(A-R!2)))/A 
i-(R*(X1-X2)*SOR(A-R!2 «))/A 


190  PRINTUSING  200.X3.Y3 

200r,TA.lGENT  PTS  X— #####.####  Y«-##k#P.#P4# 

210  PRINTUSING  220 .X4.Y4: PRINT  ‘.PRINT 

220 s x-#####.#w  y -tmt.mt 

230  SELECT  PRINT  005(64) 

240  PRINT  "TO  RUN  AGAIN.  KEY  CONTINUC.CR-LF" 

250  PRINT  "TO  STOP,  KEY  CLEAR, CR-LF" 

260  STOP 
270  GOTO  40 

200  PRINT  "POINT  IS  INSIDE  CIRCLE": PRINT  ‘.PRINT 


290  GOTO  230 

300  PRINT  "POINT  IS  ON  CIRCLE":PRINT  :PRINT 


310  GOTO  230 
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Table  A-VII.  PLANEINT  Program  Listing 


10  DIM  X(10),Y(10),Z(10),A(5,5) 

20  SELECT  D 

30  SELECT  PRINT  005(64) 

40  PRINT  -PLANEINT  PROGRAM'*: PRINT  "FINDS  INTERSECTION  POINT  OF  THREE  P 
LANES" 

50  PRINT  "EACH  PLANE  CAN  BE  INPUT  ONE  OF  THREE  WAYS:" 

GO  PRINT  " 1.  3 PTS  INDICATOR-1" 

70  PRINT  " 2.  PT,  ROT,  FB  INDICATOR-2" 

80  PRINT  " 3.  PLANE  COEFF  INDICATOR-3" 

90  FOR  1-1  TO  3: SELECT  PRINT  005(64) 

100  PRINT  "PLANE", I: INPUT  "INDICATOR-", J 

110  IF  J-1THEN  1 20 : I F J-2THEN  380:IF  J-3THEN  440:PRINT  "WRONG  INDICATO 
R TRY  AGAIN": GOTO  100 

120  PRINT  "INPUT  THE  3 PTS":FOR  K-1T0  3:PRINT  "POINT\K 
130  INPUT  "X,Y,Z«",X(K),Y(K),Z(K):NEXT  K 
140  Vl-X(l )-X(2)  :V2-Y  1)-Y(2):V3»Z(1)-Z(2) 

150  V4*X(1)-X(3):  V5«Y(1)-Y(3):V6«Z(1)-Z(3) 

150  NlaV2*V6-V3*V5:N2— V1*V6+V3*V4:N3-V1*V5-V2*V4 
170  T1-SQR(N1 !2+N2!2+N3!2) 

180  IF  Tl[jOTHrri  190:PRINT  "PTS  ARE  COLLINEAR,  TRY  AGAIN" :GOTO  120 
ICO  F-ARCSIN(II3/T1) 

200  IF  COS(F)[]OTHEN  210:16-1 :IF  N3[0THEN  250: 16-2: GOTO  250 

210  IF  N1/(T1*C05(F))]1THEN  220:IF  N1/(T1*C0S(F))[-1THEN  230:G0T0  240 

220  R»0:G0T0  250 

230  R*130:G0T0  250 

240  R-ARCCOS ( N1 / ( T1 *COS ( F ) ) ) : I F N2]-0THEN  250:R«3f»0-R 

250  A(I,1)«N1/T1  :A(I,2)-N2/T1 :A(I,3)-N3/T1:A(I,4)-A(I,1 )*X(1)+A(I,2)*Y 

(1)+A(I,3)*Z(1) 

260  SELECT  PRINT  215(80):PRINT  "INPUT  FOR  PLANE", I 
270  IF  IC]OTHEN  300 

230  PRINTUSING  290,X(1),Y(1),Z(1),R,F:G0T0  320 

2905  X-#####.####  Y-#####.####  Z-#####.####  ROT-###.###  F B— 

*##.### 

300  PRINTUSING  310,X(1),Y(1),Z(l),F 

nn^  X— #####.####  Y— #####.####  z-#####.####  ROT-*****  FB-###. 

320  FOR  L-2T0  3:PRINTUSING  330,X(L),Y(L),Z(L):NEXT  L 
3305  -#####.####  -#####.####  -#####.#### 

340  IF  IC-OTHEN  350: IF  I6-1TIIEN  345:PRINT  "*****SINCE  FB-90  ROT  NOT  UN 
IQUE":IG-0:G0T0  350 

345  16-0: PRINT  "*****SINCE  FB-90  ROT  NOT  UNIQUE" 

350  PRINTUSING  3C0,A(I,1),A(I,2),A(I,3),A(I,4) 

3C0‘»C0EFFICIENTS  A-#.#####  B-#.#####  C-#.#####  D— #####.#### 
370  GOTO  480 

300  PRINT  "INPUT  PT,  ROT,  F B": INPUT  "X,Y,Z,ROT,FB»",X(I),Y(I),Z(I),R, 
F 

390  A( 1 ,1 )-COS(F)*COS(R) :A( I,2)-C0S(F)*SIN(R) :A(I,3)»SIN(F) 
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Table  A-VII.  PLANEINT  Program  Listing  (Continued) 


400  A(I,4)W\(I,1)*X(I)+A(I,2)*Y(I)+A(I,3)*Z(I) 

410  SELECT  PRINT  21 5 (80): PRINT  “INPUT  FOR  PLANE", I 
420  PRINTUSING  290,X(I) ,Y(I),Z(I) ,R,F 
430  PRINTUSING  360,A(I.1),A(I,2),A(I,3),A(I.4):G0T0  480 
440  INPUT  "A,B,C,D •",A(I,1),A(I,2),A(I,3),A(1,4) 

450  SELECT  PRINT  21 5(80) : PRINT  "INPUT  FOR  PLANE", I 
460  PRINTUSING  470,A( I .1 ) ,A( I .2) ,A( I .3) ,A(I,4) :G0T0  480 
4702C0EFFICIENTS  A -#####.#####  B «-#####.#####  C ■-###?#.##««  D« 


wm.titm 


430  NEXT  I 

490  n»A(l,l)*(A(2.2)*A(3,3)-A(3,2)*A(2,3))-A(1,2)*(A(2,l)*A(3,3)-A(2,3 
)*A(3,1 ))+A(l , 3)*(A(2,1 )*A(3,2)-A(3,1 )*A(2.2) ) 

500  SELECT  PRINT  215(80) 

510  IF  DhOTHEN  520-.PRINT  "PLANES  DO  NOT  INTERSECT": GOTO  5C0 

520  Dl«A(1,4)*(A(2,2)*A(3,3)-A(3,2)*A(2,3))-A(l ,2)*(A(2,4)*A(3,3)-A(2, 

3)*A(3,4))+A(1,3)*(A(2,4)*A(3,2)-A(3,4)*A(2,2)) 

530  D2SA(1 ,1 )*(A(2t4)*A(3,3)-A(3,4)*A(2,3) )-A(l ,4)*(A(2,1 )*A(3,3)-A(2, 

3) *A(3,1))+A(1 ,3)*(A(2,1)*A(3,4)-A(3,1 )*A(2,4)) 

540  D3*A( 1 ,1 )*(A(2,2)*A(3,4)-A(3,2)*A(2,4))-A(1 ,2)*(A(2,1 )*A(3,4)-A(2, 

4) *A(3,1 ) )+A(l ,4)*(A(2,1 )*A(3,2)-A(3,1 }*A(2,2) ) 

550  X«D1/D:Y-D2/D:Z-D3/D 


560  PRINTUSING  F?n,X,Y,Z 

570i:illTERSECTI0N  POINT  X«-P###*.##P#  Y— ####*.####  Z —#####.#### 
580  PRINT  : PRINT  : SELECT  PRINT  005(64) 

590  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF" 

600  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 

610  STOP  :GOTO  90 


Table  A-VIII,  LINEPLAN  Program  Listing 


10  DIM  X(10),Y(10),Z(10) 

20  SELECT  D : SELECT  PRINT  005(G4) 

30  PRINT  "LINEPLAN  PROGRAM": PRINT  "FINDS  INTERSECTION  OF  A LINE  AND  PL 
ANE" 

40  PRINT  "LINE  IS  INPUT  BY  A POINT  AND  3 DELTAS" 

50  PRINT  "PLANE  CAN  BE  INPUT  ONE  OF  THREE  WAYS:" 

60  PRINT  " 1.3  PTS  INDICATOR-1" 

70  PRINT  " 2.  PT,  ROT,  FB  INDICATOR-2" 

00  PRINT  " 3.  PLANE  COEFF  INDICATOR-3" 

90  SELECT  PRINT  005(G4):PRINT  "INPUT  THE  LINE":INPUT  "X,Y,Z»",X(4),Y(4 
) *7.(4) 

100  INPUT  "DEL  X,  DEL  Y,  DEL  7-",Dl,D2,D3 
110  IF  D1[]0THEN  130:IF  D2[]0THEN  130:IF  D3[]0THEN  130 
120  PRINT  "ERROR  DELTAS  ALL  ZERO,  TRY  AGAIN": GOTO  90 
130  INPUT  "INDICATOR  FOR  PLANE  INPUT-", J 

140  IF  J-1THEN  150: IF  J-2THEN  410:IF  J-3THEN  470:PRINT  "WRONG  INDICATO 
R TRY  AGAIN": GOTO  100 

150  PRINT  "INPUT  THE  3 PTS": FOR  K-1T0  3: PRINT  "POINT", K 
160  INPUT  "X,Y,Z«",X(K),Y(K),Z(K):NEXT  K 
170  V1*X(1)-X(2)  :V2«Y(1)-Y  2):V3»Z(1)-Z(2) 

130  V4«X(1)-X(3):  V5-Y(l )-Y(3):V6-Z(l )-Z(3) 

1 90  N1 cV2*V6-V3*V5 : N2--V1 *VG+V3*V4 : N3-V1 *V5-V2*V4 
200  T1»SQR(N1!2+N2!2+N3!2) 

210  IF  T1C10THEN  220:PRINT  "PTS  ARE  COLLINEAR,  TRY  AGAIN": GOTO  150 
220  F-ARCSIN(N3/T1) 

230  IF  COS(F)[]OTHEM  240:16-1 :IF  N3[0THEN  230: 16-2: GOTO  230 

240  IF  N1/(T1*C0S(F))]1THEN  250:IF  Nl/(T1*C0S(F))[-1THEN  260:G0T0  270 

250  R-0:G0T0  230 

260  R-130:G0T0  230 

270  R-ARCC0S(N1/(T1*C0S(F))):IF  N2]»0THEN  280:R-3C0-R 
280  A-N1/T1  :B-N2/T1  :C-N3/T1  :D-A*X(1 )+B*Y(l  )*C*Z(1 ) 

290  SELECT  PRINT  215(30):PRINT  "INPUT  FOR  PLANE" 

300  IF  I6]0THEN  330 

310  PRIfITUSING  320,X(1),Y(1),Z(1),R,F:C0T0  350 

320-;  X— ?*#.####  Y-#####.####  Z-#####.####  ROT-###.###  F B— 

.<##.*## 

330  PRINTUSING  340, X(1 ),Y(1 ),Z(1 ),F 

340*  X-#####.####  Y -#####.####  z-#####.####  ROT-*****  FB-###. 

##«* 

350  FOR  L-2T0  3:PRINTUSING  3C0,X(L),Y(L),Z(L):NEXT  L 
360S  -####*.#*##  -#####.####  -#####.#### 

370  IF  I6-0TIIEN  380: IF  I6-1THEN  375:PRINT  "*****SINCE  FR-90  ROT  NOT  UN 
IOIIE":  16-0:  GOTO  330 

375  PRINT  "*****SINCC  FB-90  ROT  NOT  UNIQUE":  16-0 
380  PRIimiSING  390,A,B,C,0 

390' .COEFFICIENTS  A --#.#####  B--#.#####  C— #.#####  D--###4#.#### 
400  GOTO  520 
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Table  A-VIII.  LINEPLAN  Program  Listing  (Continued) 


410  PRINT  "INPUT  PT,  ROT,  F n":INPUT  "X,Y,Z,R0T,FB»",X(1),Y(1),Z(1),R, 
F 

420  A*COS(F)*COS(R):B«COS(F)*SIN(R):C«SIN(F) 

430  D«A*X(1)+B*Y(1)+C*Z(1) 

440  SELECT  PRINT  215(80) :PRINT  "INPUT  FOR  PLANE" 

450  PRINTUSINCi  320,X(1),Y(1),Z(1),R,F 
460  PRINTUSING  390,A,B,C,D:G0T0  520 
470  INPUT  "A  B C 0*"  A B C D 

480  IF  A[]OTHEN  490: f F*Bf ]OTI!EN  490: IF  C[]OTHEN  490:PRINT  "ERROR  ALL 

COEFFICIENTS  - 0,  TRY  AGAIN" :GOTO  470 

490  SELECT  PRINT  215(80):PRINT  "INPUT  FOR  PLANE" 

500  PRINTUSING  510,A,B,C,D 

510SC0EFFICIENTS  A«-#####.#####  B— #####.#####  C ■-#«##.##«*  D- 

520  SEI.ECT  PRINT  215(80):PRINT  SPRINT  "LINE  INPUT" 

630  PRINTUSIUG  540,X(4) ,Y(4),Z(4):PRINTUSING  550,01 ,D2,D3 
540 Y »-#####.##/»#  z— #####.#### 

550%DEL  X«-; DEL  Y«- #####.####  DEL  l—HMJW 

560  IF  A*D1+B*D2+C*D3[]0THEN  570-.PRINT  "LINE  PARALLEL  TO  PLANE,  NO  INT 

ERSECTION":GOTO  610 

570  S« ( D-A*X ( 4 ) -B*Y ( 4 ) -C*Z{ 4 ) ) / ( A*D1 +B*D2+C*D3) 

580  X-01 *S+X (4 ) : Y-02*S+Y(4 ) : Z»D3*S+Z (4 ) 

590  PRINT  : PRINTUSING  600,X,Y,Z 

GOO;;  INTERSECTION  POINT  X*-#####.####  Y -MW.####  Z«-#W#.###I 

610  PRINT  -.PRINT  :PRINT  : SELECT  PRINT  005(64) 

620  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF" 

630  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 

640  STOP  :GOTO  40 
650  STOP  :END 


t 

r t 
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Table  A-IX.  LINELINE  Program  Listing 


10  PRINT  "LINE  LINE  PROGRAM — FINDS  INTERSECTION  POINT  OF  TWO  LINES" 

20  PRINT 

30  PRINT  "LINE  1:" 

40  INPUT  "X.Y.DEL  X.UEL  Y»",X1  ,Y1  .D1  ,D2:PRINT  "LINE  2:":  INPUT  "X,Y,DE 
L X.DEL  Y»",X2,Y2,D3,D4 
50  SELECT  PRINT  215(00) 

GO  PRINT  "LINE  X Y DEL  X DEL  Y" 

70  J»1 : PRINTUSING  90,J,X1  ,Y1  ,D1 ,1)2 

00  J»2: PRINTUSING  90,J,X2,Y2,D3,D4 

90?  # -wmjw  -MMJM  -wtt.mt  -#####.#### 

ion  IF  D2C10.THEN  110: IF  D4[]0.THEN  1)0: IF  Y1-Y2THEN  240:G0TO  260 

110  IF  D1°0.THEN  210:t!l»D2/Dl  :D1»Y1-M1*X1 

120  IF  D3-0.THEN  250:N2«D4/D3:B2-Y2-M2*X2 

130  IF  M1412THOI  260:X»(D2-C1 )/(Ml-M2):Y»Ml*X+Bl 

140  PRINTUSING  150.X.Y 

1503 INTERSECTION  POINT  X—#####.####  Y— #####.#### 

1G0  PRINT  : PR I NT 

170  SELECT  PRItrr  005(G4) 

1 30  SELECT  PRINT  005(64): PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 

190  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF" :STOP 
200  GOTO  30 

210  IF  03*0. THEN  23U:X»X1  :Y«*04/U3*X+(Y2-D4/03*X2) 

220  GOTO  140 

230  IF  X1-X2THEN  240-.G0T0  2G0 

240  PRINT  "LINES  ARE  SANE": PRINT  :PRINT  :G0T0  100 

250  X»X2:Y411*X+B1:G0T0  <40 

2G0  PRINT  "LINES  00  NOT  INTERSECT": PR I NT  :PRINT  :G0T0  100 
270  STOP 


70 


Table  A-X.  RPARB  Program  Listing 


10  DIM  X(10),Y(10),Z(10) 

20  SELECT  D 

30  SELECT  PRINT  005(64) : PRINT  "RFARB  PROGRAM" 

40  PRINT  "FINDS  ARB8  DEFINED  BY  A POINT,  ROT  ANGLE," 

50  PRINT  "FB  ANGLE,  2 COORDINATES  OF  3 OTHER  POINTS," 

60  PRINT  "AND  A THICKNESS" 

70  INPUT  "X,Y,Z  OF  PT.  1,  ROT,  FB  ANGLES", X(1),Y(1),Z(1),R,F 
30  A=COS(F)*COS(R)  :G=COS(F)*SIN(R)  :C-SIN(F) 

90  D=A*X(1 )+B*Y(l )+C*Z(l ) 

100  FOR  I=2T0  4 
110  PRINT  "PT  NUM-M 

120  INPUT  "INPUT  INDICATOR  AND  2 KNOWN  COORDINATES", II  ,12,13,111  ,IJ2 

130  IF  IllOTHEN  190 

140  IF  I2J0THEN  220 

150  IF  I3=0THEN  490 

160  IF  C-OTHEN  500 

170  X(I)=I!1  :Y(I)=IJ2:Z(I)=(D-A*U1-B*U2)/C 

100  GOTO  240 

190  IF  A-OTHEN  520 

200  Y(I)-tll  :Z(I)-U2  :X(I)«(D-B*II1-C*U2)/A 
210  GOTO  240 

IF  R»OTHFN  5*10 

230  X(I)-U1  :Z( I)«U2  :Y(I)«(D-A*U1-C*U2)/B 
240  NEXT  I 

250  INPUT  "THICKNESS®", T 
260  FOR  I*5T0  3 
270  J-I-4 

230  X(I)SX(J)+(A*T) 

290  Y I -Y (J  + R*T) 

300  Z(I)®Z(J)+(C*T) 

310  NEXT  I 

320  DlaX(5)*A+Y(5)*B+Z(5)*C 
330  SELt:CT  PRINT  215(30) 

340  PRINT  "PT  X Y Z" 

350  FOR  1=1  TO  3 

300  PRINTUSING  300,I,X(I),Y(I),Z(I) 

370  NEXT  I 

-mw.mt  -#####.####  -#####.#### 

39 0 PRINT  "FACE  ARC  D ROT 

F 1." 

400  Kl=1234  :K2«5673 

410  PRINTUSING  d20,Kl ,A,B,C,D,R,F:PRINTUSING  420,K2,A,B,C,D1 ,R,F 

42 WM  -#.#####  -#.#####  -#.###44  -####.#####  -###.##  -###.#:? 

430  PRINT  "THICKNESS*", T 

440  PRINT  : PRINT 

450  SELECT  PRINT  005(64) 

460  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 


Table  A-X.  RFARB  Program  Listing  (Continued) 


470  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF" 

480  STOP  :G0T0  70 

490  PRINT  "ERROR,  NO  COORDINATE  TO  SOLVE  F0R":PRINT  :PRINT  :GOTO  110 
500  PRINT  "Z  NOT  UNIQUE  IN  THIS  PLANE": PRINT  ‘.PRINT 
510  GOTO  110 

520  PRINT  "X  NOT  UNIQUE  IN  THIS  PLANE"  : PRINT  .‘PRINT  :GOTO  110 
530  PRINT  "Y  NOT  UNIQUE  IN  THIS  PLANE"  .‘PRINT  : PRINT  .‘GOTO  110 
540  END 


Table  A-XI.  3PTARB  Program  Listing 


10  DIM  X(10),Y(10),Z(10) 

20  SELECT  D 

30  SELECT  PRINT  005(64) 

40  PRINT  "3PTARD  PROGRAM": PRINT  "FINDS  ARBfl  DEFINED  BY  3 POINTS,  2 COO 
RDINATES  OF" 

50  PRINT  "A  4TII  POINT,  AND  A THICKNESS" 

60  PRINT  "INPUT  3 POINTS" 

70  FOR  1=1  TO  3:PRINT  "POINT", I:  INPUT  "X»",X( I) : INPUT  "Y«",Y,(I): INPUT 
"Z=",Z(I):NEXT  I 

CO  V13X(1)-X(2):V2=Y(1)-Y(2):V33Z(1)-Z(2) 

90  V4=X(1 )-X(3):V5=Y(1)-Y(3):V63Z(l )-Z(3) 

100  H1=V2*V6-V3*V5:N23-V1*V6+V3*V4:N3«V1*VD-V2*V4 
110  T13S0R(N1!2+N2!2+N3!2) 

120  IF  T1P0THEN  130:PRINT  "PTS  COLLIMEAR":GOTO  390 
130  F=ARCSIN(N3/T1 ) 

140  IF  COS  F)[]OTIIEN  150:11=1  :IF  II3[0THEN  190: II  =2 .-GOTO  190 

150  IF  II1/(T1*C0S(F))]1THEN  160:IF  N1/(T1*C0S(F) )[-lTKEN  170:G0T0  100 

160  R=0:G0T0  190 


170  R31RO:GOTO  190 

ICO  R3ARCC0S(N1/(T1*C0S(F) ) ) : IF  N2]“0i'HEN  190:R*360-R 
1 90  A=N1 /T1 : B=N2/T1 : C=N3/T1 : D=A*X ( 1 ) +R*Y ( 1 ) +C*Z ( 1 ) 

200  INPUT  "INDICATOR  AND  2 KNOWN  COORDINATES3", 12,13, I4,U1 ,U2 
210  IF  I2]0THEN  240: IF  I3]0T1IEN  2G0:IF  I4[30TNEN  400 
220  IF  C=OTHEN  490:X(4)3U1 :Y(4)3U2:Z(4)3(D-A*I!1-B*U2)/C 
230  GOTO  270 

240  IF  A=OTHEN  500:Y(4)3U1 :Z(4)3II2:X(4)3(D-B*U1-C*U2)/A 
250  GOTO  270 

260  IF  B-OTIIEN  510:X(4)3U1 :Z(4)3H2:Y(4)3(D-A*U1-C*U2)/B 
270  INPUT  "THICKNESS3", T 

SCO  FOR  I3jTO  0:J3I-4:X(I)3X(J)+T*A:Y(I)3Y(J)+T*B 
290  ?:( I )»Z( J )+T*C : NEXT  I:U13A*X(5)+R*Y(5)+C*Z(5) 

300  SELECT  PRINT  215(30) 

310  PRINT  "PT  X Y 7" 


320  Fvin  I' 
330:0V*  ■ 

34-:  PRINT 


S1T0  3: PRINTUSING  330,I,X(I),Y(I),Z(I):NF.XT 

WWWW 


"FACE 


A 


I 

D 


ROT 


Fl>" 

350  K131234:K2B5670:IF  I1]0THEN  420 

360  FRINTUSING  370.K1 ,A,B,C,D,R,F:PRINTWSING  370,K2,A,B,C,D1 ,R,F 

3C0  PRINT  "THICKNESS3", T 
330  PRINT  : PRINT 


400  SELECT  PRINT  005(64) :PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  CR-LF" 
410  PRINT  "TO  STOP,  KEY  CLEAR  CR-LF": STOP  :GOTO  60 
420  SELECT  PRINT  215(30) 

430  PRIIITUSING  440.K1  ,A,B,C,D,F:PRIfiTUSING  440,K2,A,R,C,D1  ,F 

WWW  -*.WM  -#.#•?#,«#  ~PWWW*  *******  -i!r 


I ;•  :i 


Table  A-XI.  3PTARB  Program  Listing  (Continued) 


450  IF  I1-2THEN  460: I1«0:PRINT  "***  NOTE  SINCE  FB  IS  -90  ROT  IS  NOT  U 
NIQUE***":G0T0  470 

460  PRINT  "***  NOTE  SINCE  FB  IS  90  ROT  NOT  UNIQUE  ***":I1»0 

470  PRINT  "THICKNESS3", T:PRINT  ‘.PRINT  :G0T0  400 

480  PRINT  "INPUT  ERROR  NO  COORDINATE  TO  SOLVE  FOR" SPRINT  ‘.PRINT  :GOTO 

400 

490  PRINT  "7.  NOT  UNIQUE  IN  THIS  PLANE": PRINT  ‘.PRINT  :GOTO  400 

500  PRINT  "X  NOT  UNIQUE  IN  THIS  PLANE". -PRINT  : PRINT  :G0T0  400 

510  PRINT  "Y  NOT  UNIQUE  IN  THIS  PLANE": PRINT  :PRINT  :GOTO  400 


Table  A-XII.  NORMVEC  Program  Listing 


10  PRINT  "NORMAL  VECTOR  PROGRAM" 

20  PRINT  "FINDS  VECTOR  OF  LENGTH  LI  IN  DIRECTION  OF  GIVEN" 
30  PRINT  "VECTOR  V*(X,Y)  AND  VECTOR  OF  LENGTH  L2  NORMAL" 

40  PRINT  "TO  GIVEN  VECTOR  V" 

50  PRINT  "ORDER  OF  INPUT:  X,Y,L1,L2" 

CO  INPUT  "INPUT  X,Y,L1 ,L2",X,Y,L1 ,L2 
70  S=SQR(X!2+Y!2) 

80  X1=L1*X/S 

90  YlaLl*Y/S 

100  X2=-L2*Y/S 

110  Y2«L2*X/S 

120  SELECT  PRINT  215(80) 

130  PRINT  "INPUT  VECTOR:" 

140  PRINTUSING  230.S.X.Y 

150  PRINT  "VECTOR  IN  DIRECTION  OF  INPUT  VECTOR:" 

ICO  PRINTUSING  230, LI, XI, Y1 

170  PRINT  "VECTOR  PERPENDICULAR  TO  INPUT  VECTOR:" 

130  PRINTUSING  230.L2.X2.Y2 

IPO  PRINT  : PRINT 

200  SELECT  PRINT  005(04) 

210  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE.CR-LF" 

220  PRINT  "TO  STOP,  KEY  CLEAR, CR-LF" 

23rj.EIK!TH--####A.W  DEL  X«-####M#I  DEL  Y -MM.1W 
240  STOP 
250  GOTO  CO 
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Table  A-XIII.  PERPENV  Program  Listing 


10  SELECT  PRINT  005 {64) 

20  SELECT  D 

30  PRINT  "PERPENV  PROGRAM" 

40  PRINT  "FINDS  VECTOR  PERPENDICULAR  TO  2 GIVEN  VECTORS" 

50  INPUT  "VECTOR  1«".V1,V2,V3 
60  V4=S0R(V1!2+V2!2+V3!2):IF  V4[]0THEN  30 
70  PRINT  "*****  TRY  AGAIN  *****":GOTO  50 
30  INPUT  "VECTOR  2»",W1 ,H2,W3 
90  W4=SQR(W1!2+W2!2+W3!2):IF  W4[]0THEN  110 
100  PRINT  "*****  TRY  AGAIN  *****":GOTO  80 
110  INPUT  "LENGTH  OF  PERPENDICULAR  VECTOR®" ,L 
120  IF  LjOTHEN  140 

130  PRINT  "*****  TRY  AGAIN  *****":GOTO  110 
140  D3V1 *W1 +V2*W2+V3*W3 
150  T =ARCCOS ( D/ ( V4*W4 ) ) 

160  IJ1=V2*W3-V3*W2 
170  U2— (V1*W3-V3*W1) 

180  IJ3=V1  *W2-V2*U1 

190  U4=SQR(U1!2+U2!2+U3!2) 

200  IF  U4J0THEN  220 

210  PRINT  "ERROR- VECTORS  ARE  MULTIPLES": GOTO  370 
220  U1=L*U1/I]4:I12=L*U2/U4:U3“L*U3/U4 
230  SELECT  PRINT  215(90) 

240  PRINT  "INPUT  VECTORS:" 

250  PRINTUSING  2G0.Vl.V2.V3 

260  "DEL  X-#####.####  DEL  Y«- #####.###  DEL  Z-#####.#### 

270  PRINTUSING  2G0.U1 .1I2.W3 
280  PRINTUSING  290.V4.1/4 

290  LENGTHS  VECTOR  1 -####.####  VECTOR  2-####.#### 

300  PRINTUSING  310,T:PRINT 

310  WANGLE  RETWEEN  INPUT  VECTORS-###. ###  DEGREES 
320  PRINT  "VECTOR  PERPENDICULAR  TO  INPUT  VECTORS:" 

330  PRINTUSINO  350,U1,IJ2,U3 
340  PRNJTUSING  3C0.-U1  ,-U2,-l)3 

350  >;  D X-#####.####  DY— #####.####  DZ-#####.#### 

360  % OR  |)X— #####.####  DY— #####.####  DZ— #####.##:¥ 

370  PRINT  :PRINT  :PRINT  : SELECT  PRINT  005(64) :GOTO  50 
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Table  A-XIV.  AMTRACK  Program  Listing 


10  DIM  0(120), R$(12) 

20  FOR  1*1 TO  12: READ  R$(I):NEXT  I 

30  DATA  "RPP  ","D0X  "/ARB8"/RCC  ","FRNT  SLOPE  "/TRACK  BOT  "."REAR  SL 
OPE", "TRACK  TOP  ","TRCK  IDLER", "TRCK  DRIVE", "IDLER  DUM  "/DRIVE  DUM  " 
40  SELECT  PRINT  005(64) 

50  PRINT  "AMTRACK  PROGRAM" 

60  PRINT  "ADDS  TRACKS  TO  DOMESTIC  VEHICLES" 

70  PRINT  "ASSUMES  1.  IDLER  WHEEL  IS  IN  FRONT  OF  FIRST  ROAD  WHEEL" 

80  PRINT  " 2.  DRIVE  WHEEL  IS  TO  REAR  OF  LAST  ROAD  WHEEL" 

90  PRINT  " 3.  Z COORDINATE  OF  THE  ROAD  WHEELS  IS  EQUAL" 

100  INPUT  "X,Z  OF  FIRST  ROAD  WHEEL-",X1 ,Z1 
110  INPUT  "X,Z  OF  LAST  ROAD  WHEEL»",X3,Z3 

120  IF  Z1=Z3THEN  130:PRINT  "***  ERROR  - Z'S  OF  THE  ROAD  WHEELS  MUST  BE 
THE  SAME  ***":GOTO  100 

130  IF  X1]X3THEN  140:PRINT  "***  ERROR  - LAST  ROAD  WHEEL  TO  THE  FRONT  0 
F FIRST  ROAD  WHEEL  ***":GOTO  100 
140  INPUT  "RADIUS  OF  ROAD  WHEELS-\R1 

150  IF  R1]0THEN  160:PRINT  "***  ERROR  - RADIUS  [»  0 ***":GOTO  140 
160  INPUT  "X,Z  OF  IDLER  WHEEL-", X2.Z2 

170  IF  X2]X1 THEN  180:PRINT  "***  ERROR  - IDLER  WHEEL  IS  CONSIDERED  AS  I 
N FRONT  OF  1ST  ROAD  WHEEL  ***":GOTO  160 
180  INPUT  "RADIUS  OF  IDLER  WHEEL-\R2 

190  IF  R2]0THEN  200:PRINT  "***  ERROR  - RADIUS  [-  0 ***":GOTO  180 
200  INPUT  "X,Z  OF  DRIVE  WHEEL-", X4.Z4 

210  IF  X4[X3THEH  220:PRINT  "***  ERROR  - DRIVE  V/HEEL  IS  CONSIDERED  AS  T 
0 REAR  OF  LAST  ROAD  WHEEL  ***":GOTO  200 
220  INPUT  "RADIUS  OF  DRIVE  WHEEL-", R4 

230  IF  R4]0THEN  240:PRINT  "***  ERROR  - RADIUS  [-  0 ***":G0T0  220 
240  INPUT  "YMIN,  YMAX  OF  TRACK-", Y,Y1 

250  IF  V[Y1THEN  260:PRINT  "***  ERROR  - YMIN  ]-  YMAX  ***H:GOTO  240 
260  INPUT  "THICKNESS  OF  TRACK-", T 

270  IF  T]0TIIEM  2C0:PRINT  "***  ERROR  - TRACK  THICKNESS  MUST  BE  ]0  ***": 
GOTO  260 
200  D7-4 

290  X-Xl :Z-Z1 :R*R1 :K9-0:D-1 
300  M-X2:N-Z2:0-R2 
310  A1=X:B1»Z-R 
320  A1-A1+D7 

330  nosun  '01 

340  GOSUB  '02 

350  IF  D3]-0THEN  320 

360  C3»(B2-B1 )/SQR((Bl-B2) !2+(Al-A2) !2) 

370  C4-(A2-A1 )/SQR( (B1-B2) !2+(Al-A2) !2) 

380  C5-D*C3:C6— C4*D 

390  B4-D1 -T: A4- (A1 +(C5*T) ) -( ( (B1 +(T*C6) ) -B4)*C4/C3) 

400  0(1+K9)-A4:0(2+K9)-Y:0(3+K9)-B4 

41 0 0 ( 4+K9 ) - A2+( C5*T ) -A4 : 0 ( 5+K9 ) -0 : 0( 6+K9 ) -B2+( C6*T ) -B4 
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Table  A-XIV.  AMTRACK  Program  Listing  (Continued) 


420  0(7+K9)*-C5*T:0(8+K9)«0:0(9+K9)B-C6*T 
430  0(1 0+K9 ) «0 : 0(  1 1 +K9 ) ■ Y1  - Y : 0 ( 1 2+K9 ) -0 
440  IF  K9]0THEN  470 

450  0(26)«A4:0(27)«Y:0(28)»Y1 :0(29)-D1-T:0(30)«Rl 
460  K9«12:X=X3:2»Z3:D--1  :M-X4:tl-74:0»R4:D7— 07:R0T0 
470  0(25)=A4 
480 


310 


C1=(0(4)+0(1)-X2)/SQR((0(4)40(1)-X2)!2+(0(6)+0(3)-Z2)!2) 
C2*(0(G)+0(3)-Z2)/SqR((0(4)+0(l  )-X2)  !2+(0(r.)-H)(3)-7.2)  !21 


490  C2 

500  X=X2:R=R2:Z«Z2:K9«0 
510  Of 34+K9)»X+( ( R+T+300)*C1 ) 
520  0(35+K9) 

530 
540 

550  0(33+K9 
560  0(37+K9 
570  0 (40+K9 
530  0 43+K9 
590  0(46+K9 
600  0(49+K9) 

61 0 0(52+K9) 


Yl+lfl 

0(36+K9)-Z+((R+T+300)*C2) 

0 ( 31 +K9 ) -0 ( 34+K9 ) : 0 ( 32+K9 ) »0 ( 35+K9 ) 

“Z+R+T+50 

■X : 0 ( 38+K9 ) -0( 35+K9 ) : 0 ( 39+K9 ) -Z 
=:<:0(41+K0)»0(35+K9):0(42+K9)-0( 
0 ( 31 +K9 ) : 0 ( 44+K9 ) - Y-l 0 : 0( 45+K9  J 
0(34+K9) :0(47+K9)“Y-10:0(48+K9j 
0( 37+K9 ) : 0 50+K9  « Y-l 0 :0  51 +K9 
0( 40+K9 ) : 0( 53+K9 j ■ Y-l  0: 0( 54+K9J 


620  IF  K910THEH  G60 : K9-24 : X-X4 : Z-Z4 s R-R4 


33+K9) 

»0(33+K9' 

•0(36+K9 

»0(39+K9 

•0(42+K9| 


630  c!*(9?K,)+0P3)-^)/sqr{(o(16)+o(13)-x);2+(o(ib}+o( 

[15)- 


13)-X 

Z)/SQR((0(16)+0(13)-X)!2+(0(18)+0( 


640  C2«(0(18)+0( 

650  00T0  510 
660  K9=0:X=X2:Z»Z2:R«R2+T 
670  0(79+K9)"X:0(80+K9)*Y:0(31 +K9)*Z 
680  0 ( 82+K9 ) -0 : 0 ( 33+K9 ) «Y1 - Y : 0 ( 84+K9)»0 : 0( 85+K9 ) -R 
IF  K9]0THEH  700 : K9-7 : X-X4 : Z-Z4 : R-R4+T : GOTO  670 
IF  Z2+R2[]Z4+R4  THEN  730 


;i5)-Z)'2> 

;15)-Z)!2) 


690 

700 


710  0(94)-X2+10:0(93)«X4-10:0(95)»Y:0(96)-Yl 
720  O(97)=Z2+R2:0(98)-O(97)+T:L7-98:r,0T0  790 
730  C1-(X4-X2)/SQR(  X2-X4) !2+(Z2+R2-Z4-R4) !2) 

740  C2»(Z4+R4-Z2-R2  /SQR(  X2-X4)!2+(Z2+R2-Z4-R4)!2) 

750  0(93)-X2-(10*Cl):0(94)-Yl:0(95)»Z2+R2-(10*C2) 

760  0 96WX4-X2)+(20*C1):0(97)»0:0(98)»(Z4+R4-Z2-R2)+(20*C2) 

770  0(99)=T*C2:0(100)»0:0(101)-T*C1 

780  0(102)-0:0(103)»-(Y1-Y):0(104)-0:L7»104 

790  K9-0:X-X2:Z»Z2:R-R2 

COO  0(L7+K9+1 )»X:0(L7+K9+2)-Y:0(l7+K9+3)»Z 

810  0( L 7+K9+4 ) *0 : 0( L 7+K9+5 ) ■ Y1 - Y : 0( L 7+K9+6 ) "0 

320  0(L7+K9+7)-R 

330  IF  K9]0TflEfl  850 

840  K9*  7 : X * XI : Z*Z4 : R*R4 : GOTO  800 

350  SFLECT  PRINT  215(100):I-1 tPRINTUSING  860,1 ,R$(2) .0(1 ) .0(2) ,0( 3) ,0( 
4) »0(5),0(C),R5(5) 

860  r»f if  ?m  -mw.mt  -mw.tm 
-wm.w*  *mmm 


-wm.tm  -mw.tm  -ww.w 
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Table  A-XIV.  AMTRACK  Program  Listing  (Continued) 


1 


i 


S70  PRINTUSIMG  380,0(7) .0(3) ,0(9) ,0(10) ,0(11 ) ,0(12) 

080  % -Mm.m#  -#####.####  -MW.mt  -wm.tm  -mw.M 

890  Ia2: PRINTUSIMG  860, I,R$(2), 0(13) , 0(14), 0(15), 0(10), 0(17), 0(10), R$( 
7) 

900  PRINTUSING  080, 0(1 9), 0(20) ,0(21 ) ,0(22) ,0(23) ,0(24) 

910  I a3: PRINTUSIMG  0G0,I,R$(1 ) ,0(25) ,0(26) ,0(27) ,0(20) ,0(29) ,0(30) ,R$( 
6) 

920  Ia4:PRINTUSING  860,I,R$(3) ,0(31 ) ,0(32) ,0(33) ,0(34) ,0(35) ,0(3G) ,R$( 

ID 

930  PRINTUSIMG  880,0(37) ,0(38), 0(39), 0(40) ,0(41), 0(42 
940  PRINTUSING  000,0(43), 0 44  ,0  45  ,0  46), 0 47) ,0  48 ) 

950  PRINTUSING  330,0(49), 0(50  ,0  51), 0(52  ,0  53  ,0(54) 

960  1=5 : PR I flTUS I MG  060, 1, R$(3) ,0(55) ,0( 56) . 0(57), 0(50) ,0(59), 0(60), R$( 

12) 

970  PRINTUSING  080,0(61 ), 0(62 ),0(G3), 0(64) ,0(65), 0(66) 

900  PRINTUSIMG  380,0  67) ,0  68  ,0  69  ,0  70  ,0  71 ) ,0(72 
990  PRINTUSING  000,0(73), 0(74). 0(75), 0(76). 0(77). 0(70) 

1000  I *6: PRINTUSING  060, I,R$(4),0(79),0(80), 0(01), 0(02), 0(03), 0(84), R$ 

(9) 

1010  PRINTUSING  080,0(05) 

1020  Ia7:PRINTUSING  060, I, R$(4), 0(06), 0(07) ,0(00) ,0(09), 0(90) ,0(91 ) ,R$ 

(10) 

1030  PRINTUSIMG  080,0(92) 

1040  1=0: IF  7.2+R2aZ4+R4  THEN  1070 

1050  PRINTUSIMG  060,I.R$(2),0(93), 0(94) ,0(95), 0(96), 0(97), 0(90) ,R$(8) 
1060  PRINTUSING  000,0(99), 0(100), 0(101 ) ,0(102) ,0(103) ,0(104) :G0T0  1000 
1070  PRINTUSIMG  3G0,I,R$(1),0(93),0(94  ,0(95), 0(96), 0(97), 0(9P),R$(0) 
liJOO  I =9: PRINTUSING  060,1 ,R$(4) ,0(L7+1 ) ,0(L7+2) ,0(L7+3) ,0(L7+4) ,0(L7+5 
) ,0(L7+6),R$(11 ) 

1090  PRINTUSIMG  080,0(L7+7) 

1100  I=10:PRINTUSING  060,I,RS(4),0(L7+0) ,0(L7+9) , 0(17+10), 0(L7+11),0(L 
7+1 2), 0(17+13) ,R$(12) 

1110  PRINTUSIMG  000,U(L7+14) 

1120  PRINT  : PR I NT  :PRINT  "REGION  TABLE" 

1130  PRINTUSING  1140,1 ,1 ,-4,0,R$(5) 

1140  % w -w  -m  -###  #*####### 

1150  PRINTUSING  1140,2,2,-5,0,R$(7) 

1160  PRINTUSIMG  1140,3,3,-1 ,-2,R$(6) 

1170  PRINTUSING  1140,4,6,4,-9, R$(9) 

1180  PRiriTUSIfIG  1140,5,7, 5, -10,RS(10) 

1190  PRINTUSIMG  1140,G,0,-4,-5,R$(0) 

1200  GOTO  40 

1210  DEFRPOl 

1220  K2a(Al  -11)  !2+(Dl-fl)  !2 

1230  Glafl+(((Q!2)*(A1  -II))-(Q*(B1-N)*SQR(K2-Q!2)))/K2 
1 240  G2afl+(  ( (0!? )*(A1  -f I)  )+(Q*(Bl -N)*SQR(IC2-Q!2) ) )/K2 
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Table  A-XIV.  AMTRACK  Program  Listing  (Continued) 


1250  fi3-fl+(((n!2)*(C1-N))+(Q*(Al-H)*SQR(K2-Q!2)))/K2 
1260  r4«!l+i((Q!2)*(B1-H))-(Q*(A1-M)*SQR(K2-Q!2)))/K2 

1270  IF  G4[G3THEH  1280:B2»G3:A2«G1 :G0T0  1290 

12B0  B2«G4:A2«G2 

1290  RETURN 

1300  DEFFN'02 

1310  D1-A2-A1 :D2=B2-01 

1320  S«D2/D1:Y7»B2-S*A2:A»1+S!2 

1330  B»2*S*(Y7-Z)-2*X:C»X!2+(Y7-Z)!2-R!2 

1340  D3*F!!2-4*A*C 

1350  RETURN 
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Table  A-XV.  SOLIDROT  Program  Listing 


t > 
♦ 


v 

y 

f 

h 


w 


» t 


10  Din  A$(20),D(25),B(4,6),S(25),R(25),W(6) 

20  PRINT  "SOLIDROT  PROGRAM" 

30  PRINT  "ROTATES  IN  THE  XY,  XZ,  YZ  PLANES  SELECTED": PRINT  "COH-GEOM  S 
OLIOS  ABOUT  ANY  POINT" 

40  PRINT  "POSITIVE  ROTATION  IS  FROM  POSITIVE  AXIS  TO  POSITIVE  AXIS" 

50  FOR  1=1  TO  20 : READ  A$(I):NEXT  I 

GO  DATA  "RPP  "."BOX  "."RAW  ","SPH  Y'RCC  \"REC  Y’TRC  "."TEC  "."TOR  " 
,"ELL1","ELL  "."ARnO"."ARD7","ARRG"."ARB5"."ARB4",HBLK"."XY%"XZ"."YZ" 
70  SELECT  D 

00  SELECT  PRINT  005(G4) : INPUT  "PLANE  OF  ROTATION (1 -XY,2=XZ ,3=Y7 ) " .P 

90  IF  P[1THEN  1 00 : 1 F P]3TI!EN  100:GOTO  110 

100  PRINT  "*****TRY  AGA I tl*****" : GOTO  00 

110  INPUT  "ANGLE  OF  ROTATION(DEG)-" ,A 

120  INPUT  "X.Y.Z  OF  PT  TO  ROTATE  SOLID  ABOUT-", X5.Y5.Z5 

130  PRINT  "THE  FOLLOWING  SOLID  TYPES  ARE  AVAILABLE" 

140  PRINT  TAB(4) ;"RPP(1 ) , B0X(2),  RAW(3),  SPH(4).  RCC(5),  REC(G) , TPC( 
7)" 

15.)  PRINT  TAP (4); "TEC (3),  TOR(P),  ELL1(10),  ELL(II).  ARrB(l2).  AP.B7(13 
)" 

1G0  PRINT  TAB(4) ;"ARBG(14),  ARB5(15),  ARB4(16)" 

170  INPUT  "YOUR  SOLID  TYPE  IS",T 

180  IF  T[OTIIEN  190:011  T GOTO  210,340.400,550.500,020,600,730,790.330,0 

70,910 ,9GQ ,930 ,1 000 , 1 020 

190  PRINT  "******TRY  AGAIN******" .‘GOTO  170 

290  REM  ***SOLID  TYPE  IS  RPP*** 

210  INPUT  "XI1III,Xf!AX,YMIN,^lAX,ZMIN,ZnAX-",D(l),D(2),D(3),D(4),D(5),D( 
G) 

220  FOR  1*1  TO  5 STEP  2:IF  D(I)]D(I+1  )THEN  230:NF.XT  I:GOTO  2G0 

230  IF  1*1 THEII  240: IF  I-3THEN  250:PRIMT  "ERROR  ZHIH]ZMAX,  TRY  AGAIN": 

GOTO  210 

240  PRINT  "ERROR  XMINlXMAX,  TRY  AGAIN" :GOTO  210 
250  PRINT  "ERROR  YniNjYHAX,  TRY  AGAIN" :GOTO  210 
2G0  FOR  1=1  TO  6:l’(I)«n(I):NEXT  I 
270  R2*D(?,):K3-D(5) 

230  I!1=D(2)-[)(1):W1-D(4)-D(3):D1-D(G)-D(5) 

290  0(2K2:l)(3  -K3 
300  0(4  =lll:n(5)=0:D(G)»0 
310  D(7)=0:D(,3}-W1 : D ( 9 ) *0 

320  n(10)B0:D(ll)«0:D(12)-Dl:FG-fi:F«10:F7»3:GOTO  390 
330  REM  ***SOLID  TYPE  IS  BOX*** 

34ii  INPUT  "X.Y.Z  OF  VERTEX-", D(1),D(2),D( 3) 

350  PRINT  "NOTE — THE  ORDER  OF  FOLLOWING  VECTORS  MAY  VARY" 

330  INPUT  "HEIGHT  VECTOR-", D(4),D(5),D(G) 

370  INPUT  "WIDTH  VECTOR-", D(7) ,D(3) ,D(9) 

330  INPUT  "DEPTH  VECT0R=\D(10),D{11 ),D(12):F-10:F7-3:FG»12 
39(1  FOR  1=1  TO  25:S(I)-D(I):NEXT  I 

400  FOR  1=1  TO  FSTEP  3:IF  I]F7TNEN  410:S(I)-D(I)-X5:S(I+1)-D(I+1)-Y5:S( 
1+2)  =0(1+2) -7.5 


I 
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Table  A-XV.  SOLIDROT  Program  Listing  (Continued) 


410  IF  P-1THEN  420: IF  P-2THEN  430:IF  P« 
RY  AGAIN" ‘.GOTO  HO 


3THEN  440: PRINT  "INPUT  ERROR,  T 


420  GOSUB 
430  GOSUB 
440  GOSUB 


,S(I+2) 
.S  1+2 
.S  1+2)) 


e 

. e 
e 

, e 


GOTO  450 
GOTO  450 


450  IF  I]F7TUEN  460:R(I)-R(I)+X5:R(I+l)«R(I+l)+Y5:R(I+2)«R(I+2)+Z 
460  NEXT  I :G0T0  1030 

470  REM  ***SOLID  TYPE  IS  Ml*** 

4C0  INPUT  "X.Y.Z  OF  VERTEX-", 0(1), 0(2), 0(3) 

490  PRINT  "NOTE-HEIGHT  AND  WIDTH  VECTORS  MAY  BE  INTERCHANGED" 

500  INPUT  "HEIGHT  VECTOR-", D(4) ,D(5),D(6) 

510  INPUT  "WIDTH  VECTOR-" ,D(7) ,D(R) ,D(9) 

520  INPUT  "DEPTH  VECT0R-",D(10),D(11),D(12) 

530  F-10:F6-12:F7-3:G0T0  390 
540  REM 
550  INPUT 
560  INPUT 
570  REM 
500  INPUT 
590  INPUT 
GOO  INPUT 
610  REII 
620  INPUT 
630  INPUT 
640  INPUT 


"X.Y.Z  OF  CENTER-", 


***SOLID  TYPE  IS  SPH*** 

D(1),D(2),D(3) 

"RADIUS-", D(4):R(4)-D(4):F-1:F6-4:F7-3:G0T0  390 
***SOLID  TYPE  IS  RCC*** 

"X.Y.Z  OF  VERTEX-", D(1),D(2),D(3) 

"HEIGHT  VECTOR-", D(4),D{5),D(6) 

"RADIUS  OF  DASE-",R(7):F-4:F6-7:F7-3:D(7)-R(7):G0T0  390 
***SOLID  TYPE  IS  REC*** 

"X.Y.Z  OF  VERTEX-", D(l) ,D(2),D(3) 

HEIGHT  VECTOR-" ,D(4) ,D(5) ,D(C) 


"VECTOR  DEFINING  SEMI-MAJOR  AXIS«".D(7),D(8) ,D(9) 

650  INPUT  "VECTOR  DEFINING  SEMI-MINOR  AXIS-" ,D(1 0) ,D( 11 ) ,D(1 2) 
660  F-10:F6-12:F7-3:G0T0  390 


670  REM  ***SOLID  TYPE  IS  TRC*** 

600  INPUT  "X.Y.Z  OF  VERTEX-", D(1).D(2),D(3) 

690  INPUT  "HEIGHT  VECTOR-", D(4),D(5),D(G) 

700  INPUT  "LARGE  RADIUS-", R(7) : INPUT  "SMALL  RADIUS-",R(0) 

710  F-4:FC-B:F7-3:D(7)-R(7):D(0)-R(8):G0T0  390 
720  RD1  ***SOLID  TYPE  IS  TEC*** 

730  INPUT  "X.Y.Z  OF  VERTEX-", 0(1 ) ,D(2) -D(3) 

740  INPUT  "HEIGHT  VECTOR-" ,D(4) ,D(5) ,D(C) 

750  INPUT  "VECTOR  DEFINING  SEMI -MAJOR  AXIS-",D(7) ,D(8) ,D(9) 

760  INPUT  "VECTOR  DEFINING  SEMI-MINOR  AXIS-",D(10),D(11),D(12) 

770  INPUT  "RATIO-" ,R(13):F-10:F6-13:F7-3:D(13)-R(13):G0T0  390 
700  REM  ***SOLID  TYPE  IS  TOR*** 

790  INPUT  "X.Y.Z  OF  VERTEX-". D(1).D(2),D{3) 

OuO  INPUT  "NORJIAL  VECT0R-",D(4) ,D(5),D(6) 

010  INPUT  "R1-",R(7):IHPUT  "R2-",R(n):F-4:FG-n:F7-3:D(7)«R(7):D(0)«R(n 
):GOTO  390 

020  REM  ***SOLID  TYPE  IS  ELLl *** 

030  INPUT  "X.Y.Z  OF  VERTEX-",D(1 ) ,D(2) ,D(3) 

040  INPUT  "VECTOR  DEFINING  SEMI-MAJOR  AXIS-",D(4),D(5),D(6) 

050  INPUT  "RADHIS»",R(7):F»4:FG-7:F7»3:D(7)»R(7):G0T0  3<V) 


m 


Table  A-XV.  SOLIDROT  Program  Listing  (Continued) 


860  REM  ***S0LID  TYPE  IS  ELL*** 

870  INPUT  "X.Y.Z  OF  FOCI  1-".D(1).D(2).D(3) 

880  INPUT  "X.Y.Z  OF  FOCI  2a",D(4),D  5),D(G) 

890  INPUT  "LENGTH  OF  MAJOR  AXISa",R(7) :Fa4:FG»7:F7a6:D(7)aR(7) :GOTO  39 
0 

900  REM  ***SOLID  TYPE  IS  ARBB*** 

910  FI =22 : F=22 : FG=24 : F7aF6 

920  J«0:F0R  I=1T0  F1STEP  3:JaJ+l  : PRINT  "POINT*,.! 

930  INPUT  ''X,Y,Za",n(I),D(I+l),D(I+2):NEXT  I 
940  GOTO  390 

950  REM  ***SOLID  TYPE  IS  ARB 7*** 

9C0  Fa19:Fl=19:F6a21  :F7aF6:G0T0  920 

970  REM  ***SOLID  TYPE  IS  ARB6*** 

930  F=16:Fl=lG:FG=13:F7aFG:G0T0  920 

990  REM  ***SOL ID  TYPE  IS  AR85*** 

1000  F=1 3 : FI  =1 3 : F6al  5 : F7aF6 :GOTO  920 

1010  REM  ***SOLID  TYPE  IS  ARB4*** 

1020  F=10:Fl=10:FG=12:F7aFG:G0T0  920 

1030  SELECT  PRINT  215(85):PRINTUSING  1040,A,A$(P+17) 

1040  WANGLE  OF  ROTATIOH—I**#.##*#  DEG  IN  THE  PLANE 
1050  PRIIITUSING  1060,X5,Y5,Z5:G6aFG 

10G0  -;PT  AROllill)  WHICH  SOLID  HAS  ROTATED  Xa-fW#.M»#  Ya-’?iW.iW  Za- 


1070  K=0 : FOR  I-lTO  4: FOR  JalTO  G:KaK+l :D(I,J)aD(K):NEXT  JsNFXT  I 
1030  IF  THITHER  1090:FOR  I-lTO  6:0(1 , 1 )«!.'( I) .‘NEXT  I 
1090  PRINT  "INPUT  SOLID":IF  T[14THE!I  11 00: PRIIITUSING  11 10,A$(T) ,0(1 ,1 ) 
,R( 1 ,2) .0(1,3) ,B(1 ,4) :GOTO  1200 

1100  PRIilTUSING  ino,A$(T),B(l,l).B(l,2),B(l,3),R(l ,4),B(1 ,5),B(1  ,G) 
mo  %mt  -m«jM  -mt.wt  -tw.tw  -tw.w*  -« 

ii2o  % -mmm  -wt.mt  -m 


M.Wf 

1130  GGaGG-G:IF  GG[6TIIEII  1150:F0R  I-2T0  4:PRINTUSING  1120,B(I,1),R(I,2 
) .6(1,8)  ,B(I,4),l>(I,r»),n(I,f>) 

1140  G6»GG-G:IF  GufGTJIEN  llGO:fOT  IsGOTO  1200 
1150  1=1 


11C0  IF  GO-ITHEN  1170: IF  GCa2TNCM  1180: IF  GG-3THEN  1190: IF  GG-OTHFN  12 

00 

1170  PRINTUSING  1120,B(I+l,l):r,0T0  1200 

11.0  PRINTUSING  1120,3(1+1,1  ,B(I+1,2):G0T0  1200 

1190  PRINTUSIMG  112O,B(I+l,l),n(I+l,2),B(I+l,3):GOT0  1200 

1200  IF  E9-0TIIEN  1 21 0:E9a0: PRINT  :PRINT  :GOTO  12G0 

1210  PRINT  "ROTATED  SOLID" 

1220  IF  T[]4Tl!OI  1230:PRi;iTUS!flG  1110,A$(T),R(1),R(2),R(3),R(4):PRINT 
: PRINT  :GOTO  12GD 

1230  l>0:F0K  I=1T0  4: FOR  J-1TO  C:K«K+1  :R(I,J)-R(K):iOT  J:f!EXT  I 
12*'0  IF  T[]lTfiOI  1250:T»2:F;>12 


Table  A- XV.  SOLIDROT  Program  Listing  (Continued) 


12f>0  EO-1  :Q6-FG:G0T0  1100 
1260  SELECT  PRINT  005(G4) 

1270  PRINT  "TO  RUN  AGAIN,  KEY  CONTINUE  EXECUTE" 
1200  PRINT  "TO  STOP,  KEY  CLEAR  EXECUTE" 

1290  STOP  : GOTO  00 
1300  DEFFN'Ql (SI ,S2,S3) 

1310  R(I)»$1*C0S(A)-S2*SIN(A) 

1320  R 1+1  »S1*5IN(A)+S2*C0S{A) 

1330  R(I+2)»S3 
1.340  RETURN 

1350  DEFFN'02(S4,S5,Sf>) 

1360  R(I)»S4*C0S(A)-SG*SII!(A) 

1370  r; j T+i  )«r>n 

1 300  R( I+2)«S4*SIi!(A)+S0*C0S(A) 

1.390  RETURN 

1400  nEFFN'03(S7,SfJ,S9) 

1410  R(I)“S7 

1420  R 1+1 )«S0*C0S{A)-S9*SIN(A) 

1430  R( I+2)-S0*SIN(A)+S9*C0S(A) 

1440  RETURN 


Table  A-XVI.  PLOTSOL  Program  Listing 


10  Dll!  A(120),B(120),C(120),D(120),I!(120),F(8) ,G(8) ,M(8) ,KG(30) 

20  DIM  11(240), 0(240), P(240) .0(240) ,R(240) 

30  DIM  S{  1 20) ,T{  120) ,IJ(  120) ,V(  120) ,X(5),Y(5),Z(G) 

40  PRINT  "PLOTSOL  PROGRAM" 

50  PRINT  "PLOTS  SELECTED  COM-GEOM  SOLIDS  AT  ANY  ASPECT" 

60  PRINT  "HAS  SCALE  ADJUST  CAPABILITIES" 

70  PRINT  "SUM  OF  RCC'S  AND  TRC'S  MUST  BE  [ = 5" 

CO  PRINT  "MAX  HUM  OF  SPHERES  = 5" 

00  SELECT  I; 

100  INPUT  "HUM  OF  SOLIDS  TO  PL0T«",K9:G9=K9 

110  IF  K9]0THEM  120:PRIHT  "***  TRY  AGAIN  ***": PRINT  IIEX(07)tG0T0  100 
120  INPUT  "AZ,  EL  FOR  THIS  VIEtl-",A,E 

130  X0«0:X9-n:Z«>0:Z9-0:Kl*l:K2-0:K3«0:KS(l)»l:L9»0:A9*0:E9-l:S9=" 

140  FOR  I-E9T0  GO 

150  PRINT  "YOU  ARE  ON  SOLID'M 

160  PRINT  i;EX(07) 

170  PRINT  "THE  FOLLOWING  SOLIDS  ARE  AVAILABLE" 

100  PRINT  TAB(4);"RPP{1)  R0X(2)  RAlf(3)  ARR8(4)" 

11jO  PRINT  TAB (4  ;"ARC7(5)  ARBC(O)  ARB5(7)  ARB4(0)" 

200  PRINT  TAP.(4);"RCC(9)  TRC(IO)  SPH(ll)" 

210  INPUT  "SOLID  TYPE  IS",T 

220  OH  TGOTO  330,540,800,240,720,740,760,780,1000,1310,1520 
220  PRINT  "***  TRY  AGAIN  ***":GOTO  210 
240  REM  ***  A R B 3 *** 

250  K2»K1+7:IF  K2[-120THEH  260:PRIHT  "HO  MORE  ROOM  FOR  ARBS' S":K2»K2-7 
:GOTO  1600 

260  J-Q:FOR  IC-K1T0  K2:J»J+1  .‘PRINT  "POINT", J 

270  INPUT  "X,Y,Z»",A(K),ff(K),B(K):NEXT  K 

200  INPUT  "INPUT  OK  YES-1  H0«0",P 

230  IF  P-ITIIEN  300: GOTO  260 

300  PRINT  "***  THINKING... THINKING  ***" 

31  FOR  K-K1T0  K2 

32n  C(K)-(i!  K)*COS(A) )+(A(K)*SIH(A) ) 

330  I)(i;)=(l-(K)*CbS(E))+(A(K)*COS(A)*SIH(E))-(H(K)*SIN(A)*SIfi(E)) 

340  IF  i;r]lTHEN  350:IF  X8[]X9TIIEH  350:XC»C(K):X9»C(K):ZC«D(K):79*D(!:): 
GOTO  360 

350  GOSUB  ,07(C(K),D(K)) 

3C0  NEXT  K 

370  K1 -K2+1  :K3-K2:J«I-LS-S9+1  :KG(J)»K1  .-GOTO  1680 

p£;  *****  |>  p p ***** 

3:0  K2*K1+7:IF  K?r»120THEfl  40G: PRINT  "NO  MORE  ROOM  FOR  RPP'S":K2»K?-7: 

GOTO  1683 

403~ INPUT  ”XMIN,XMAX,YI1IN,Y!lAX,ZniN,Z!iAXa",ni  ,R2,R3,R4,R5,R6 
413  IF  P.l ]*P2TI!EI!  420: IF  P3]«R4THEN  430: IF  R53«RfT»'C::  440:00T9  45^ 

420  PRINT  "***  ERROR  M!IH  > TM  ***":GOTO  400 

433  PRINT  "***  E.RROR  YMIfl  ]■  YMAX  ***":GOTC  400 

44«  PRINT  "***  ERROR  ZMI'l  J»  ZiiAX  ***":GOTO  400 
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Table  A-XVI.  PLOTSOL  Prograa  Listing  (Continued) 


450  INPUT  "INPUT  OK  YES-1  N0-0",P 

4G0  IF  P-1 THEN  470: GOTO  400 

470  PRINT  *****  THINKIIIG... THINKING  ***** 

400  FOR  K-1T0  4:F(K)=R1  :J-K+4:F(»1)-R2:NEXT  K 

490  FOR  K-1T0  2:G(K)=R3:J»K+4:G(J)-R3:J»K+2:G(K+2)»R4:J-K+6:G(J)-R4:I!E 
XT  K 

500  M(l)=R5:ll(2)=R6:M(3)»RC:rK4)»R5:M(5)-R5:M(6)-RG:M(7)»RC:M(8)-R5 

510  FOR  K-Kl TO  K2:J-K-K3 

520  A(K)-F(J):H(K)=G(J):B(K)-;i(J):NEXT  K 

530  GOTO  310 

540  REM  ***  COX*** 

•550  K2-K1+7:IF  K2[-120THEN  5G0: PRINT  "NO  MORE  ROOM  FOR  B0X'S":K2-K2-7: 
GOTO  1G30 

5G0  INPUT  "X.Y.Z  OF  VERTEX-", F(1),G(1),M(1) 

570  INPUT  "HEIGHT  VECTOR-", R1 ,R2,R3:G0SUC  'n5(Rl ,R2,R3) 

530  IF  C]OTHFN  570 

590  INPUT  "WIDTH  VECTOR-", R4,R5,RG:G0SUD  ,05(R4,R5,R6) 

COO  IF  CjOTHEN  590 

G10  INPUT  "DEPTH  VECT0R«",R7,R0,R9:G0SUB  ,05(P.7,R3,R9) 

C20  IF  C]OTIIEH  CIO 

630  INPUT  "INPUT  OK  YES-1  N0-0",P 

640  IF  P-1  THEN  650: GOTO  560 

G50  PRINT  *****  THINKING... THINKING  ***" 

CGO  F(2)-F(l )+Rl :G(2)-G(1)+R2:M(2)-M(1)+R3 
670  F(3)-F  2}+R4:G(3)-G(2)+R5:f1(3)-fl  2}+R6 
600  F(4)-F(l)+R4:G(4)-G(l)+RD:ri(4)»M(l)+R6 
690  FOR  K-5T0  3:J-K-4 

700  F(K)-F(J)+R7:G(K)-G(J)+R0:M(K)-M(J)+R9:IIEXT  K 
710  GOTO  510 

720  REM  ***  ARC7  *** 

730  K2»K1+6:IF  K2[-120THEN  2G0:K2»K2-G:PRINT  "NO  MORE  ROOM  FOR  ARB7'S" 
:GOTO  1600 

740  REM  ***  ARR6  *** 

750  K2»K1+5:IF  K2[-120THEN  2G0:K2-K2-5:PRINT  "NO  MORE  ROOM  FOR  ARBu'S" 
.•GOTO  1630 

760  REM  ***  AR05  *** 

770  K2-K1+4:IF  K2[-120THni  2G0-.K2-K2-4: PRINT  "NO  MORE  ROOM  FOR  ARDS'S" 
:G0T0  1680 

760  REM  ***  ARC4  *** 

790  K2»i:i+3:IF  K2[»120THEfl  2G0:K2»K2-3:PRINT  "NO  MORE  ROOf!  FOR  AREA'S" 
- :GOTO  1650 

* > 800  REM  ***  RAU  *** 

* . 310  K2-K1+5.-IF  K2[-1?07HEN  020 :K2-K2-5: PRINT  "NO  MORE  ROOM  FOR  RAW'S": 

* :,oto  lo  in 

320  INPUT  "X,Y,Z  OF  VERTEX-" ,F( 3) ,G(3) ,M(3) 
i 830  IflPtlT  "HEIGHT  VECTOR- ", R1  ,R2,R3:GOSUB  '05(R1  ,R2,R3) 

5 340  IF  C]0T!O  030 
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Table  A- XVI.  PLOTSOL  Program  Listing  (Continued) 


850  INPUT  "WIDTH  VECTOR* ",R4,R5,R6:G0SUB  '05(R4,R5,RG) 

3G0  IF  C]0THEN  350 

370  INPUT  "DEPTH  VECTOR* ",R7,R3,R9:G0SUB  '05(R7,R3,R9) 

880  IF  C]0THEN  870 

390  INPUT  "INPUT  OK  YES=1  N0=0",P 

900  IF  P=lTHEfl  910 .‘GOTO  820 

910  PRINT  "***  THINKING... THINKING  ***" 

920  F(4)=F(3)+R1:G(4  *G(3)+R2:M(4)*M(3)+R3 
930  F G)=F  3)+R4:G(G  =G(3)+R5:M(G)-M  3)+RG 
940  F(l)»F(4)+R7:n(l)=G(4)+R8:M  1)=M(4)+R9 
950  F(2)=F(3)+R7:G  2 =G(3  +R8:M  2 «M(3)+R9 
9G0  F(5)=F(G)+R7:G(5)<i(6)+R«8:M(5)-M(G)+R9 
970  FOR  K*K1T0  K2:J*K-K3 
9G0  A(K)=F(J):H(K)=G(J):B(K)=N(J):NEXT  K 
990  GOTO  310 

1000  RBI  ***  RCC  *** 

1010  L9=L9+1 :IF  L9[*5TNEN  1020:L9»L9-1  :PRIIIT  "flllll  RCC'S  AND  TRC'S  ] 5" 
:G0T0  1C30 


1020  INPUT  "X,Y,Z  OF  VERTEX**" .V3.V4.V5 

1030  INPUT  "HEIGHT  VECTOR*", 111 .H2 .1 13: GOSUB  '05011,112,113) 

1040  IF  C]OTHEN  1030 

1050  INPUT  "RADIUS*", H4:fl-15:G0SUB  '0G(H4) 

10G0  IF  C]OTHEII  1050 

1070  INPUT  "INPUT  OK  YES-1  N0*0",P 

10.50  IF  P*1THLII  1090: GOTO  1020 

1090  PRINT  "***  THINKING... Til  INKING  ***" 

1100  PRINT  "RCC'S  AND  TRC'S  TAKE  A WHILE" 

1 1 1 0 L5= ( (L9-1 )*720/N) : LG*L5+( 360/11  )+l :H5«S0R(fll 12+H2 !2+H3 ! 2) 

1120  L7*SQR(H1!2+II.!!2) 

1130  FOR  K*1T0  3G0/N:N(K+L5)»0:J»K+L5:0(J)-H4*C0S(U*(K-1)):P(J)=II4*SIN 

D) 

1140  S1-N(J):S2*0M):S3=P(J) 

1150  IF  L7[]0THEN  1160:N(J)*-SGN(H3)*S3+V3:0(d)»S2+V4:P(J)*V5:GOTO  1190 


1160  N| 
1170  01 
1180  PI 

1 1 93  l)| 
1200  P.| 


«S1*H1/H5-S2*H2/L7-S3*H1*H3/(L7*H5)+V3 
■SI *H2/H5+S?*H1 /L7- S3*H2*H3/ ( L7*H5 ) +V4 
■S1*H3/H5+S3*L7/H5+V5 


0(J)*C0S(A 
P(J)*COS(E 


1M 


J)*SIN(A)) 

J)*COS(A)*SIN(E))-(0(J)*SIN(A)*Si::(F)) 


1210  IF  T-10THFN  1230 
1220  02*J+(  160/5!)  :N(02)*N(J)+H1 :0(02)*0(0)+H2:P(02)*P(J)+II3 
1230  IF  0!"J1TI!CN  1240:IF  X3[?X9THEN  1240:Xn-Q(J):X9-X3:73-R(J):Z9-Z8:G 
OTO  1250 


1240  0US1  IK  '07(0(J),R(J)) 

125)  : :L XT  K:  IF  T-lOTliCfl  1430 
1 2u  J J1  *I.6:02*L9*720/H:  FOR  K-OITO  02 
1270  l;(K)*(0(K)*C0S(A))+(ri(K)*SIN(A)) 

121:1;  r:(K)={P(:0*C0S(E))+(N(K)*C0S(A)*SIN(E))-(0(K)*SIN(A)*SIN{F)) 


Table  A-XVI.  PLOTSOL  Program  Listing  (Continued) 


1290  GOSUD  '07(Q(K),R(K)) 

1300  IIEXT  KsGOTO  1G80 
1310  REM  ***  TRC  *** 

1320  L9-L9+1 : IF  L9[«5TI!Efl  1330:LP»L9-1  :PRINT  "NIJM  RCC'S  AND  TRC’S  ] 5" 
:GOTO  1G80 

1330  INPUT  "X,Y,Z  OF  VERTEX-", V3,V4,V5 

1340  INPUT  "HEIGHT  VECTOR-", 111  ,H2,H3:G0SUB  '05(111,112,113) 

1 350  IF  ClOTHEil  1340 

1360  INPUT  "RADIUS  OF  CASE»",H4:G0SUB  '0 6(l!4) 

1370  IF  CjOTHEfl  1360 

1300  INPUT  "RADIUS  OF  T0P«",HG:N«15:G0SUB  '06(116) 

1300  IF  C]OTIIEM  1330 
1400  inn.'T  " INPUT  OK  YES-1  N0«O",P 
1410  IF  P-1  TIIEIi  1420: GOTO  1330 
1420  GOTO  1090 

1430  IF  A9-1TKFN  1440:L5«L6-1 : LG-L5+25 : 1 14-: IG : A9-1 :fi070  1130 
1440  A9«0:T4»II1  :T5-l!2:Tu-ll3 

1450  l!1“(T4*COS(A)yfCOS(E))-(TD*SIi:(A)*COS(E))-(TG*SIN(E)) 

14GC  U2-  T5*C0S(a))+(T4*SIN(A)) 

1470  H3-(T6*C0S(E) j+(T4*C0S(A)*SIN(E))-(T5*SlN(A)*SIfl(E)) 

1430  J1-L5+1 :J2-LG-1 :FOR  K-J1T0  02 

1490  N(K)-N(K)+T4:O(K)-O{K)+T5:P(K)-P(K)+T6:Q(K)»0(K)+H2:R(K)-R(K)+U3 
1500  GOSUD  '07(Q(K),R(K)) 

1510  NEXT  K:GOTO  1G00 
1520  REM  ***  SPH  *** 

1530  5?»S9+1:IF  S9[«5TilEfl  1GG0:PRIflT  "ERROR-HUM  SPH’S  ] GH 

1540  SO- 59-1: GOTO  1G00 

1GF.N  INPUT  "X,Y,Z  CENTER-", X(S9) ,Y(S9) ,Z(S9) 

1 5G0  INPUT  "RADIUS-", SO:GOSUB  'Of. (SO) 

1570  IF  CjOTliEII  1560 

1580  INPUT  "INPUT  OK  YES-1  N0»0H,P 

1590  IF  P-1  THEN  1G00;G0T0  1550 

1500  PRINT  "***  THINKING.. .THINKING  ***" 

1G10  r.G«(S9-l  )*24:N-15 

1G20  S3«(Y(S9)*C0S(A))+(X(S9)*SIN{A)):S4-(7.(S9)*C0S(E))+(X(S9)*C0S(f )* 
5i;:(E))-(Y(S9)*SIN(A)*5IN(E)) 

1b 30  FOR  K-1T0  360/ri j J«K+SG:!l(J)»Sn*C0S(N*(K-1 ) ):V(J)»S8*SIH(il*(K-l ) ) 
1G4C  S(J)-U(J)+S3:T(J)»V(J)+S4 

1650  IF  JHITHCII  1GC0: IF  XG[]X9THEK  1 GGO : XO-S ( J > : X9-S { J ) : 70-T ( J ) : Z9-T{ 
J):G0T0  1670 

1G60  GOSUD  '07(5(J),T(J)) 

1G7G  NEXT  I'. 

If  GO  IIEXT  I 
1C93  i:r»K9-L9-59 
1703  GOSUD  '04 

1713  PRINT  "***^*  START  PLOT  *****" 

1720  IF  I.9-0TUE:I  20G0 
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Table  A-XVI.  PLOTSOL  Program  Listing  (Continued) 


1730  PRINT  "WHEN  READY  TO  PLOT,  KEY  CONTINUE" ’.PRINT  HEX(07):ST0P 

1740  GO SUR  'lG(X8,X9,Za,Z9):G0SUB  '01 (A.E):G0SUB  '02:G0SUB  '03 

1750  FOR  KalT0  L9:G0SUB  '18(XG,X9,Z8,Z9):MW(K-1)*720/N)+1 

1760  M2»m+(360/N)-l:G0SUB  '11(Q(M1),R(M1)1:F0R  J-M1T0  M2 

1770  GOSUB  '12(Q(J),R(J)):NEXT  J:G0SUB  '12(Q(M1  ),R(M1)):H3-M2+1 :M4-M2+ 

(360/M) 

1780  GOSUB  '11(Q(M3),R(M3)):F0R  J»I13T0  M4:G0S1JB  '12(Q(J),R(J)):NEXT  J 
1790  GOSUB  'lZCQCfOJ.RtnOjjrGOSUr.  '12(Q(M1),R(M1)) 

1800  J 1 a!  11 + ( 1 80/ II ) : 0 2-M3+  ( 1 80/ N) : GOSUB  a 11  (Q(J1  ).R(J1 )) 

1010  GOSUB  1 12(Q(J2),R(J2)) 

1020  Jl=Ml+(90/N):J2-M3+(90/N) 

1030  GOSUB  *11 (Q(J1 ),R(J1)):G0SUB  ' 12(Q(J2) ,R(J2) ) 

1040  J1=J1+(180/N):J2®J2+(180/N) 

1050  GOSUB  'll (Q(J1 ),R(J1 )):G0SUB  ' 1 2 (Q( J2) ,R( J2) ) 

10C0  PLOT  [,,R]:11EXT  K 

1070  IF  S9]0THF.N  2000: IF  K9]0TIIEN  2300 

1800  INPUT  "SAME  SOLIDS  DIFF  VIEW  1=YFS  0»N0",J 

1090  IF  J]0THEN  1940 

1900  INPUT  "DO  YOU  WISH  TO  ADD  MORE  SOLIDS  1-YFS  0«N0",J 
1910  IF  JHTHEN  10 

1920  INPUT  "HOW  MANY  SOLIDS  DO  YOU  WISH  TO  ADD",J 
1930  E9-G9+1 : G9-G9+J : K9-G9  :GOTO  140 
1840  INPUT  "HEW  AZ,  EL  »",A,E 

1950  PRINT  "******  I'M  WORKING  HARD  - YOU  JUST  RELAX  ******" 

I960  X8-0 :X9«0: Z3-0 : Z9-0 

1970  f]laL9*720/N:F0R  K*1T0  111 

1580  8(K)-(0(K)*C0S(A))+(H(K)*SIfl(A)) 

1980  R(K)a(p(K)*C0S(C))+(H(K)*C0S(A)*SIN(F.))-(0(K)*SIN(A)*SIN(E)) 

2000  IF  K[]17HEN  2010:IF  X3[)X9THEN  2010:X8-Q(K):X9-X3:ZR-R(K):Z9-Z8:G 
OTO  2020 

2010  GOSUB  '07(0(10, R(K)) 

2020  NEXT  K 

2030  PRINT  "WHEW  READY  TO  PLOT,  KEY  CONTINUE" ’.PRINT  IIFX(07) ’.STOP 
2040  GOSUB  '10(X8,:'I  ,Z8,Z9):G0SUB  '01  (A,E):GOSUP  '02 ’.GOSUB  '03 
2050  i, OTO  1750 

2060  IF  S'-OTNFN  2200:PRINT  "WHEN  READY  TO  PLOT,  KEY  CONT!N!'E":PRINT  H 
LX (07): STOP 

2070  GOSUB  1 1 8 ( X0 , X9 , Z8 , Z9 ) : GOSUB  '01 (A,E):GOSUR  '02:G0SUB  '03 
2080  FOR  I>1T0  SO: GOSUB  'lG(Xa,X9,Z8,Z9):Sl-(K-l )*24+l :S2-Sl+23 
2090  KjSlir.  '11(S(S1),T(S1)):F0R  .1-51  TO  52 
180  f'C5!)!'  'lO’sfJJ.TtJJJ.’NEXT  J:GOSUB  '12(S(S1),T(S1)) 

2110  PLOT  r,,;il:.;EXT  K 
2120  IF  l^iOTNEil  2300 

2130  INPUT  "SAME  SOLIDS  DIFF  VIEW  1CYFS  0«f!0",J 
214  • IF  JlUTHEN  2190 

2150  i.lFIIT  HDU  YOU  !,'IS1!  TO  ADD  MORE  SOLIDS  1-YFS  0»N9" , J 

21  :o  if  .i[-07sri  lo 
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Table  A-XVI . PLQTSOL  Program  Listing  (Continued) 


INPUT  "HOW  MANY  SOLIDS  DO  YOU  WISH  TO  ADD\J 
E9=G9+1:G9aG9+J:K9aG9:G0T0  140 
INPUT  "HEW  AZ,  EL  a",A,E 

PRINT  "******  I'M  WORKING  HARD  - YOU  JUST  RELAX  ******* 

X8a0 : X9=0 : Z8a0: Z9a0 

FOR  KalT0  S9:S3a(Y(K)*C0S(A))+(X(K)*SIN(A)) 

S4a(Z(K)*C0S(E))+(X(K)*C0S(A)*SIN(E))-(Y(K)*SIN(A)*SIN(E)) 

SGa(K-l )*24+l :S7°SG+23 

FOR  J=S6T0  S7:S(J)aU(J)+S3:T(J)aV(J)+S4 

GOSUB  ' 07 (S( J) »T( J) ) 

NEXT  J:NEXT  K:IF  L9]0T1IEN  1 970: GOTO  2%0 
PRINT  "WHEN  READY  TO  PLOT,  KEY  CONTINUE": PRINT  HEX(07):ST0P 
GOSUB  ' 1 8 ( XO , X9 , ZG ,Z9 ) : GOSUB  *01 (A.E) rGOSUB  '02:G0SUB  '03 
FOR  1=1 TO  K9 

J=KU(I+1)-K0(I):IF  J-8THEN  2340: IF  Ja7THEN  2810 
IF  J-STHEII  2520:  IF  Ja5THE!l  2720:  IF  J-4THEN  2800 
PRINT  "ERROR--CANNOT  PL0T":ST0P  .-GOTO  10 
GOSUB  '10(X0,X9,ZB,Z9) 

M1=K3(I) 

GOSUB  ,ll(C(Hl),n(Ml)) 

FUR  J=M1T0  Ml +3: GOSUB  '12(C(J),D(J)):NEXT  J 
GOSUB  '11 (C(Ml+4),D(Ml+4)) 

FOR  J*M1+4T0  M1+7:G0SUB  '12(C(J),D(J)):NEXT  J 
GOSUB  *11(C(M1  ,D  Ml  ) 

GOSUB  '12(C  M1),I)(M1)) 

GOSUB  '12(0(111+3} ,D(Ml+3)) 

GOSUB  ,12(c(m1+7},D(M1+7)) 

GOSUB  '12  C(Ml+4j,D(ni+4)) 

GOSUB  '12  C(M1),D(M1)) 


I : GOSUB  '12(C(Ml+2),D(Ml+2)) 


GOSUB  '12  C(Ml),D(f11)) 

GOSUB  'll  (uMl+l ) tD(Ml+l  )):GOSUB  '12(C(m+l),D(m+l)) 

gosub  •i2(cnn-K>).ncm+D}} 

GOSUB  'll(o(Ml+t'.),D(Ml+0)) 

GOSUB  '12(c(:il+<;),D(Ml+G)}:G0SUB  '12(C(M1+2),D(M1« 
PLOT  [,,RJ 
GOTO  2BG0 

:il«KB( I ) : GOSUB  • 18(XC,X9,ZC,Z9) 

GOSUB  'll (C(Ml).D(m)) 

TOR  J=illTO  Ml+3:G0SUn  '12(C(J),D(J)):NEXT  j 
GOSUB  '12(C(i11),D(Ml)):G0SUB  h2(c(Ml+4),D(Ml+4)) 
GOSUB  '12(c(m+5),D(m+5)}:G0SUB  '12(C(M1+3),D(MH 
gosub  'ii!c(m+2),n(m+2)):Gosun  '12(c(;ii+5),d(mih 
gosub  'ii (c(in+i ).n(rn+i)):Gosun  'iztccm+AhDtnn 

PLOT  [„R] 

GOTO  23G0 

: 11  =. ( I ) : GOSUD  ' 1 8 ( XO , X9 ,Z3,Z9 ) 

gosub  'll (c(rii).D(m)) 

FOR  J=M1T0  111 +3: GOSUB  '12(C(J),D(J)):NEXT  J 


GOSUB  '12(C(Ml+5),D(m+5)):G0SUB  ' 12(C(Ml+3> ,D<m+3) ) 
GOSUB  'll (C(Ml+2), 0(111+2  :GOSUB  '12(c(;il+5  ,d(M1+5)) 
gosub  'll (c(mi+i ),d(mi+i)):gosub  •i2(c(rn+^) ,n(rn+4) ) 


X 
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Table  A-XVI.  PLOTSOL  Program  Listing  (Continued) 


2C.40  gosub  ,n(c(m+/i},i)(:n+4)) 

2G50  FOR  J-M1+4T0  M1+6;G0SUB  ' 1 2 (C( J ) ,P( J) ) rflEXT  J 
2GG0  GOSUD  '12(C(m+4),D(Ml+4)):G0SllD  '12(C(M1),D(?I1)) 

2G70  G0SUD  ,12(C(fil+3j,D(lil+3)):nosilB  '12(0(01+4), D 111+4)1 
2r.no  GOSUD  '11(C  l11+G),0(!11+u)):G0SUn  ,12(C(m+2),D(tn+2)j 
2090  GOSUB  '11  (C(M1+1  ),D(ril+l )):GOSUH  '12(C(l!l+5)  ,D(Ml+f>}) 

2700  PLOT  [„R] 

2710  GOTO  2860 

2720  m-Ks(i):nosun  'lofx.wa.zo) 

2730  10SUD  '11 (C(ni),D(Ml)) 

2740  FOR  J-M1T0  M1+4:G0SIJB  '12(C(J),D(J)):flD:T  .1 
27.00  GOSUC  '12(C(M1),D(M1)):G0SUB  ' 1 2 (C(M1  +3) ,rj(*!l  +3 ) ) 

2700  GOSUC  'll  C(M1+4),D  l!l+4)):G0SUB  '12(C(M1+1),D(M1+1)) 

2770  GOSUD  'll  (C(II1+2),D(M1+2)):G0SIIB  'lJ»(C(rn+4) »D(m+4)) 

27SG  PLOT  [„R] 

2790  GOTO  2TG3 

2300  ?ll-i;r,(I):GOSIIB  1 13(X0,X9,Z3,Z9) 

2310  GOSUC  '11  (C(lll).D(m)) 

2320  FOR  J**l  11T0  M1+3:G0SUB  '12(C(J),D(J)):IIEXT  J 
2330  GOSUD  '12(C(I11),D(M1)):G0SUR  '12(C(Ml+2)t D(Hl+2)) 

2340  GOSUB  ' 1 1 (C(Ml+3) ,D(Ill+3) ) rGOSUB  '12(C(M1+1),D(M1+1)) 

2330  PLOT  [,,R] 

2300  NEXT  I 

2870  INPUT  "SAME  SOLIDS  DIFF  VIEW  1-YES  0«N0\J 
2330  IF  J]OTHEN  2930 

2390  INPUT  "DO  YOU  WISH  TO  ADD  MORE  SOLIDS  1-YES  >NO\J 
2900  IF  J[— OTlIEri  10 

2910  INPUT  "NOW  MANY  SOLIDS  DO  YOU  WISH  TO  ADD\J 
2920  E9-GC+1 : G0-G9+J : K9-G9 :GOTO  140 
2930  INPUT  "NEW  A7„  EL  «",A,E 

2940  PRINT  "******  I'M  WORKING  HARD  - YOU  JUST  RELAX  ******" 

2930  X3-0 : X9-0 : ZB-0 : Z2-0 

2930  K2-K3(K.1+1  }-l : FOR  K-1T0  K2 

2970  C(K)-(N(K)*COS(A))+(A(K)*SIN(A)) 

2930  l)(K)-(B(K)*COS(E))+(A(K)*COS(A)*Siri(r))-(H(K)*SIN(A)*SIN(F)) 
2990  IF  KfjlTUEN  3000:X3-C(1):X9-X3:Z3-D(1):29-Z3:G0T0  3010 
3000  GOSUD  '07(C(K),n(K)) 

3013  NEXT  K 

3020  IF  S9]0THEN  2220:  IF  L9]l)THEI!  1970 

3033  GOTO  2230 

3340  OEFFIi' 13(1)1  ,II2,V1,V2) 

3U?:o  FI -900/ (UN-Ill} 

3030  F2-900//V3-V1) 

3070  IF  Fl’INTWCN  3090 
3r;>0  F1-F2:G0T0  3100 
3‘jOO  F2aFl 

31 UO  SELECT  PLOT  414 
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Table  A- XVI.  PLOTSOL  Program  Listing  (Continued) 


3110  PLOT  [„R],[100f100,U] 
3120  GOSUD  'IO(UI.VI) 

3130  RETURN 
3140  nEFFN'10(X,Y) 

31  59  XlcINT(Fl*X-X0+,5) 

3160  Yl*INT(F2*Y-Y0+.5) 


3170  X0-X0+X1 
3180  Y0-Y0+Y1 
3190  RETURN 
3200  DEFFIHKX.Y) 
3210  G0S1JB  ‘lO(X.Y) 
3220  PLOT  [Xl.Yl.lJ] 
3230  RETURN 
3240  DEFFN,12(X,Y) 
3200  OOSIJC  'IO(X.Y) 
3260  PLOT  TXl.Yl.n] 
3270  RETURN 
3280  DEFFN*01(A,E) 


3290  SELEC 


PRINT  414 


3300  PLOT  [„R] 

3310  PLOT  [2,,C],[23,,S],[320,20,U3 
3320  PRINT jSINfi  3330.A.E; 

3330  5,A«- 

3340  SELECT  PRINT  305(64) 

3350  RETURN 


3360  OEFFN'02 
3370  PLOT  [,,R] 

3380  SELECT  PRINT  414 

3390  PLOT  n00,70,U], [0,5,0], [150,0,Dl,[0,-5,D] 
3400  PLOT  [1»»C].[14,,S],[,,RJ,[95,50,UJ 
3410  PRINTIlSIln  3420,0; 

3420 

3430  PLOT  Ll..C],n4„S],[n5,0,U3 
3440  PRIliTUSIllO  3450,150/Fl; 

3450  iijW.? 

3460  PLOT  :„R] 

3470  SELECT  PRINT  005(64) 

34*. . RETliRi! 


3490  DEFFN'03 
85  i ) PLOT  " ,R] 

3518  PLOT  C\905,0],r„R], [999,0, D] 

3523  RETURN 
3581.  :tFFi:'04 
8540  r*-io 

3550  INPUT  "AXES  LENGTHS  ADJUSTED  1-YES  <>[!0",R 

3560  IF  n]-lTi:EN  3G20.-F-F+10 
357  ' GOSUR  ,lU(X.r;,X9,ZL,Z9) 


92 


Table  A-XVI.  PLOTSOL  Program  Listing  (Continued) 


3530  PLOT  r„R),rF,F#U],[50O,0tD]t[#tR] 

3590  PLOT  [FtFtU]t[0,500,D]i[ttR] 

3600  PRINT  "IF  THESE  LINES  ARE  OF  EQUAL  LENGTH  THEN  SCALE  IS  OKAY" 

3610  PRINT  “IF  NOT  THEN  READJUST  WITH  THE  SCALE  ADJUST  BUTTON": GOTO  35 
50 


3620  RETURN 

3630  DEFFN'05(G,B,D) 

3640  C«0: IF  GI2+B!2+D!2]0THEN  3660 

3650  C-1  jPRItlT  "***  ZERO  VECTOR  TRY  AGAIN  ***":PRINT  IIEX(07) 
3660  RETURN 
3G70  DEFFN'Oo(R) 

3600  C«0:IF  RjOTUEN  3700 

3690  C-1 SPRINT  "***  RADIUS  [-  0 TRY  AGAIN  ***":PRIHT  HEX(07) 

3700  RETURN 

371o  DEFFN'07(C,D) 

3720  IF  C]XPTHEN  3740: IF  C[X3THEN  3730:G0T0  3750 
3730  XO«C:GOTO  3750 
3740  XO-C 

3750  IF  D[Z3T1IEN  3760: IF  D]Z9THEN  3770:G0T0  3730 
3760  Zr>D:G0T0  3700 
3770  Z?W) 

3730  RETURil 
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Table  A-XVII.  D ARB IN  Program  Listing 


10  DIM  X{8).Y(0),7.(0),F(4,C),T(6),A(G),R(6),C(f»),R(3,3),Fl(G), 
X1(3),Y1(8),Z1(0) 

20  &]-## 

30  SFACE  MM 

40 

50  %F0UR  POINTS  NOT  CO-PLANAR  IN  FACE  if  if  if  if  ON MM 

oo  %wm 


70  C9-. 00005 

30  INPUT  MtlO.  OF  VERTICES", V: INPUT  "MO.  OF  FACES", H 
PO  ! 1-0 : f i»0 :0-0 : SELECT  PRINT  000(04) 


100  FOR  J-1T0  V: PRINTUSING  20,J 

110  INPUT  "X(J)",X(J):INPUT  "Y(J)",Y(J):INPiJT  "Z(U)",Z(J) 

120  M-!1+X(J):N«II+Y(J):0«0+Z(J):X1(J)-X(J):Y1(J)-Y(0):Z1(J)-Z(0) 
130  IF  ::(J)ho  THEM  140:X<J)-.00001 
1 40  IF  Y(J)  “o  THEN  150:Y  J)«.00001 
130  IF  Zijjtlo  THEN  lf.0:Z(J)«.00rX)l 

100  NEXT  J.-SELECT  PRINT  215(80) {PRINT  {PRINT  " INPUT" {GOSUD  470 


170  FOR  J-1T0  1 1 { C— 0 


130  INPUT  "FACE  !iUMnCR5M,G{IF  G-0  THEN  210{F1(J)-G 

190  FOR  K-1T0  4{P-10!(4-lO 

200  F(K,J)«i::T(n/P){G-r-F(K,J)*P{HEXT  It 

210  PRIilTUSI.IG  CO, FI  (J)  ;{NEXT  J.-PRIHT  iE-0 

220  FOR  J-1T0  1 1 { FOR  It-ITO  3*.L=F(K.J) 

2?0  R(K»l)BX(L){R(lt,2)-Y(L){R(K,3)-Z(L) 

240  ixt  i; {Gosiid  r>io 


250  A(J)»A1{D(U)-D1{C(J)-C1 
2C0  J1-F(4,J){IF  Jl-0  THEN  300 

270  n>(A(J)*::(Jl)+ri(J)*Y(Jl)+C(J)*Z(Jl)-l)/S0R(A(J)!2+P.(J)!2 
+C(o)!2) 

230  IF  ADS(n2)[.u01  THEN  300{E-1 

290  PRINTUSIIIG  50,F(1 ,U),F(2,J),F(3,U),F(4,J),D2 

3C0  :JEXT  J{IF  E-0  THEN  310 .‘GOSUD  390.-PRINT  "SOLUTION" .-GOSUD  470 

310  INPUT  "INSIOE  ARO  (YES-1 ,NO«n)",Ji IF  J-0  THEN  00{ 

PRINT  "THICKNESS  ARE" 

32)  1-’  1/ V {ll-M/V { O-O/V { FOR  J-1T0  IUSELECT  PRINT  005(04) 

330  PRINTUSING  30,F1(J){ INPUT  "THICKIIESS",T(J) 
L«i-A(o)*n-D(j)*n-c(v-'*0{Ti«T(J) 

350  IF  c'o  THE!!  300{T1— Ti 

3CC  r.»l/(l-Tl*SOR(A(,l)!2+0(J)!2+C(U)!2)) 

370  A(J)-A(J)*E{B(J)»D(0)*E{C(J)-C(U)*E 
3N0  SELECT  PRINT  :>15(nO){PRINT!)SING  40,T(J);{HEXT  J.-PRIHT  { 
GOSUD  3P0.-PRINT  "INSIDE  ARU"{GOSUR  470{GOTO  00 


393  FOR  J-TTO  V:P-0 


400  FOR  K-1T0  lUFOR  L-1T0  4 
410  IF  F(L,K)-J  THEN  420{NEXT  L{GOTO  440 
420  P=P+l{R(P,l)-A(lt){R(P,2)-D(K){R(P,3)-C(lt) 
430  IF  P-3  THEN  450 
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Table  A-XVII.  DARBIN  Program  Listing  (Continued) 


440  MEXT  1C:  ROTO  460 

450  ROSUD  510:X1(J)«A1+SGN(Al)*C9:Yl(J)*Rl+SRH(m)*C9 
:Z1(J)»C1+SGN(C1)*C9 
460  NEXT  J: RETURN 

470  FOR  J-1T0  V STEP  2: IF  J«V  THEN  490:K«J+1 

480  PRINTUSINR  40,X1(J);YHO);Z1(J);X1(IC)jY1(K);Z1(K):ROTO  600 

490  PRINTUSING  40,Xl{J);Yl(J);:i(0) 

500  NEXT  J: RETURN 

510  D1“R(1 ,1)*R(2,2)*R(3,3)+R(3,1 )*R(1 ,2)*R(2,3)+R(2,1)*R(3,2)* 
R(1 »3/-R(3,l J*R(2,2)*R(1,3)-R(2,1 }*R(1 ,2)*R(3,3)-R(1 ,1 )* 
R(3,2)*R(2,3) 

520  IF  D1[]0  THEN  530:01-.00001 

530  A1-(R(1,2  *(R(2,3)-R(3,3))+R(2,2)*(R(3,3)-R(1,3))+R(3,2)* 

R 1,3  -R(2,3)  )/Dl 

540  B1»(R{1,3  *(R(2,1  -R(3,1))+R{2.3)*(R(3,1)-R(1,1))+R(3,3)* 

R 1,1  -R(2,l)  )/Dl 

550  C1*(R  1,1  *(R(2,2  -R(3,2})+R(2,1)*(R(3,2)-R(1,2))+R(3,D* 
(R(1,2)-R(2,2)))/01 
560  RETURN 
570  END 


i 
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Table  A-XVIII.  BOXIN  Prograa  Listing 


10  DIii 

2a  zmm  and  mm  vectors  are  hot  normal,  ahole  —mm 

40  SELECt'pRINT  215(64} :SELECT  D :B$(1)»"LENGTH": 

B$(2)«"WIDTH  " :B£»(3)***IIEinilTH 
50  INPUT  "VERTICE",P(1),P(2),P(3) 

60  INPUT  "LENGTH  VECTOR" ,V(1 ,1 ) ,V(1 ,2),V(1 ,3) 

70  INPUT  "WIDTH  VECTOR", V(2, 1 ),V(2,2),V(2, 3) 

80  INPUT  "HEIGHT  VECTOR" ,V(3,1 ),V(3,2),V(3, 3) 

90  PRINT  : PR I NT  "IHPUT":G0SUB  300 

100  FOR  0*1  TO  3: FOR  K-1T0  3:U(K)*V(U,IC):IIEXT  K 

110  GOSIIB  3l0:M(0)-;il:F0R  K-1T0  3:W(0,K)*tl(IC):NEXT  K: NEXT  0 

120  E-0:K*l:L-2:n-3:G0SUn  330 


130 

140 

150 

ICO 

170 

100 

190 

200 

210 

220 

230 

240 

230 

260 

o 7 i 

<./  >1 

200 

200 

300 


310 

820 

330 


35o 

860 


K*3:L*l::i*2:G0SUB  330 
!>2:L*3:M*1  :G0$IIB  330 
IF  E-0  THEN  220 

FOR  J-1T0  3:A(1,0)«P(J):A(2,J)*P(J)+V(1,J): 
A(3,J)*A(2,U)+V(2,J):A(4,J)*P(J)+V(2,J) 

FOR  K-1T0  4:L-K+4:A(LiJ)-A(K,J)+V(3,J):NEXT  K:NEXT  J 
PRINT  "ARC  SOLUTION": FOR  J-1T0  G STEP  2:K«J+1 
PRINTUSING  30,A(J,1);A(J,2);A(J,3);A(K,1);A(K,2);A(K,3): 
NEXT  J 

FOR  K-1TO  3: FOR  J-1T0  3:V(K,J)-ll(K,J)*tl(K):NEXT  JrllEXT  K 
PRINT  "ROX  SOLUTION" :G0SUB  300 
INPUT  "INSIDE  COX  (YES-1 ,N0«0}"  I:IF  >0  THEN  50 
INPUT  "THICKNESSES",T(1 ) ,T(2) ,T(3) ,T(4) ,T(5) ,T(6) 

PRINT  "TIIICKNCSS  ARE": 

PRINTUSING  30,T(l);T(2)iT(3):T(4)iT(5);T(6) 

FOR  O-ITO  3:P1(J)*P(0)+W(1 ,J)*T(1)+W(2,J)*T(2)+W(3,J)*T(3) 
FOR  K-1T0  3:L»K+3:Vl(i;tJ)-V(K,J)-'.;(K,j)*(T(K)+T(L)) 

NEXT  K.’NEXT  0 
PRINT  "INSIDE  COX" 

PRINTi'SIfiG  30, PI  (1 ) jPl  (2)  ;P1  ( 3)  ;Vl  M ,1 ) • VI  (1 ,2)  ;V1  (1 ,3) , 
V1(2,1);V1(2,2):V1(2.3);V1(3,1 );Y1(3,2);V1(3,3):G0T0  50 
PRINTIISING  30,P  l);P(2);P(3)iV(l,l);V(l,2);V(l,3),V(2,l); 
V(2,2);V(2,3) iV(3,l ) ;V ( 3,2) ;V( 3,3): RETURN 
Ml*SOn(lj(l)!2+U(2)!2+U(3)!2) 
U(l)-U(l)/ni:l!(2)-ll(2)/ni:U(3)*U(3)/m:RETllRN 
D-i'(K,l)*U(L,l)+U(K,2)*U(L,2)+U(K,3)*W(L,3): 

IF  ACS (0)].  1)007  THEN  340 .’RETURN 
!!(l)-'/(L,2  *U(M,3)-U(L,3  *W(:!,2): 

U(2)-U  L ,1 5*VJ(fl,3)-VJ(L.3)*H(!l,l ): 

U(3)*t'(L,l Mm,2  -N(L.2)*H(M.l ) 

S-l : IF  U(1;*W(K,1 5+U(2)*W(K,2)+U(3)*VI(K,3)]0  THEN  360:S— 1 
FOR  0*1  TO  3:H(K,0)-S*tl(J):NEXT  J:D-ARCCOS(D): 
E*l:rRINTUSING  2C,r.S(K),D$(L},D:RETURH  :END 
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Table  A-XIX.  RAWIN  Program  Listing 


Vi  Oir:  P(r,),V(3,3),T(G)tU(3)#li{3,3),H(3)fA(C#3),P1(3)#Vl{3,3)f 


20  wtm  Ann  mm  vectors  are  mot  horjial,  aixjle 

40  SELECT*  PR  HIT  21 5 (34): SELECT  0 :P,S(1)»"!!EIG:;T": 
li$(2)="l!r;Tl!  ":f.$(0)a"PEPTII  " 

30  I.'IPUT  "VK;VnCF.",P(l),P(2},P(3) 

GO  INPUT  "HEIGHT  VECT0R",V(1 ,1 ) ,V(1 ,2) tV(l ,3) 

70  IilPtrr  "WIDTH  VECTOR", V(2,1),V(2,2),V(2, 3) 

80  ir-Pl!T  "DEPTH  VECTOR", V(3,1),V(3,2),V(3,3) 

30  PR  HIT  : PRINT  "INPUT" :G0SUR  300 

100  FOR  J-1T0  3: FOR  K*1T0  3:U(K)-V(J,K):HEXT  K 

110  GOSUB  310::i(J)*ni :F0R  K-1T0  3:l!(J,KH(K):f!EXT  KrllFXT  0 

123  E=0:K»1:L*2:If3:G0SU[;  330 

130  i;»3:L-l:il»2:GUSUU  330 

140  K»2:L»3:H«1:G0SUG  330 

150  IF  E*0  THEM  220 

ICO  FOR  J-11'0  3:A(1,J)»P(J):A(2,J)«P{J)+V(1,J): 

A( 3 , J ) *A(2 ,0 )+V{ 3, J) 

170  A(4,J)-P{0)+V(3,J):A{5,J}»P(il)+V(2,0):A(u,Jl-A(5,JHV(3fJ): 
NEXT  J 

130  PRINT  "ARB  SOLUTION": FOR  J*1T0  6 STEP  2:K»J+1 
100  PRIHTUSIHG  30,A(J,1 ) ;A(J,2) ;A(J,3) ;A( K,1 ) ;A(K,2) ;A(K,3) : 
NEXT  J 

200  FOR  K-1TG  3:F0R  J*1T0  3 : V ( K, J ) »l .*{ K , J ) *M { K) : N’EXT  J:NFXT  K 

210  PRINT  "RAW  SOLUTION"  :GOStlD  300 

28)  INPUT  "INSIDE  RAM  (YFS-1 ,N0«O)",I:IF  1-0  THEN  SO 

230  INPUT  "THICKNESSES" fT{l),T(2),T(3),T(4),T(r.) 

240  PRINT  "THICKNESS  ARE": 

PRINTUSING  30,T{1 ) ;T(2) ;T(3) ;T(4) ;T(5] 

250  R«(:i{l)/i:(2))!2:Tl»T 
T:T*SQP.(T1  +T1  *R)*$D+T 


u 

d 


!2:Sn-SGN(T(r.) ) :T3«T(3)+7(*) : 
1):TG»S0R(T1+T1/R)*S9+T{?) 
r,n  FOR  3«1T0  3:P1{J)»P(J)+U(1,J)*T(1)+W(2,J)*T(2)+I!(3,.l)n(3) 
17 V VI  ;i , J -V(l , J -i;(1  ,J)*TS:V1  (2,J  «V(2, J)-U(2, J)*T0: 
Vl(3f4)«V(3,J}-U{3fJ)n3:f!EXT  J 
800  PRINT  "INSIDE  RAH" 

800  PP.INTUSIfIG  30, PI  (1 ) ;P1  (2) ;P1  (3)  ;V1  (1 ,1 ) ;V1  (1 ,2)  ;V1  (1 ,3) , 

VI (2,1 ) ;V1  (2,2) ;V1  (2,3  ;V1  (3,1 ) }V1  3,2) ;V1  (3.3) -.GOTO  50 
300  PRIIITUSING  30,P(1);P(2);P(3);V(1,1);V(1,2);V(1,3),V(2,1); 

V(2,2);V(2.3);V  3.1);V(3.2)sV(3,3):Rrn)R;i 
I!  1 0 : !l  =nnn(.i(  1 ) \ 2+u( ?.)  !2+<J( 3)  !2) 

3.80  ,!(l)»t!(l)/:il:U(3)-U(2)/in:U(3)«ll(3)/m .‘RETURN 
330  ;)»:'(K,l}*U(L,l)+i:(K,2)*U(L,2)+lI(K,3)*H(L,3): 

IF  A;!5(D)3.00u7  THEM  340:RETllRri 


i*?J  l!(l}»H{L,3)*U(ti,3)-U(L,3)*U(M,2): 
"(2)=l  (L,l  *U(M,3  -1!(L,3  *H(I1,1  : 
U(3)<;(L>l}*ti(M,2)-L,(L.2)*U((lall 
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Table  A-XIX.  RABIN  Prograa  Listing  (Continued) 


350  S»1:IF  U(1)*U(K,1)+U(Z)*W(K,2)-HJ(3)MK,3)]0  THEN  360:S-1 
3C0  FOR  J-1T0  3:W(K,J)*S*ll(J):NEXT  J:D«ARCCOS(n): 
E-1:PRIIITUSING  ,?0,B$(K),B$(L),D:  RETURN  :EHO 
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Table  A-XX.  TRCIN  Prograa  Listing 

10  DIM  P{3),V(3),T{3),M(3),Q(3)»U(3) 

30  SUARNINC,',  TOP  RADIUS  OF  INSIDE  TRC  — #*##.*# 

40  SRADIUS  RESET  TO  .0001  AND  THICKNESS(2)  - 0 
50  SELECT  PRINT  215(30) 

60  INPUT  "VERTICE",P(1),P(2),P(3) 

70  INPUT  "HEIGHT  VECTOR" ,V(1 ) ,V(2) ,V(3) 

00  INPUT  "RADIUS  OF  RASE", R1 : INPUT  "RADIUS  OF  T0P",R2 
90  INPUT  "THICKNESSES" ,T(1),T(2),T(3):PRINT  SPRINT  "INPUT" 

100  PRINTUSIHO  20,P{1);P(2);P(3);V(l)iV(2);V(3),R1;R2 
110  PRINT  "THICKNESSES  ARE": PRINTUSING  20,T(1);T(2);T(3) 

120  ll«SQP.(V(1)!2+V(2)12+V(3)I2):G-H/(R1-R2) 

130  T2»T(2):F"T(3)/SIN(ARCTAfl(n)):Sl«Rl-F-T(l)/fi:S2»R2-F+T(2)/G 
140  IF  S2>.0001  THEN  150:PRIHTUSING  30,S2:PRINTUSING  40: 
T2»(F-R2)*C.:S2-.0001 

13.)  FOR  J-1T0  3:U-V( J)/ll:Q(J)-P(J)-HJ*T(l ) 

ICO  U(J)-V(0)-U*(T(1)+T2):HEXT  J 
170  PRINT  "INSIDE  TRC" 

mo  PRINTUSIHG  20,0(1 ) ;Q(2) ;Q(3) ;ll(l ) ;l.'(2) ;ll(3) ,S1  *,S2 
190  GOTO  00: CUD 
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Table  A-XXI.  RECIN  Program  Listing 


10  OIM  P(3),V{3.3),T{3),U(3),M{3,3),ft(3),Pl(1),Vl(3,3),B5(3) 

20  *##########  AND  ##########  VECTORS  NOT  NORMAL,  ANGLE-###.## 
30  S-####.#### 

40  SELECT  PRINT  21 5(04) {SELECT  0 :B$(1  ^"HEIGHT": 
B$(J?)"**SET1I-MAJ0n*'  tB$(3)*“SEf1I-MinOR'' 

50  INPUT  "VERTICE",P(1),P(2),P(3) 

CO  INPUT  "HEIGHT  VECT0R",V{1,1 },V(1,2),V(1,3) 

70  INPUT  "SEMI-MAOOR  VECTOR" ,V(2,1),V(2,2),V(2,3) 

00  INPUT  "SEMI-MINOR  VECT0R",V(3,1),V(3,2),V(3,3) 

90  PRINT  ‘.PRINT  "INPtlT":GOSIin  270 

100  FOR  ,1-1  TO  3:F0R  K-1T0  3:U(K)«V(0,K):NEXT  K 

110  GOSUB  2G0:U{J)-m:F0R  K-1T0  3:W(0,K)-U(K):NEXT  KsHEXT  J 

120  E-0:K-l:L-2:H-3:G0SI)n  300 

130  K-3:L-1:!1-2:G0SUB  300 

140  K-2:L»3:M«1:G0SHB  300 

150  IF  E-0  THFN  180 

ICO  FOR  K-1T0  3: FOR  J-1T0  3:V(K,J)«U(K,i1)*!(K):NE.?r  JsNFXT  K 

170  PRINT  "SOLUTION" :GOSUR  270 

130  INPIfT  "INSIDE  REC  (YES-1  ,N0-0)",I: IF  1-0  THEN  SO 

150  INPUT  "TU ICKNCSSES" ,T  ( 1 ) ,T  ( 2 ) ,T(  3 ) 

200  PRINT  "THICKNESS  ARE": 

PRINTUSING  30,T(1);T(2)jT(3) 

210  T3-T(3):IF  M(3)/n{2)}.8  THEN  220: 

T3-T3/( 1 ,016447*(M( 3)/M( 2) ) l .071834) 

220  FOR  J-170  3:P1  (J j-P(J)"H/(l  .J j*T(l ) 

230  V1(1,J  -V(liji-ll(l,,m(T(l)+T(2)) 


240 


VI  1,0  -V  1,0 J-MJU )*(T(  M 2)) 

Vl(2,J)-V(2,il)-!/(2,.l)n(3):Vl(3,J)-V(3,J)-U(3,,l)*T3:NEXT  J 
2jU  PRINT  "INSIDE  REC" 

2C3  PRIIJUSING  3s).ri(l)‘,Pl(2);Pl(3)iVl(l,l);Vl(l,2)iVl(l,3), 

VI (2,1 ) jVl (2.2);V1 (2.3  ;V1 (3,1 ) ;V1  (3,2) -,V1  (3.3):GOTO  50 
27U  PRINTUSING  30,P  1);P(2  ;P(3);y(l,l)‘,V(l,2)‘,V(1.2),V(2,l); 
V(2,2);V  2,3);V  3,1);V3.2);V(3,3):PXTHRN 
0 M1-S0R(U(1 ) !2+U(2)  !2+U(3)  12] 


* 

kU 


230  t?(n-H(l)/ill:U(2)-U(2)/?n:U(3)-H(3)/»ll:RETHRN 
300  0-H(K,l)*t/(L,l)+l'(K,2)*tl(L,2)*Mt(K,3)*tl(L,3): 
IF  ARS(D) ].00H7  THEN  310:RETURI! 

310  U(1 - 

U(? 


H L,2)*W 
H(L,1  )*U 


M,3)-ll(L,3)*H(H,2): 

M,3)-II(L,3  *V  M,l): 

, . ( ...M,2  -U  L.2)*U(M.1  j 

rn  S-1:IF  U(1)*W  K,1  ♦U(?)*H(K,2)+<J(3)*t,(K»3)3»  THEN  330:S-1 
133  FOR  ,1-1T0  3:W(K,J)-S*U(J):NEXT  J:0-ARCCOS(0): 
E-l:f*niflTUSirrt  20,B$(K),r‘$(l.),n:RETURN  :END 
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Table  A-XXII.  TECIN  Program  Listing 


1C  DIM  P(3),V(3.3  ,T(3).U(3),W(3,3),ri(3),Pl(3),Vl(3,3),fj(3), 

A1 (2) ,A2(2) ,D1 (2) ,B2(2) 

20  SSEMI-MAJOR  AND  SEMI-MINOR  VECTORS  NOT  NORMAL,  ANGLE  ■-#?.#? 

30 

40  SRATIO  OF  MAJOR  AXIS  -M.###  RATIO  OF  MINOR  AXIS  -##.•?**# 
50  «TOP  MAJOR  AXIS  -####.#*«  TOP  MINOR  AXIS  -MW.mf 
CO  SELECT  PRINT  21 5 (DO) .'SELECT  D 
70  INPUT  "VERTICE",P(1),P(2),P(3) 

00  INPUT  "HEIGHT  VECTOR" ,V(1,1),V(1,2),V(1,3) 

DO  INPUT  "SDH-MAJOR  VECTOR", V(2,1),V(2,2),V(2,3) 

IOC  INPUT  "SEMI-MINOR  VECTOR", V( 3,1 ),V(3, 2), V(3,3): 

INPUT  "RATIO", RO 

110  PRINT  :PRINT  "INPUT" :GOSUB  460 


120 

130 

140 

150 

160 

170 

mo 

190 
200 
21 0 
220 
230 
240 
250 
2o0 
270 
200 
290 

one 


FOR  J-1TU  3:  FOR  K-1T0  3sU(K)-V(J,IC) tNEXT  K 

GOSUf;  470:M(J)-M1:F0R  K-1T0  3:H(J,K)»U(K):HEXT  KrftFXT  J 

1)-U(2,1 ) *11(3,1  )+H(2,2)*l!(3,2)+U(2.3|  *W(3,3) : 

IF  AP,S(D)[.00fi7  THEN  260:D*ARCC0S(D):PRINTtlSING  20,D 

J1-2:J2*3:IF  M(2)]M(3)  THEN  160: Jl-3: J2-2 

C»O:F0R  J-1T0  3: IF  V(Jl,J)-0  THEN  170:C-C+1 

NEXT  J: IF  C]1  THEN  100,K-J1:J1-J2:J2-K 

F-0:F0R  J-1T0  3:IF  F]A0S(V(J1,J))  THEN  200 

IF  V(J2,J)»0  THEN  200:F»ARS(V(J1,J)):K1*J 

NEXT  J 

!'(J1 ,K1 }-0:F-U(Jl ,1 )+U(Jl ,2)+H(Jl ,3) 

U(Jl,Kl)-(U(2,l|*U  3,1)4«(2,2)*W(3,2)+H(2,3)*H(3,3))/F 

FOR  J-1T0  3:U(J)-W  Jl.JhflEXT  J:G0SUB  470 

FOR  J-1T0  3:V(J1 ,J)«U(J)*M(J1 ) :W(J1 ,J)«ll(J):NEXT  J 

PRINT  "SOLUTION":UOSUB  4C0 

INPUT  "INSIDE  TEC  (YES-1 ,MO-O)"tl : IF  1-0  THEN  70 

IliPUT  "THICKNESSES  ,T(1),T(2),T(3) 

PRINT  "THICKNESS  ARE":PRIMTUSING  30,T(1);T(2);T(3) 
T4-T(3): IF  M(3)/M(2)].n  THEN  300: 

T4-T4/1.01G447*(?I(3)/1 1(2))  1.071034) 
cn)-!'(2,2)*U(0,3)-U(2,3)*l/(3,?j: 


0 ? «N  3,1  *11  2,3  -III 


*»  i 

t.,i 

*>  r* 


*H(3,3): 


."’1(2,1  *11(3,2  -N  2.0*11(3.1) 

C1-ARS(0(1  )*H(1 .1  +0(2)*U(1  »2)'*0(3)*tJ(l,3)) 
01»T(1)/C1:02-T(2)/C1:I!-C1*M(1):N-1 


310 
320 

330  FOR  K-1T0  2:M4-M(2)/RO:;KH1(3)/RO 


350 

300 

370 

300 

:;oo 


FOP.  J«1T0  3::i4-M4+N*H(2,J)*V(l,J):M5-M5««Wl(3,J)*V(l,J): 
NEXT  0 

h»M(  2} -.  F»T(  3)  *SQR(  E ! 2+H ! 2)/l  UG-E/II 
Al(K)-l(2)-F-T(l)*G:A2(i;)-M4-F+T(2)*G 
E»;!(3)-nri:Fl»T4*SQR(E!r+H!2)/H:Cl-E/ll 


n(K)«iI(3)-Fl-T(l)*fil:C2(K)-MG-Fl+T(2)*Gl:N— 1:I!EXT  K 
A1».5*(AK1)+A1(2)):A2«.5*(A2(1)+A2(2] 
:r.2».n*(R2(l)*n2(2)) 


2)):R1-.5*(31(1)*B1(2)) 
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Table  A-XXII.  TECIN  Program  Listing  (Continued) 

400  R2-A1  /A2 : P.3-P1  /H2:R1  ■ ( R2*D2+R3*A2) / (A2+R2) . 

410  FOR  J-1T0  3:P1 (J)*Pfd)+W(l ,J)*D1+.5*((A1 (1 l-Al (2) )*W(2»J)+ 
(D10)-D1(2))*U(3,J)J:V1  l,J)-V(l,J)-M(l,J)*(Dl+02) 

420  VI (2,J)«A1*W(2,J) :V1 (3,J  ‘ “ ‘ 

430  PRINT  fcIHSIDC  TEC" 


M*W(3,J):tOT  J 


440  PRIHTUSIMG  30,P1(l);Pl(2);Pl(3hV1(l,1hV1(l ,2hVl(l,3), 
VI  (2,1 ) ;V1  (2,2)  ;V1  (2,3)  -,V1  (3,1 ) ;V1  (3,2)  ;V1  (3,3)  ,R1 
450  pniirrusinr,  4o,r2,R3:prhitusiiig  50^,52:0^0  70 
460  PRINTUSIHG  30,P(1 ) jF(2) gP(3) ;V(1  ,f) gV(l ,2) jV(l , 

V{2,1 ) ;V(2.2) ;V(2,3) ;V  3,1 ) ;V(3,2) ;V(3,3) .ROtRETflRR 

470  ! 11  -Sf»r.  (H  (1 ) ! 2-H ’ 2) ! 2+U  ( 3)  1 2 j 

4C0  ll(l)»U(l)/i!l:»l(2)-H(2)/Hl:U(3H'(3)/in 

4PO  KETllP..fl  :C!1D 
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Table  A-XX1II.  PARS  Prograa  Listing 


1U  Dili  A(2C)tB(3),C(2r.),P(3),A2(2G),n2(25)tC2{25)iR(3t3) 

;’o 

30  % w m # -mw.mi1  -mw.tm  -mw.tm 

/n 

50  iilHPUT  FOR  FACE  H 
60  XAZIiUTN  » -*#.«  ELEVATIOfi  ■ 

70  rsponiT » -###.«##  -tm.mp 

go  aiane  e»».  -f.fww  x+-< ?.####  Y+-#.####  z ■ immmw? 

.?* 

P0  R1-1F.5:  SELECT  D :SELECT  PRINT  2ir.(30):C3-.WnO:Cf)».nnoO'i 

100  INPUT  "NO.  OF  FACES  (N[26)",F 

110  FOR  0*1 TO  FifiOSCR  340: NEXT  J 

120  PRINT  " PLAilES  POSSIBLE  SOLUTION" 

130  F2*F-2:F1*F-1:N-0:F0R  J-1T0  F2:J1-J+1:F0R  K-J1T0  FI 

140  D-A2(J)*A2(N)+C2(J)*B2(K)+C2(J)*C2(K):IF  ABS(n)]C3  THEN  240 

ISO  Kl-K+l :F0R  L-klTO  F 

100  D*A2(J)*A2(L)+B2(J)*B2(L)+C2(J)*C2(L):IF  ABS(D)]C0  THEM  230 
170  D»A2(K)*A2(L)+B2(k)*B2  L)+C2(i:)*C2iL):IF  Ar.S(0)]C8  THEN  230 


ion  n 
lr.j  n 
:oj  r 


1.1) -A(Jj: 
2.1  -A  K): 

3.1) »A(L 


:R(1.2 
R 
:R 


« <» 
•»  <i 


-3(J):R(1,3)-C(J 
■B  K :R  2,3  -C  K 
«C(L):R(3t3)-C(L] 


210  OOSUr.  230:II*1H-1  :D-SQP.(A1 12+B1  I2+C1 12) 

220  IF  0]R1  THEM  230:Al-Al+SCN(Al)*CP:Dl-r»l+?QH(Bl)*CP: 

ci-cnsoii(ci )*C9:PRinTusnr.  3o,j,k,l ,ai  ,r<i  ,ci 
230  NEXT  L 
240  NEXT  KrUEXT  J 
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•0  THEN  120 


250  LIPtfT  "NEW  ARB  (YES-1  )",J: IF  JHO  THEN  1 
20!)  INPUT  "FACE  IMSIER  (NO  H0RE-O)h,J:IF  J-«' 

270  GOSUC  340: GOTO  2C0 
230  Dl-r:(l,ll*R(2,2}*R(3,3)+R(3,l)*R(1.2)*R(2,3)+R(2,l)*R(3,?)* 

R(la3)-n>vlm2v2H(lv3}-R(2tlHO»2M3»3H(l»l)* 

R(3,2)*n(2,3) 

223  ir  01 PO  Tip  300:01-. 00001 

303  A1-(3{1,2)*(R(2,3)-R(3,3))+R(2,2)*(R(3,3)-R(1,3))+R(3,2)* 

(K  1,3  -R(2,3)  )/01 

310  ill- (P.  1,3  *(n(2,l  -R(3,1))+R(2,3)*(R(3,1  )-R(l ,1  ))+R(3,3)* 
(Rl.l  -11(2,1)  )/Dl 

320  Cl-lRi'1,1  *(R(2,2.  -R(3,2))+R(2,1)*(R(3,2)-R(1,2))+R(3,1)* 
M l,2)-3(2,2)))/Dl 
330  RETURi* 

V.O  INPUT  "A.E-PT-O,  3-PTS-l,  E0-2",0:PRINT  : PR  INTI’S  I NO  DO.J: 

IF  0[]1  THEN  400 

350  SELECT  PRIIIT  00S(C4):F0R  K-1T0  3:PRINTHSINC  4%K 
3uO  iriPMT  "X(!0",R(K.l): INPUT  "Y(Kr,R(K,2):iriPUT  "ZW.RtK.S) 
370  FOR  L-1TU  3: IF  R(K,L)[]0  THEN  380:R(K,L)-.U0001 
3!iO  NEXT  l : NEXT  K:GOSW  200:A(J)-A1:P(J)-B1:C(J)-C1: 

SELECT  PT.Iinr  215(00) 
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Table  A-XXXIX.  PARS  Prograa  Listing  (Continued) 


390  PRINT  "3-P0IHTS"iPRINTUSING  20,R(l,1)sR(1,2);R(l,3);R(2,1h 
R(2,2) jR(2,3):PRINTUSIHG  20fR(3J);Ri3.2);R(3.3):60T0  490 
400  IF  0[]2  THEN  410:INPIIT  "A\A(J):INPUT  SR\B(J): 

INPUT  "C",C(J)i INPUT  "D".D:60T0  470 
410  INPUT  "AZIMUTH  (ROTATION) ",G: INPUT  "ELEVATION  (FALLBACK) "aH 
420  INPUT  "P(X)".P(l)s INPUT  *P(Y)\P(2):  INPUT  “P(Z)",P(3) 

430  FOR  K-1T0  3: IF  P(K)[]0  THEN  440:P(K)«.00001 

440  NEXT  K:A(J)*CO$(0)*COS(H)sB(J)*SIN(G)*COS(H):C(J)»SIM(ll) 

450  D-A(J)*P(1)+B(J)*P(2)+C(J)*P(3) 

400  PRIimiSINP.  00»G,H:PRIIITHSINC  70,P(1),P(2),P(3) 

470  IF  n[]0  THEN  400:D«.00001 
4T.0  A(J)»A(»1)/D:B(0)»B(0)/D:C(U)»C(J)/D 
490  r,-1/sqR(A(J)l2+B(J)!24C(J)!2):A2(J)-A(J)*r.:B2(JH(i1)*Cj 
C2(J)*C(J)*fi 

500  PRINTUSING  30aA2(J)aB2(U)aC2(J)aG 
51')  RETURN  :END 
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